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EXPERIMENTAL  AND  THEORETICAL  STUDIES  OF  HYDROGENATED 
AMORPHOUS  SEMICONDUCTOR  ALLOYS  AND  SUPERLATTICES 

A.  STATEMENT  OF  THE  PROBLEMS  STUDIED 

The  electronic  properties,  band  structure,  and  the  distribution  of  localized  states  in 
semiconductors,  thin-film  hydrogenated  amorphous  Si  and  its  alloys,  and  disordered 
systems. 

B.  SUMMARY  OF  THE  MOST  IMPORTANT  RESULTS 

Study  of  Bulk  and  Surface  Recombination  in  a-Si:H 

Using  our  combined  photoconductivity  and  diffusion  length  measurement 
technique  we  have  derived  bulk  mobility-lifetime  products  and  the  surface  recombination 
velocity  of  both  types  of  carriers  in  undoped  hydrogenated  amorphous  silicon  (a-Si:H) 
deposited  by  glow  discharge  decomposition  of  silane.  The  study  was  based  on 
comparison  of  the  experimental  data  with  the  classical  theory  of  semiconductors.  The 
measurements  were  carried  out  as  a  function  of  light  intensity  with  two  types  of 
illumination;  a  He-Ne  (7328  A)  laser,  and  a  He-Cd  (4420  A)  laser.  The  results,  at  a  light 
intensity  of  25  mW/cm^,  have  shown  that  for  the  same  photon  flux  the  He-Ne  laser  yields 
a  photoconductivity  15  times  larger  than  the  one  achieved  with  the  He-Cd  laser.  On  the 
other  hand,  the  corresponding  ratio  for  the  minority  carrier  mobility-lifetime  product  was 
lower  only  by  a  factor  of  two.  These  results  were  interpreted  assuming  surface 
recombination  velocities  of  the  order  of  10-100  cm/sec  and  a  moderate  improvement  in 
the  quality  of  the  material  with  film  thickness.  We  have  also  found  that  the  surface 
recombination  velocity  may  be  quite  different  for  the  two  surfaces  of  a  film,  indicating  that 
growth  conditions  play  an  important  role  in  the  elearonic  structure  of  the  surface.  From 
the  illumination  intensity  dependence  we  were  able  to  derive  information  regarding  the 
states  distribution  on  these  surfaces. 

Measurements  of  the  above  photoelectronic  properties  as  a  function  of 
temperature  has  revealed  information  on  the  temperature  shift  of  the  electrons'  and  holes' 
quasi-Fermi  levels.  For  example,  we  have  found  that  while  the  photoconduaivity  is 
almost  insensitive  to  temperature  variations  around  room  temperature,  the  mobility-life 
time  product  of  the  minority  carriers  has  an  activation  energy  of  E^^.  1 1  eV.  Under  some 
assumptions  on  the  material  parameters,  these  results  may  yield  the  density'  of  state?  in 
some  energy  imervals  of  the  pseudogap  of  the  materials. 
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Within  the  framwork  of  the  present  study  we  have  developed  the  theory  of  the 
photocanier  grating  (PCG)  technique  and  utilized  it  to  interpret  our  experimental  results 
on  a-SiiH.  This  has  enabled  the  derivation  of  the  mobility-lifetime  (pt)  products  of  both 
types  of  carriers  in  undoped  and  doped  a-Si;H.  Previous  experimental  works  using  this 
technique  were  limited  to  the  derivation  of  the  px  products  of  the  minority  carriers  only 
Furthennore,  we  were  able  to  develop  a  criterion  under  which  the  application  of  the  PCG 
technique  yields  indeed  the  value  of  the  minority  carriers  tix  product.  We  also  showed  that 
this  criterion  is  not  always  met.  Hence,  we  made  the  PCG  a  routine  tool  for  the  evaluation 
of  the  microscopic  photoelectronic  properties  of  a-Si;H. 

In  our  studies  we  have  also  used  the  photocarrier  grating  technique  for  the  first 
time  in  the  high  electric-field  regime.  The  results  were  shown  to  confirm  the  so  far 
unproven  theoretical  predictions  for  this  regime.  We  have  demonstrated  that  by 
application  of  the  PCG  technique  in  the  high  field  regime,  accurate  values  for  the  ratio  of 
the  two  carriers  |ix  products  can  be  deduced.  This  is  in  contrast  with  the  fact  that  their 
sum  cannot  be  derived  accurately  fi'om  the  measurement  of  photoconductivity  because  one 
cannot  determine  accurately  the  carrier  generation  rate.  Combining  the  ratio  determined  by 
the  high-field  PCG,  with  the  low  field  data,  yields  then  accurate  values  for  the  mobility  - 
lifetime  produrts  of  both  carriers. 

In  the  low  electric-field  regime  the  many  reports  concerning  the  utilization  of  the 
PCG  technique  have  assumed  that  ambipolar  transport  takes  place  in  the  PCG  when  it  is 
applied  to  a-Si;H.  This  assumption,  which  is  decisive  in  the  interpretation  of  the 
experimental  results  in  terms  of  the  ambipolar  diffusion  length,  has  not  been  tested  so  far 
In  our  work  we  developed  the  first  corresponding  ambipolarity  criterion  and  have 
demonstrated  that  while  ambipolarity  is  maintained  in  device  quality  a-Si:H  it  does  not 
apply  to  lightly  doped  or  poor  quality  a-Si:H  materials.  The  results  were  shown  to  confirm 
our  theoretical  suggestion  that  the  observed  ambipolarity  in  a-Si:H  is  due  to-trapping 
effects. 

The  relation  between  the  inhomogeneity  along  the  direction  of  the  film  growth  and 
the  surface  recombination  velocity  has  been  studied  in  device  quality  a-Si:H  films.  It  is 
found  that  the  material  is  electronically  inhomogeneous  "deep"  into  the  bulk  up  to  a 
typical  depth  of  the  order  of  1000  A.  Hence,  the  surface  recombination  velocity  in  a-Si:H 
is  not  a  true  surface  property  as  in  crystalline  semiconductors.  This  result  is  important  for 
modeling  electronic  devices,  since  these  arc  typically  only  a  few  thousand  angstroms  thick 

Using  the  deposition  temperature  as  a  control  parameter  for  variation  of  the 
dangling  bond  concentration  and  the  relative  position  of  the  dark  Fermi  level,  the 
photoconductivity  and  the  ambipolar  diffusion  length  were  measured  on  samples  prepared 
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by  sputtering  and  glow  discharge  at  different  temperatures.  This  has  enabled  us  to 
differentiate  between  different  suggested  models  for  the  energetic  position  of  the 
recombination  centers,  as  weU  as  the  recombination  kinetics  in  a-Si;H.  It  was  found  that 
the  "defect-pool”  model  yields  the  best  description  of  the  system.. 

Finally,  we  were  able  to  achieve  the  first  combined  application  of  the  Metal  Oxide 
Semiconductor  and  the  Photocarrier  Grating  configurations.  This  combination  enabled  us 
the  first  simultaneous  study  of  the  two  carriers  mobility-lifetime  products  and  their  light 
intensity  exponents  as  a  fimction  of  the  position  of  the  Fermi  level,  in  undoped  a-Si;H 
We  found  that  anticorrelations  and  correlations  prevail  between  these  two  sets  of 
quantities.  The  conclusion  we  derive  from  these  behaviors  is  also  that  the  defect  pool 
model  accounts  for  the  phototransport  data  much  better  than  any  other  model. 

We  have  studied  the  Raman  spectrum  in  a-Si  :H  as  a  function  of  film  thickness  in 
order  to  follow  variations  in  the  degree  of  order  in  the  amorphous  network.  The  results 
indicate  an  improvement  of  the  order  with  sample  thickness.  There  are  considerable 
differences  between  the  optical  phonon  modes  at  the  surface  and  in  the  bulk  which  we 
attribute  to  the  fact  that  the  surface  layer  is  less  ordered  than  the  bulk.  These  results  are  in 
agreement  with  our  former  observation  of  improvement  of  the  photoelectronic  properties 
with  the  film  thickness.  More  recently  we  have  studied  a-SixCj.x  H  alloys  and  found  quite 
surprisingly  that  the  silicon  network  order  increases  with  carbon  content  Furthermore, 
due  to  the  possibility  of  bandgap  tuning  by  carbon  content  and/or  variation  of  laser 
wavelength  excitation,  we  found  a  resotumt  Raman  effect. 

The  Semiconductor/Electrolyte  Interface 

Non-Parabolicity  in  the  Lowest  Conduction  Band  of  CdS. 

Pulse  measurements  on  the  semicondurtor/electrolyte  system  were  used  to  induce 
and  study  space-charge  layers  at  the  semiconductor  surface.  Measurements  on  CdS  show 
that  the  free-electron-like  lowest  conduction  band  in  hexagonal  CdS  is  highly  non¬ 
parabolic.  A  simple  model  for  the  structure  of  this  band,  which  assumes  parabolicity  up  to 
0. 125  eV  above  the  band  edge  and  a  linear  dependence  of  the  energy  on  the  wavevector  at 
higher  energies,  accounts  well  for  the  experimental  results.  The  density-of-statcs  function 
that  emerges  fi’om  these  results  is  compatible  with  reported  theoretical  calculations. 

In  addition,  the  measurements  indicate  that  surface  states  are  practically  absent  at 
the  CdS  surface  in  contact  with  the  electrolyte.  They  also  shed  light  on  the  process  of 
charge  leakage  acros:  the  CdS/electrolyte  interface,  leakage  that  occurs  mostly  when 
strong  accumulation  layers  are  induced. 


Surface  States  at  the  Silicon/Electrolyte  Interface. 

Pulse  measurements  on  the  Si^lectrolyte  interface  were  used  to  study  surface 
states  at  the  silicon  surface.  We  find  that  for  CP-4  etched  silicon,  surface  states  exist  in 
the  upper  half  of  the  energy  gap,  with  a  total  density  of  about  10^2  cm"2.  Addition  of  a 
minute  amount  of  hydrofluoric  acid  to  the  electrolyte  reduces  their  density  by  one  order  of 
magnitude.  It  has  been  reported  that  the  surface  recombination  velocity  on  an  oxidized 
surface  is  considerably  lowered  by  an  HF  treatment,  suggesting  a  corresponding  reduction 
in  the  density  of  surface  recombination  centers.  Our  results,  however,  constitute  the  first 
direct  evidence  for  such  'a  reduction.  It  appears  that  Si-H  bonds  are  formed  at  the 
interface,  just  as  has  been  reported  for  an  HF-treated  free  Si  surface.  Very  likely,  the 
formation  of  such  bonds  leaves,  on  both  types  of  surface,  considerably  fewer  dangling 
bonds  and  hence  considerably  fewer  surface  states. 

Study  of  Density  of  States  in  a-Si:H. 

Pulsed  measurements  on  the  solid-electrolyte  system,  which  proved  very  useful  in  the 
study  of  crystalline  semiconductors,  have  been  found  to  be  equally  effective  when  applied 
to  hydrogenated  amorphous  Si  films.  Here,  as  well,  the  a-Si;H/electrolyte  interface  is 
essentially  blocking  to  current  flow  and,  as  a  result,  surface  space-charge  layers,  ranging 
from  large  depletion  to  very  strong  accumulation  conditions,  have  been  induced  and 
studied.  Measurements  in  the  depletion  range  under  illumination  yielded  directly  the  total 
density  of  occupied  states  in  the  entire  energy  gap.  This  is  very  useful  in  obtaining  a  quick 
and  reliable  assessment  of  the  quality  of  the  amorphous  films.  In  high-grade  films  we  find 
that  the  total  density  of  occupied  states  is  around  cm*^  The  data  in  the  accumulation 
range,  on  the  other  hand,  provide  useful  information  on  unoccupied  states  near  the 
conduction  band  edge.  The  blocking  nature  of  the  amorphous  Si/electrolyte  interface  was 
utilized  also  to  apply  a  sweep-out  technique  for  an  accurate  determination  of  px,  the 
product  of  the  electron  mobility  and  lifetime,  even  when  this  value  is  very  low.  In  a  rather 
poor-quality  film,  for  example,  we  find  px  to  be  5x10'* 

Cermets  and  Metallic  Colloidal  Suspensions 

Tunneling  and  Percolation  Behaviour  in  Granular  Metals. 

The  nature  of  the  percolation  process  in  granular  metals  was  examined  for  the  first 
time  by  a  computer  simulation  of  a  system  of  metallic  grains  embedded  in  an  insulating 
matnx.  Assuming  that  the  intergrain  conduction  is  due  to  quantum  mechanical  tunneling  it 
was  found  that  a  percolation-like  critical  behavior  of  the  conductivity  is  obtained,  but  that 
a  percolation  universal  behavior  will  be  found  only  in  a  very  special  case.  In  contrast,  the 


behavior  of  the  electrical  noise  does  not  deviate  substantially  fi-om  the  universal  one. 
Comparison  of  these  results  with  experimental  observations  suggests  that  in  the  metallic 
range,  both  transport  properties  are  controlled  by  the  continuous  metallic  network  rather 
than  by  intergrain  tunneling.  We  propose  that  the  metallic  network  resembles  the 
previously  studied  system  of  "inverted  random  voids". 

Local  Fields  Around  Clusters  of  Prolate  Spheroids. 

T_-matrix  formalism  was  used  to  calculate  local  electric  fields  around  clusters  of 
prolate  spheroids  in  the  long  wavelength  regime.  The  calculations  are  performed  as  a 
function  of  interparticle  distance  as  well  as  the  angle  of  orientation.  The  observed  red 
shifts  in  the  resonant  wavelengths  of  the  charaaeristic  peaks  are  shown  to  obey  an 
exponential  relationship  as  a  function  of  interparticle  separation  and  a  sinusoidal 
relationship  as  a  function  of  angle  of  rotation  of  the  spheroid.  The  behavior  of  the  cluster 
was  examined  and  the  two  effects,  of  separation  and  rotation,  were  compared. 

The  results  of  these  calculations  were  applied  to  cermet  materials  and  metallic 
colloidal  suspensions  in  alkali  halide  crystals  and  their  effect  on  enhanced  Raman 
scattering  was  discussed. 
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NON-PARABOLICITY  IN  THE  LOWEST  CONDUCTION 

BAND  OF  CdS 
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(Xeetived  36  Ftbnary  1990;  aeetpud  IS  tiareh  1990) 

Abctract — Pulee  measurements  on  the  CdS/electrolyte  system  ate  used  to  induce  and  study  spitce<harge 
layers  at  the  CdS  surface.  Such  measurements  show  that  the  free-dectron'like  lowest  conduction  band 
in  hcxaional  CdS  is  highly  non-parabolic.  A  simple  model  for  the  structure  of  this  band,  which  assumes 
parabolldty  up  to  0. 12S  eV  abow  the  band  edp  and  a  linear  dependence  of  the  energy  on  the  wavevector 
at  higher  energiet,  aoeounu  well  for  the  experimental  reeulu.  TIm  density  of  sums  function  that  emerges 
from  these  results  is  compatible  with  reported  theoretical  calculations. 

In  addition,  the  measurements  indkatt  that  surface  sutes  are  practically  absent  at  the  CdS  surface  in 
contact  with  the  electrolyte.  They  also  shed  light  on  the  process  of  charge  leakage  across  the 
CdS/electrolyte  interface,  leakage  t^t  occun  mostly  when  strong  accumulation  layers  are  induced. 


Keywords:  CdS,  band  structure,  surface,  quantum  wells. 


INTRODUCTION 

Band  structure  calculations  for  hexagonal  CdS  [1.2] 
indicate  the  presence  of  a  nearly-spherical,  free- 
electron-like  lowest  conduction  band  extending  up  to 
about  1.8  eV  above  the  band  edge.  While  these 
calculations  provide  a  broad,  overall  picture  of  the 
band  structure,  little  detail  can  be  gleaned  about  the 
precise  shape  of  the  lowest  conduction  band.  Experi¬ 
mental  results  such  as  those  derived  from  optical 
measurements  are  also  lacking.  In  this  paper  we  make 
use  of  the  semiconductor/electrolyte  (S/E)  system  to 
study  the  surface  space<harge  layer  of  CdS.  The 
higher  reaches  of  the  lowest  conduction  band  become 
occupied  in  very  strong  accumulation  iayen.  Such 
layers  can  be  induced  in  the  S/E  system  provided 
pulsed  [3-5]  rather  than  d.c.  biases  are  employed. 
Electrochemical  processes  are  then  largely  circum¬ 
vented,  permitting  the  derivation  of  the  physical 
characteristics  of  the  surface  with  minimal  involve¬ 
ment  of  the  ratht-i  complex  surface  chemistry. 
Similarly  to  the  case  of  ZnO  [6, 7]  and  Si  [5],  accumu¬ 
lation  layers  of  surface  electron  densities  as  high  as 
I0'*cm~’  can  be  induced  and  studied  at  the  CdS 
surface.  Analysis  of  the  data  obtained  for  this  range 
indicates  that  the  lowest  conduction  band  is  highly 
non-parabolic,  and  provides  a  fairly  good  idea  as  to 
the  shape  of  the  band  up  to  about  1.5  eV  above  the 
conduction-band  edge. 

The  S/E  interface  is  not  perfectly  blocking  and 
some  charge  leaks  between  the  two  phases,  especially 
when  strong  accumulation  layers  are  induced  at  the 
CdS  surface.  The  pulse  measuremenu  yield  a  detailed 
characterization  of  the  leakage  process.  They  also 
indicate  that  surface  states  are  practically  absent  at 


the  CdS  surface  when  in  contact  with  the  electrolytes 
used. 

THE  SPACE-CHARGE  LAYER  AND  THE 
SEMICONDUCTOR/ELECTROLYTE  INTERFACE 

In  this  section  we  review  briefly  the  characteristics 
of  the  semiconductor  surface  si»ce-charge  layer  [8] 
and  consider  several  features  of  the  semiconductor/ 
electrolyte  system,  stressing  those  points  that  are 
directly  related  to  our  present  work.  We  consider  a 
homogeneous,  non-degenerate,  n-type  semiconduc¬ 
tor  having  fully  ionized  donors.  For  depletion  and 
accumulation  layers  (one-carrier  system),  the  surface 
space-charge  density  (per  unit  area)  can  be  expressed 
as 

(1) 

where  q  is  the  absolute  magnitude  of  the  electronic 
charge  and  N,  is  the  surface  electron  density.  By 
definition,  N,  is  zero  under  flat-band  conditions, 
positive  for  accumulation  layers  and  negative  for 
depletion  layers.  The  space-charge  capacitance  is 
given  by 

C«-|G,/F.|.  (2) 

where  V,  is  the  potential  barrier  height  at  the  surface. 
The  effective  charge  distance  L,,  which  represents  the 
effective  width  of  the  space-charge  layer,  is  given  by 
the  expression 

C,-ia,/L,.  (3) 

where  k  it  the  relative  dielectric  constant  of  the 
temiconductor  and  Co  is  the  permittivity  of  free  space. 
First  integration  of  Poisson's  equation  leads  to  m 

(4) 
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where  (he  positive  sign  applies  to  accumulation 
layers,  the  negative  sign  to  depletion  layers;  is  the 
electron  bulk  concentration,  kg  is  Boltzmann's 
constant  and  T  the  absolute  temperature. 

In  the  general  case,  the  function  F,  is  given  by  (3] 


1  ■+■  exp(£/k,  T  +  Wt-  0,) 
1  +exp(Elk,T  +  w*) 


(5) 


where  N(E)  is  the  density  of  states  in  the  conduction 
band,  **>»  =  W„lktT,  IF,  being  the  energy  separation 
between  the  Fermi  level  and  the  conduction>band 
edge  in  the  bulk  (positive  for  the  n-type  sample 
being  considered),  and  v,-qt',lkfT  is  the  sosalled 
dimensionless  barrier  height. 

For  parabolic  conduction  bands  the  function  F, 
assumes  the  simpler  form 


F,-y2|i 


(4/3)n"''’e*p(H'jj' 

•jin 

(1 +exp(x +  H>4-r,)l-' dx -r,- IV  . 


(6) 


For  strong  accumulation  layers,  however,  in  which 
the  conduction-band  edge  lies  well  below  the  Fermi 
level,  taking  the  conduction  band  u  parabolic  may 
not  always  be  a  good  approximation.  This  is  defin¬ 
itely  the  case  in  CdS  and  one  should  use  eqn  (5), 
substituting  for  N{E)  the  actual  density  of  states  in 
the  (non-parabolic)  conduction  band. 

Another  imporunt  factor  that  needs  to  be  consid¬ 
ered  in  strong  accumulation  layers  is  associated  with 
the  spatial  confinement  of  the  electron  gas.  When  the 
effective  charge  distance  L,  becomes  comparable  to  or 
less  than  the  de  Broglie  wavelength  of  a  conduction 
electron,  quantization  effects  set  in.  The  classical 
calculations  reviewed  above  may  no  longer  be  ade¬ 
quate,  the  less  so  the  stronger  the  accumulation  layer 
(larger  N,,  smaller  L,).  One  should  then  solve  sclf- 
consistently  Poisson's  and  Schrodinger’s  equations 
[6.7,9,10].  Unfortunately,  such  calculations  have 
not  been  carried  out  for  CdS.  However,  the  similarity 
of  the  values  of  the  dielectric  constant  and  the 
effective  mass  (at  the  bottom  of  the  band)  in  CdS 
{k  -9.02  and  m/m«“  0.205.  where  m,  is  the  free- 
electron  mass)  and  ZnO  (x  »  8.5.  m/m*  «  0.25)  per¬ 
mits  the  use  of  the  ZnO  calculations  (6,  7]  for  CdS 
without  introducing  too  large  an  error. 

For  non-degenerate  surface  conditions,  corre¬ 
sponding  to  depletion  and  vreak  accumulation  layers, 
only  the  bottom  of  the  band  is  involved  and  this  can 
safely  be  uken  as  parabolic.  Equation  (6)  then 
reduces,  to  a  good  approximation,  to  the  simpler 
expression 


(7) 


The  semiconductor/electrolyte  system  considered 
here  consisu  of  a  semiconductor  surface  in  contact 
with  an  indifferent  electrolyte  such  as  an  aqueous 
solution  of  a  salt  [KCl  or  CafNO,)}].  By  biasing  the 
electrolyte  with  respect  to  the  semiconductor,  charge 
can  be  induced  electrostatically  at  the  semiconductor 
surface.  This  is  made  possible  by  the  essentially 
blocking  nature  of  the  S/E  interface,  which  largely 
prevents  charge  transfer  between  the  two  phases. 
Such  blocking  is  due  to  the  fact  that  in  the  electrolyte 
the  negative  ions  (such  as  Cl~)  usually  lie  well  below, 
and  the  positive  ions  (such  as  K.*)  well  above,  the 
semiconductor  band  e>''£es  [11].  Hence,  in  an  n-type 
semiconductor  being  considered,  for  a  negatively 
biased  electrolyte,  the  negative  ions  pile  up  against 
the  interface  without  being  able  to  transfer  electrons 
into  the  conduction  band,  and  a  depletion  layer  is 
induced.  In  the  reverse  polarity,  positive  ions  pile  up 
against  the  interface  and  an  accumulation  layer  is 
induced.  The  induced  charge  is  in  general  distributed 
between  the  surface  space-charge  region  and  surface 
states  (if  present).  In  the  absence  of  a  foreign  layer 
(such  as  an  oxide)  at  the  semiconductor  interface,  the 
applied  bias  drops  across  three  adjacent  regions 
[8, 12];  the  space-charge  layer  at  the  semiconductor 
surface,  and  the  Helmholtz  and  Gouy  layers  in  the 
electrolyte.  Of  these  three  layers,  the  latter  two  are 
much  thinner  than  the  first,  at  least  under  pulsed  bias 
conditions,  so  that  practically  the  entire  applied  bias 
drops  across  the  semiconductor  space-charge  layer. 
This  is  evidenced  by  space-charge  capacitance 
measurements  on  metal/electrolyte  interfaces  [3]. 
Hence  the  measured  bias  represents  to  a  very  good 
approximation  the  change  6y,  in  barrier  height  V, 
across  the  semiconductor  space-charge  Uyer. 

In  practice,  the  interface  is  not  perfectly  blocking 
and  to  a  lesser  or  greater  extent  some  charge  leaks 
across  the  S/E  interface.  The  leakage  charge  is  usually 
insignificant  when  depletion  layers  are  induced 
(electrolyte  biased  negatively  with  respect  to  the 
n-type  semiconductor),  it  becomes  appreciable,  how¬ 
ever.  under  accumulation  conditions,  the  more  so  the 
stronger  the  induced  accumulation  layer.  It  is  this 
leakage  that  prevents  the  formation  of  strong  accu¬ 
mulation  layers  by  d.c.  biases.  By  employing  our  pulse 
techniques  [3-5].  on  the  other  hand.. it  is  possible  to 
induce  extremely  strong  accumulation  layers  not  only 
on  CdS  but  also  on  a  number  of  other  semiconduc¬ 
tors.  such  as  Ge  [3],  Si  [5]  and  ZnO  (6.  7],  without 
being  hampered  by  charge  leakage.  Moreover,  these 
techniques  enable  a  straightforward  separation  of  (he 
different  components  of  the  induced  charge,  so  as  to 
obtain  the  variation  of  the  free  space-charge  density 
as  well  as  of  the  leaked  charge  and  of  the  occupancy 
of  the  interface  states  with  the  surface  potential 
barrier.  More  important  still,  the  use  of  a  limited 
number  of  short-duration  pulses  minimizes  electro¬ 
chemical  reactions  which  may  otherwise  greatly  com¬ 
plicate  the  interpretation  of  the  data.  By  the  same 
token,  surface  damage  that  usually  occurs  under 


F,’‘^2{exp(v,)-v,-\y'\ 
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prolonged  application  of  large  d.c.  biases  is  largely 
eliminated. 

EXPERIMtiVTAL 

The  samples  were  cut  from  single  crystal  (n-type) 
CdS  ingots  ranging  in  resistivity  between  about  0.1 
and  lOohm-cm.  The  samples'  dimensions  were  typi¬ 
cally  2  X  3  X  0.5  mm\  with  the  large  faces  perpen¬ 
dicular  to  the  (-axis.  Following  mechanical  polishing, 
each  sample  was  first  etched  in  HCI  and  then  in  a 
solution  of  K;Cr;07  in  dilute  H;S04,  at  90‘'C,  for 
10  min.  This  resulted  in  a  smooth  surface  with  a 
negli^‘ble  density  of  surface  states  and  low  leakage 
currents.  An  ohmic  contact  was  applied  to  one  of  tne 
large  surfaces.  The  wire  lead,  contact  area  and  the 
entire  sample  were  masked  by  epoxy  cement,  except 
for  a  small  area  ( ~  1  mm-)  to  be  exposed  to  the 
electrolyte.  The  exposed  area  was  usually  the  (OOOT) 
face  (the  “sulfur"  face).  The  CdS  sample  and  a  large 
platinum  electrode  were  immersed  in  a  CalNOj)}  or 


KCI  electrolyte.  These  are  represented  by  the  equiv¬ 
alent  circuiu  (to  be  discussed  b^ow)  between  points 
A  and  B  in  Fig.  1(a)  which  is  a  schematic  diagram  of 
the  experimental  arrangement.  The  pulse  generator 
(Hewlett-Packard  type  214A)  of  internal  resistance  r, 
(  =  500),  can  provide  positive  and  negative  pulses  of 
up  to  100  V  in  amplitude  and  0.05  to  200  #is  in 
duration.  It  is  connected  to  the  platinum  and  semi¬ 
conductor  electrodes  (points  A  and  B)  via  a  diode  D 
and  a  large  (0.1-1  uF)  series  capacitor  C.  The  diode 
polarity  is  such  as  to  permit  rapid  charging  of  the 
system  while  the  pulse  is  on  and  to  prevent  its 
discharge  through  the  low-impedance  pulse  generator 
after  the  termination  of  the  pulse.  An  electronic 
switch  S  can  short-circuit  the  Pt  electrode  to  ground 
at  any  prescribed  time  following  the  termination  of 
the  pulse.  The  residual  resistance  r,  of  the  shorted 
switch  is  a  few  ohms.  Short-duration  (0.1-10  #ts) 
pulses  are  applied  singly,  one  pulse  per  data  point 
taken,  in  order  to  minimize  any  electrochemical 


Fig.  I.  (a)  Circuitry  of  the  experimental  arrangement.  The  sample  is  rimulated  (between  points  A  and 
B)  by  the  space-charge  capacitance  C,  and  a  resistance  r  in  series,  and  by  the  leakage  and  surface-state 
equivalent  circuits,  (b)  Schematic  representation  of  expected  pulse  responses  of  the  voltage  drops  K.  and 
y,.  (c)  Actual  pulse  response  of  F.  (upper  trace.  I  Vcm'')  and  F,  (lower  trace,  lOOmVcm-')  for  a 

CdS/electrolyte  interface. 
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reaction  that  may  take  place  at  the  CdS  surface.  A 
Tektronix  type  7407  CRO  is  used  to  monitor  the 
voltage  drop  V,  between  the  Pt  electrode  (pr  int  A) 
and  the  ohmic  contact  of  the  sample  (point  B),  and 
the  drop  K,  across  the  series  capacitor  C.  Strktly 
speaking,  the  relevant  oscillogram  is  that  of  the 
voltage  drop  V,  between  a  reference  electrode  and  the 
sample,  rather  than  V,.  However,  the  space-charge 
capacitance  at  the  platinum  surface  is  orders  of 
magnitude  larger  than  that  at  the  semiconductor 
surface.  Hence,  only  a  negligible  fraction  of  the 
applied  bias  is  expected  to  drop  across  the 
Pt/electrolyte  interface,  so  that  the  two  traces  should 
yield  practically  identical  results.  This  was  indeed 
checked  to  be  the  case  by  the  occasional  use  [6. 7]  of 
a  reference  electrode.  The  magnitude  of  C  is  chosen 
so  as  to  be  large  compa^  to  the  semiconductor 
surface  capacitance,  typically  by  a  factor  of  10-100. 
The  voltage  developed  across  it  (y,)  is  thus  smaller 
(by  the  same  factor)  than  the  voluge  developed 
across  the  space-charge  layer  (f'^). 

In  most  cases  a  space-charge  layer  already  exists  at 
the  semiconductor  surface,  before  applying  any  bias. 
It  is  characterized  by  an  equilibrium  barrier  height 
and  a  space-charge  density  Q^.  If  surface  states  are 
present,  there  may  also  be  an  equilibrium  surface- 
state  charge  of  density  Q^.  The  voltage  pulse,  of 
amplitude  anywhere  between  0  and  100  V,  charges  up 
the  semiconductor  surface  region.  As  will  be  shown 
in  a  moment,  the  voluge  drop  V,,  measured  just  after 
the  termination  of  the  pulse,  represents  to  a  very  good 
approximation  the  change  if',  ■=  f',  -  f'jo  tn  barrier 
height  induced  by  the  applied  pulse.  The  voluge  drop 
y,  across  the  series  capacitor  C,  again  measured  just 
after  the  termination  of  the  pulse,  yields  the  overall 
charge  Cm  ( **  Cf',)  passing  through  the  systems.  The 
charge  density  Q^„  (per  unit  area)  is  obuined  by 
dividing  Cm  by  the  area  of  the  semiconductor  surface 
exposed  to  the  electrolyte. 

In  general.  Cm  nude  up  of  three  components; 

0m“i0»  +  i0»i  +  Ql.»  (8) 

where  C»  ~  Quo  •*  *be  change  in  space-charge 
density,  6Q„  ■■  Cn  ~  Cuo  the  change  in  surface- 
state  charge  density,  and  Cl  the  density  of  the 
charge  leaked  across  the  S/E  interface.  In  terms  of 
electron  densities,  this  can  be  rewritten  as 


N„<-‘SN.  +  dN„  +  N^, 


where  is  the  toul  electron  deruity  passing 
through  the  system,  SN,  -  A/,  -  A^,o  i*  the  change  in 
surface  electron  density,  6N„  ■  is  the  den¬ 

sity  change  in  surface-suu  occupation,  and  the 
density  of  electroru  that  have  leaked  acrou  the 
interface. 

We  shall  now  discuss  the  experimenUl  procedure 
used  in  determining  the  magnitudes  of  the  three 
components  in  eqn  (8)  or  (9).  This  is  best  done  in 


terms  of  the  equivalent  circuits  shown  in  Fig.  1(a). 
The  ideal  case  of  a  perfectly  blocking  interface  and  no 
surface  sutes  is  simulated  amply  by  a  series  combi¬ 
nation  of  the  space-charge  capacitance  C„  and  the 
combined  resistance  r  of  the  electrolyte  and  the 
underlying  semiconductor  bulk  below  the  surface. 
Typically,  r  is  about  500,  approximately  half  of 
which  represents  the  electrolyte  resistance  and  the 
other  half  the  sample  resisunce.  Charge  leakage 
across  the  interface,  a  prominent  process  in  strong 
accumulation  layers  of  CdS,  is  simulated  by  adding 
a  parallel  branch  across  C^,  consisting  of  a  diode  in 
series  with  a  parallel  combination  of  a  capacitor  C« 
and  resistor  /{„.  As  will  be  shown  below,  such  a 
circuit  is  prescribed  by  the  experimenUl  observation 
that  the  leakage  occurs  very  fast  and  remains  stored 
at  or  near  the  interface  for  a  long  time  (R„C„).  Once 
the  charging  pulse  is  terminated,  the  stored  charge 
does  not  flow  into  C„,  hence  the  diode.  If  deep 
surface  states  are  also  present,  another  equivalent 
circuit  is  called  for,  shown  added  by  the  dashed  lines 
in  Fig.  1(a).  The  circuit  consists  of  a  resistor  R„ 
shunted  by  a  diode,  both  in  series  with  a  capacitor 
C„.  On  the  application  of  a  positive  pulse,  such  sutes 
charge  up  very  quickly,  hence  the  shunting  diode. 
Their  discharge  following  the  termination  of  the 
pulse  is  usually  much  slower,  the  more  so  the  deeper 
they  are  located  energetically  and  the  lower  the 
temperature  [8]. 

It  should  be  noted  that  the  equivalent  circuits 
shown  in  Fig.  1(a)  serve  merely  as  a  means  of 
classifying  the  different  processes  at  the  S/E  interface 
on  the  one  hand,  and  as  aids  in  the  understanding  of 
the  method  of  measurement,  on  the  other  hand.  They 
do  not  represent  accurately  the  actual  conditions  at 
the  interface.  In  particular,  all  elements  of  the  equiv¬ 
alent  circuits  (except  r)  are  in  reality  not  passive  but 
voltage  dependent. 

Figure  1(b)  illustrates  the  time  variations  of  y,  and 
y,  during  and  following  the  application  of  a  positive 
pulse  of  amplitude  V^.  All  the  values  marked  in 
Fig.  1(b),  as  well  as  in  the  following  analysis,  corre¬ 
spond  to  unit  surface  area.  The  F,  scale  is  expanded 
by  a  factor  of  100  or  so  relative  to  the  y,  scale.  The 
initial  rise  in  V,  (at  (^0  +  ),  up  to  the  value 
F,r/(r  +  r,),  is  very  fast  (about  0.01  ^s)  and 
represents  the  ohmic  voltage  drop  across  r.  Subse¬ 
quently,  V,  rises  much  more  slowly  as  C^.  C„  and  C„ 
charge  up.  At  the  termination  of  the  pulse  (r  >  T),  F. 
drops  rapidly  to  the  level  6  F, ,  where  S  F,  is  the  voltage 
developed  across  i  e.  the  induced  change  in  the 
potential  barrier  height.  As  for  F,  it  rises  monotoni- 
cally  up  to  r  «  T,  to  the  level  IC,  where  Q„  is  the 
total  charge  density  induced  by  the  pulse,  consisting 
of  6Qk,  Ql  Md  shorting  switch  S  is 

activated  at  r  >■  r  -f  dT,  where  the  interval  dT is  very 
short  (usually  0.1  fit),  sufficient  to  permit  accurate 
readings  of  F.  and  F,  just  after  the  termination  of  the 
pulse.  At  this  point  F,  drops  abruptly  to  practically 
zero  (F,  being  very  small  compared  with  F,),  and 
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dinife  redis^bution  imoiig  C«.,  C«.  C.  nnd  C 
begins  to  take  place.  In  the  first  stage,  the  diarge 
stored  in  C,  and  its  equal  counterpart  in  C  discharge 
relatively  fast  throuih  r-hr,.  The  decay  constant 
associated  with  this  proce^  ••  approximately 
(r  +  r,)C«.  (since  C«.  is  small  compared  with  C)  and 
is  typically  several  microseconds.  In  the  absence  of 
surface  states  (as  is  very  nearly  the  case  for  QIS),  Kr 
decays  to  the  value  QJC,  Ql  being  the  charge 
remaining  in  C  after  the  fast  decay  process.  This 
charge  will  be  referred  to  as  a  “lots”  since,  in  a  sente, 
it  detracts  from  the  induced  free  space  charge.  There¬ 
after,  K,  remains  practically  constant  (solid  curve  in 
the  bottom  diagram  of  Fig.  1(b)].  Actually,  it  decays 
very  slowly  (O.I-l  s)  through  the  input  resistance  of 
the  CRO.  If,  on  the  other  hand,  surface  states  are 
present  and  there  it  no  leakage,  F,  decays  with  the 
surface-sute  time  constant  of  approximately  R„Ca 
(dashed  curve  in  the  bottom  diagram  of  Fig.  1(b)].  In 
this  process,  electrons  trapped  in  the  surface  sutes  by 
the  charging  pulse  are  thermally  re-emitted  into  the 
conduction  bend.  If  both  leakage  and  surface  states 
are  present,  V,  first  decays  rapidly  to  the  value 
(dQ„  +  Qt)/C.  Subsequently  it  decays  with  the  time 
constant  RuC„  to  the  level  of  Qi_/C  and  thereafter 
remains  practically  constant.  This  behavior  enables 
the  separate  determination  of  the  three  components 
in  eqn  (8)  that  make  up  the  total  induced  charge 
density  .  As  pointed  out  above,  the  surface-sute 
component  is  essentially  absent  in  CdS. 

Typical  oscillograms  illustrating  the  actual  pulse 
response  of  the  CdS/electroiyte  interface  for  an  ap¬ 
plied  positive  pulse  (induced  accumulation  layer)  are 
displayed  in  Fig.  1(c).  The  upper  and  lower  traces 
depict  the  time  variation  of  K.  and  respectively. 
They  are  seen  to  be  very  similar  to  the  corresponding 
schematic  diagrams  in  Fig.  1(b),  indicating  that  the 
models  discussed  above  do  indeed  provide  a  realistic 
description  of  the  S/E  interface.  As  expected,  how¬ 
ever,  the  fast  decay  process  of  F,  following  the 
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Fig.  2.  Change  M,  in  surface  electron  density  against  F.ia 
depletion  and  moderate  accumulation  layers.  Straight  line 
(solid)  corresponds  to  y^  >  — 0.06S  V,  dashed  curves  to 
y^  >  —0.065  V  0.025  V.  Corresponding  values  of  the  bar¬ 
rier  height  are  marked  off  at  the  upper  abscitsa. 


Fig.  3.  Eketron  density  vs  barrier  height.  Stars  represent  the 
free  electron  density  JV,  for  a  sample  with  a  relatively  low 
electron  bulk  concentration  n^.  The  classical  curve  (dashed) 
is  calculated  on  the  assumption  of  a  parabolic  conduction 
band  throughout  (eqns  (4)  and  (6)].  The  solid  curve  (“quan¬ 
tized”)  is  based  on  self-consistent  calculations  carried  out 
for  ZnO.  (In  the  depletion  range,  F,  <  0,  A',  is  negative  and 
the  plot  is  that  of  — /V,.)  The  squares  represent  the  leakage 
electron  dei^ty 

activation  of  the  shoiiing  switch  is  not  exponential, 
since  C„  in  the  actual  S/E  interface  is  not  a  passive 
element  but  voltage  dependent.  This  does  not  inter¬ 
fere  at  all  with  the  measurement  of  ^  f ,  Sfn', ,  SN„  and 
Ni,  and  the  results  to  be  presented  below  were 
obtained  precisely  as  describ^  in  connection  with 
Fig.  1(b).  It  should  be  noted  that  the  fut  decay  in  V, 
is  characterized  by  a  decay  time  of  2  to  3  ps,  just  the 
value  to  be  expected  from  the  discharge  of  the  free 
electrons  in  the  accumulation  layer  previously  in¬ 
duced  by  the  applied  pulse  ((r -I- r,)S5C„/5F,,  S 
being  the  area  of  the  CdS  surface  in  contact  with  the 
electrolyte].  Accordingly,  we  identify  the  amplitude  of 
the  fast  component  as  representing  the  induced 
change  6Q^  in  free  space-charge  density.  The  subse¬ 
quent,  much  slower  decay  in  the  F,.  trace  (decay  time 
in  the  millisecond  range)  has  a  low  amplitude,  signi¬ 
fying  the  low  density  of  surface  sutes  (less  than 
10'* cm"’).  This  is  followed  by  an  essentially  non¬ 
decaying  component  representing  the  leakage  charge. 
The  non-decaying  component  cannot  be  seen  on  the 
scale  of  Fig.  1(c). 


RESULTS  AND  DISCUSSION 

The  results  to  be  reported  are  based  on  the  analysis 
of  oscillogram  traces  of  the  type  shown  in  Fig.  1(c) 
above.  Under  equilibrium  conditions,  a  depletion 
layer  usually  exists  at  the  CdS  surface  in  contact  with 
the  electrolyte.  Since  the  measurements  described 
above  yield  directly  only  the  changes  6N, 
in  the  surface  electron  density  as  a  function  of 
dF,-  V,~  F«,  the  first  task  is  to  determine  the 
equilibrium  barrier  height  F^.  Using  eqn  (4)  one  can 
write 

iN.  -  ±(*c£,n,k,  Tlqy^Ffy^  +  Sy.)  -  (10) 
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Fig.  4.  Electron  density  vs  barrier  height.  Stars  represent  the 
Tree  electron  density  N,  for  a  sample  vrith  a  high  electron 
bulk  concentration  n^.  The  classical  curve  (dashed)  is  calcu¬ 
lated  on  the  auumption  of  a  parabolic  conduction  band 
throughout  (eqns  (4)  and  (6)].  The  solid  curve  (“quantised”) 
is  based  on  self-consistent  calculations  carried  out  for  ZnO. 
(In  the  depletion  rang:.  1',  <0,  N,  is  negative  and  the  plot 
is  that  of  -N,.)  The  squares  represent  the  leakage  electron 
density  Nt,. 

For  depletion  and  weak  (non-degenerate)  accumula¬ 
tion  layers,  Fiiy^a  +  SV,)  can  be  evaluated  from  the 
simpler  expression  given  by  eqn  (7).  A  plot  of  the  SN, 
data  against  the  calculated  values  of  F,  should  yield 
a  straight  line  only  if  the  correct  value  of  is  chosen 
in  the  calculation  of  F,.  A  simple  computer  program 
has  been  set  up  to  so  determine  F,,.  The  results  are 
shown  in  Fig.  2,  where  the  choice  of  *=  —0,065  V 
is  seen  to  lead  to  the  required  linear  relation.  The 
accuracy  of  this  method  is  about  ±0.02  V,  the 
limits  for  which  deviations  from  linearity  become 
detectable.  This  is  illustrated  by  the  two  dashed  lines 
in  Fr.  2  which  correspond  to  F,,*— 0.065V 
±  0.025  V,  for  which  marked  deviations  from  the 
linear  relations  are  clearly  apparent.  The  slope  of  the 
straight  line  yields  the  electron  bulk  concentration 
[see  eqn  (10)].  The  value  so  obtained  is  marked  in 
Fig.  2  and  agrees  well  with  Hall  effect  measurements. 

is  derived  from  the  intercept  of  the  straight  line 
with  the  ordinate  (F,  ~  0). 

Once  F,o  and  have  been  determined,  the  com¬ 
plete  N,  vs  F,  curve  (depletion  and  accumulation 
ranges)  can  be  evaluated  from  the  measured  data 
through  the  relations  F,  =  y^  +  5F„  N,^  +  SN,. 

Typical  results  obtained  in  this  manner  are  displayed 
by  the  full  points  (stars)  in  the  semi-log  plots  of  Figs  3 
and  4.  corresponding  to  the  two  extremes  of  the  bulk 
carrier  concentration  in  the  available  resistivity  range 
of  the  samples.  In  the  depletion  range  ( F,  <  0).  N,  is 
negative  but.  because  of  the  logarithmic  scale  used, 
the  plot  is  that  of  its  absolute  magnitude.  The  dashed 
curves  in  both  figures  (referred  to  as  “classical”) 
represent  eqns  (4)  and  (6),  with  the  appropriate 
values  of  and  IF^.  They  assume  parabolicity 
of  the  conduction  band  throughout.  The  solid 
curves  (“quantized")  take  quantization  effects  into 
account.  As  pointed  out  above,  they  are  based  on 
self-consistent  calculations  carried  out  for  ZnO  [6, 7) 


which,  to  a  fair  approximation,  are  adequate  also  for 
CdS.  As  expected.  tHe  classical  and  quantized  curves 
coincide  for  depletion  and  moderate  accumulation 
layers.  In^hese  ranges,  they  are  seen  to  account  well 
for  the  experimental  data.  Such  agreement,  obtained 
with  no  adjustable  parameters,  shows,  among  other 
things,  that  there  is  no  significant  density  of  shallow 
surface  sutes  at  the  CdS  surface.  After  tlw  activation 
of  the  shorting  switch,  such  states  might  release  the 
trapped  elections  induced  by  the  pulse  back  into  the 
conduction  band  very  fast,  in  a  time  comparable  to 
that  of  the  discharge  of  the  free  electrons  in  the 
space-charge  layer.  The  two  processes  would  then  be 
indistinguishable,  with  the  result  that  the  measured 
value  of  SQ,r  in  Fig.  1(c)  would  be  higher  than  the 
change  in  the  free  space-charge  density.  The  good 
agreement  of  the  data  with  the  theoretical  curves, 
which  represent  the  free  charge  only,  rules  out  this 
possibility. 

For  strohg  accumulation  layers,  on  the  other  hand, 
the  experimenul  points  in  Figs  3  and  4  are  seen  to  lie 
well  above  both  of  the  theoretical  curves.  This  can 
only  mean  that  the  density  of  states  in  the  upper 
reaches  of  the  CdS  conduction  band  is  considerably 
larger  than  that  associated  with  a  parabolic  band. 
This  point  will  be  discussed  in  detail  below.  The 
squares  in  Figs  3  and  4  represent  the  electron  leakage 


Fig.  5.  Model  for  the  structure  of  the  CdS  lowest  conduc¬ 
tion  band  (solid  curves),  (a)  Energy  E  vs  wavevector  k.  (b) 
Density  of  states  N{E)  vs  energy  E.  The  dashed  curves 
depict  the  extensions  of  the  parabolic  region  for  £(k)  and 
N{E). 
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Fig.  8.  (a)  Schematic  repieaentation  of  the  expected  te- 
sponse  of  the  voltage  drop*  y,  and  F,  to  a  double  pute.  (b) 
Decay  characteristics  of  the  leakage  charge  for  different 
induced  barrier  heights  F,. 

0.2-10^$.  In  other  words,  the  leakage  takes  place 
practically  instanuneously  following  the  onset  of  the 
pulse  and  is  not  augment^  during  the  pulse.  As  will 
be  shown  below,  it  remains  stored  at  or  near  the 
interface  for  a  considerable  time  (hundreds  of  micro¬ 
seconds).  It  is  gratifying  that  for  both  electrolytes  and 
for  the  different  pulse  durations,  the  results  of  AT,  vs 
y,,  which  take  into  account  the  electron  leakage,  are 
identical.  These  observations  are  compatible  with  the 
equivalent  circuit  in  Fig.  1(a),  and  in  fact  led  to  its 
construction.  The  dashed  curve  represents,  as  previ¬ 
ously,  the  classical  approximation  on  the  assumption 
of  a  parabolic  band  throughout  (constant  effective 
mass).  The  solid  curve  is  based  on  the  non-parabolic 
mode)  (Fig.  S),  and  is  again  seen  to  account  well  for 
the  experimenul  data. 

As  pointed  out  above,  electron  leakage,  the  domi¬ 
nant  loss  mechanism  in  strong  accumulation  layers  of 
CdS,  is  associated  with  the  non-decaying  component 
of  Qm  [Fig.  1(b)].  Hence,  apart  from  the  magnitude 
of  the  electron  leakage,  the  circuit  of  Fig.  1(a)  is 
unable  to  provide  any  information  on  the  decay 
characteristics  of  the  leakage  charge.  To  that  end,  we 
employed  two  pulse  generators  connected  in  parallel 
so  as  to  provide  two  pulses  separated  by  a  variable 
delay  time  /^.  Except  for  the  diode  D  and  the  shotting 
switch  S.  which  have  been  omitted,  the  rest  of  the 
circuit  is  identical  to  that  in  Fig.  1(a).  The  variatioiu 


of  V,  and  V,  with  time  are  Khematically  displayed  in 
Fig.  8(a).  The  first  pulse  charges  up  the  interface 
capacitance  (consisting  of  the  space^charge  capaci¬ 
tance  Ca.  and  the  leakage  capacitance  C^).  At  the 
termination  of  the  pulse  (r «  D,  Kt  doe*  not  yield 
directly  the  value  of  SV,  because  in  the  absence  of  the 
diode  D,  SQ^  now  discharges  through  the  internal 
resisunce  of  the  two  pulse  generators  (r^/2).  Thus,  at 
i^T,V,  assumes  the  value  of  [(r,/2)/(r  +  r,f2)]Sy„ 
from  which  6V,  can  be  evaluated.  The  discharge  time 
r  is  given  approximately  (since  C  >  C^)  by 
(r  r^/2)C^,  which  is  typically  of  the  order  of  several 
microseconds.  As  for  F,,  it  rises  monotonically  dur¬ 
ing  the  first  pulse,  up  to  the  value  of  Q^{T)IC,  Q^{T) 
being  the  total  charge  induced  by  the  first  pulse. 
Thereafter  it  decays  with  the  same  time  constent  of 
T  to  the  level  Qu/C,  Qu  being  the  leakage  charge 
incurred  by  the  first  pulse.  After  a  delay  time  taken 
to  be  larger  than  r  such  that  practically  the  entire 
charge  in  has  decayed  (F,-*0),  the  second  pulse 
is  applied.  Its  amplitude  is  adjusted  so  as  to  rec^rge 
to  the  same  value  of  SF,.  If  is  smaller  than 
the  storage  time  of  the  leakage  charge,  then 
C„  would  still  retain  some  of  its  charge  at  r  k  and 
the  leakage  charge  Cu  measured  after  r^+T-i-r 
would  be  less  than  that  {Qu)  after  T  -t-  r.  By  measur¬ 
ing  Qu  -  Ql\  >  the  charge  retained  at  the  interface  up 
to  the  application  of  the  second  pulse,  as  a  function 
of  the  delay  time  between  the  two  pulses,  one  can 
derive  the  decay  characteristics  of  the  leakage  charge. 
In  Fig.  8(b),  -  Wii  [»  “(Q«  “  Qti )/?)  >*  plotted 

against  on  a  log-log  scale  for  different  values  of  the 
induced  barrier  height  V,.  The  decay  characteristics 
are  most  interesting.  For  a  considerable  length 
of  time  (~500|is  for  F,  — 0.71V,  — 200us  for 
F, «  1.04  V  and  ~70  /is  for  V, »  1.35  V),  the  leakage 
charge  is  almost  fully  retained  at  the  interface 
(A^4,  »  0).  Subsequently,  the  charge  leaves  the  inter¬ 
face  and  its  decay  with  time  is  nearly  linear.  The 
unique  features  of  the  leakage  charge  are:  that  the 
leakage  occurs  in  a  time  shorter  than  0.2 /is  (see 
Fig.  7),  immediately  upon  the  application  of  the 
pulse;  that  no  further  leakage  takes  place  through  the 
pulse  duration  (up  to  at  least  10 /is);  and  that  the 
charge  remains  stored  at  or  near  the  interface  for 
hundreds  of  microseconds.  From  the  experimental 
point  of  view,  these  features  are  most  fortunate  in 
that  the  stored  leakage  charge  can  be  accurately 
determined  and  taken  into  account  in  the  determi¬ 
nation  of  the  N,  vs  V,  relation,  which  is  of  primary 
interest  in  this  paper. 

CONCLUSION 

The  space-charge  characteristics  of  strong  accumu¬ 
lation  layers  of  hexagonal  CdS  show  considerable 
non-parabolicity  of  the  lowest  conduction  band  in 
this  material.  A  simple  model  for  the  structure  of  the 
lowest  conduction  band  (Fig.  S)  accounts  well  for  the 
experimental  data.  It  assumes  a  parabolic  band  up  to 
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an  energy  of  0.12S  eV  above  the  band  edge  and  a 
linear  dependence  of  the  energy  on  the  wavevector  at 
higher  energies.  Such  a  model  is  no  doubt  an  oversim¬ 
plified  represenution  of  the  actual  situation,  but  it 
provides  a  good  idea  as  to  the  structure  of  the  lowest 
band;  and  indeed,  the  density  of  states  function  of  the 
model  agrees  well  with  that  calculated  by  Huang  and 
Ching  [2]  using  a  minimal  basis,  orthogonalized 
LCAO  method. 

It  should  be  noted  that  most  of  the  electrons  in 
strong  accumulation  layen  occupy  states  in  the  non¬ 
parabolic  section  of  the  band.  As  such  they  are 
characterized  by  considerably  higher  effective  masses 
than  at  the  bottom  of  the  band.  Consequently,  quan¬ 
tization  effects  are  expected  to  be  insignificant,  so  that 
our  using  the  classical  calculations  [eqn  (5)]  for 
the  space-charge  characteristics  in  the  comparison 
between  theory  and  experiment  is  justified. 

The  measurements  presented  were  made  possible 
by  the  pulse  techniques  employed  in  the  study  of  the 
S/E  interface.  These  techniques  enable  both  the 
creation  and  study  of  space-charge  layers  at  the  semi¬ 
conductor  surface,  ranging  from  large  depletion  to 
huge  accumulation  conditions.  What  is  more  impor¬ 
tant,  they  permit  a  straightforward  separation  of  the 
different  components  of  the  induced  charge  at  the  S/E 
interface,  so  as  to  obtain  the  required  variation  of  the 
free  electron  density  N,  with  the  barrier  height  V,.  At 
the  same  time,  the  measurements  yield  the  electron 
leakage  across  the  interface,  which  becomes  promi¬ 
nent  in  strong  accumulation  layers.  While  a  determi¬ 
nation  of  such  leakage  is  essential  for  deriving  the 
experimental  N,  vs  V,  curve,  the  leakage  process  is  of 
considerable  interest  in  itself  The  observation  that 
the  leakage  increases  nearly  exponentially  with  bar¬ 
rier  height  and  hence  with  interface  field  suggests  that 
electrons  tunnel  from  the  accumulation  layer  into 
some  ionic  species  (such  as  H*)  in  the  electrolyte, 
possibly  at  sharp  points  and  other  imperfections  at 
the  interface.  The  leakage  lakes  place  practically 
instantaneously  (in  less  than  0.2  ps)  following  the 
onset  of  the  pulse  and  no  further  leakage  occurs 


during  tlw  pulse  (up  to  at  least  10  ps  duration — see 
Fig.  7).  This  suggesU  that  the  density  of  the  ionic 
species  is  relatively  low:  the  tunneling  electrons 
appear  to  saturate  completely  all  the  active  ionic 
species  within  a  tunneling  distance  corresponding  to 
the  induced  interface  field.  Another  obs^ation  is 
that  the  leakage  remains  stored  in  the  electrolyte,  for 
a  considerably  long  time  (see  Fig.  8).  The  storage  time 
probably  represenu  the  time  it  takes  for  diffusion  to 
replenish  the  exhausted  ionic  species  near  the  inter¬ 
face.  The  storage  time  of  several  hundreds  of  micro¬ 
seconds  is  not  unreasonable  for  such  a  diffusion 
process. 
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ABSTRACT 

The  nature  of  the  percolation  process  in  granular  metals  is  examined  for  the 
first  time  by  a  computer  simulation  of  a  system  of  metallic  grains  embedded  in  an 
insulating  matrix.  Assumiim  that  the  intergrain  conduction  is  due  to  quantum 
mechanic^  tunneling  it  is  found  that  a  p«rcoiation-like  critical  behavior  of  the 
conductivity  is  obtained,  but  that  a  percolation  utuvmal  behavior  will  be  found  only 
in  a  very  special  case.  In  contrast,  the  behavior  of  the  electrical  noise  does  not 
deviate  substantiallv  from  the  universal  one.  Comparison  of  these  results  with 
recent  ejqjerimental  observations  suggests  that  in  the  metallic  range,  both  transport 
properties  are  controlled  tty  the  continuous  metallic  network  rather  than  by 
mtergrain  tunnelii^  We  propose  that  the  metallic  network  resembles  the  previously 
studied  system  of  Averted  random  voids". 


INTRODUCTION 

About  twenty  years  ago  Abeles  and  coworkers  [1]  initiated  a  comprehensive 
study  of  the  electronic  properties  of  granular  metals,  by  studying  the  me^  content 
and  tempnature  dependence  of  their  electrical  conductivity,  o(x,T).  As  a  result  of 
these  studies  they  have  attributed  the  electrical  conductance  in  the  materials  of  low 
metal  content  to  quantum  mechanical  tunneling  under  the  constraint  of  metal  grain 
charging  [2],  and  uie  conductance  of  the  materials  of  high  metal  content,  to  simple 
metallic  transport  in  the  network  formed  by  the  coalescence  of  the  metallic  grains 
[3].  Correspondingly,  in  the  latter  case  a  percolation -like  U]  behavior  is  to  be 
expMted  [5]  where  me  percolation  threshold,  x^  ,  is  associateo  with  the  onset  of  a 
continuous  metallic  network.  While  the  more  specific  nature  of  the  conduction 
mechanism  for  x<Xg  was  never  settled  [6],  the  conduction  in  the  x>Xc  range 
seemed  until  recently  to  be  well  accounted  for  by  the  percolation  process  mentioned 
above.  The  most  convincing  piece  of  evidence  to  support  the  percolation 
conduction  was  the  power— law  de^ndence: 

u  «  (X-Xg)*  (1) 

found  {3]  in  the  W— AI2O3  system,  with  an  exponent  t=1.9,  which  is  very  close  to 
the  universal  value  of  tg  »  2.0  predicted  for  three  dimensional  resistor  networks  by 
percolation  theoiy  [44|].  That  a  percolation  geometty  is  formed  Ity  the  coalescence 
of  metallic  particles  u  supported  by  two  observations.  The  first  is  that  in  a 
composite  made  of  macroscopic  silver-coated  glass  spheres  which  are  embedded  in 
Teflon  (U.  where  physical  contact  rather  than  tunneling  between  the  "metallic" 
spheres  takes  place)  the  same,  universal,  behavior  was  found  [7].  The  second 
observation  is  that  a  transition  firom  thermally  activated  behavior  to  a  metallic- 
temperature  dependence  of  the  conductivity  takes  place  [1]  at  metal  contents  dose 
to  Xq  ,  udien  x^  is  determined  by  fitting  the  e]q>erimental  data  to  Eq.  (1).  Even  the 
problem  of  "unexpected"  high  x^  values  obtained  eroerimentally  [3j,  wmch  seem^ 
to  cast  some  question  on  the  percolation  description,  was  remo^  recently  [8]. 
This  removal  is  based  on  the  argument  that  such  metallic  networks  are  better 
described  by  a  structure  which  approximates  a  randomly  dosed-padced  tystem  of 
conducting  spheres  than  by  a  three  dimensional  (Scher  and  Zallen-Iike)  system  of  a 
mixture  of  ccmducting  and  insulating  spheres,  or  an  effective-medium -like  [3,9] 
mixture  of  the  two  phases. 

Recently  there  was,  however,  an  attempt  by  Mantese  et  aL  [9,10]  to  describe 
the  metallic  conductivity  in  P{-Al203,  for  the  range  x>x^  ,  tty  an  effective  medium 
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theory  which,  a  priori,  should  yield  for 
the  exponent  in  Eq.  (1).  It  was  also  argued 
[9]  on  the  basis  of  that  theory  and  the 
ea^rimentalty  found  [10]  x  depepdence  of 
the  resistivity  nonnalixed  noise  [11],  in 
the  above  material,  that  even  tot  t>>i^  , 

Sb  is  dominated  by  intergrain-tuimeimg 
conduction  (rather  than  by  metallic 
transport).  In  view  of  the  contradiction 
between  the  effective  medium  and  the 
percolation  approaches,  the  first  question 
that  arises  u  whether  the  eroe^ental  g  io* 
power-law  behavior  reported  by  Abeles  et  aL 
[3]  is  specific  to  gi^ular  W-AI2O3  or  o 
whether  it  is  a  general  phenomenon  in 
.  manuiar  metals.  Fiamining  the  literature  we 
did  not  find  any  other  presentation  of 
conductivity  data  of  manular  metals  that 
indicate  a  power— law  oehaviour  in  general, 
and  a  percolation  universal  behavior,  of  t^tg 
in  paracular.  Following  then  the  possibiiity 
that  such  a  percolation-like  behavior  was  a 
unique  observation  we  have  carried  out 
experimental  measurements  cm  another 
granular  tystem,  Le.  Co-Si02.  The  results 
of  the  correnKmding  measurements  are 
shown  in  Fig.  1  and  are  ccmipared  to  the 
results  [3]  cm  W-AI2P3.  Our  results  clearly 
confirm  the  universal  oehavior  for  x>x^  , 
and  thus  the  suggestiems  of  Abeles  et  ui.  |3] 
of  a  percolatum-like  rather  than  an 
effective-medium-like  behavior.  More 
details  of  the  Co-SiO?  granular  tystem  and 
the  temperature  dependence  of  its 
conductivity  were  presented  previously  [12]. 
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conductivity  of  granular  Co-Si02 
(as  obtains  both  at  78  K  and  300 
K)  as  a  function  of  (x-Xg)/Xg  , 
where  x  is  the  volume  fraction  of 
metal  and  x.  was  obtained  from  a 
best  fit  of  the  data  to  Eq.  (1).  The 
dashed  line  represents  the  data  of 
Abeles  etaL,  ret.  3. 

Since' the  ab^  eaqierimentai  e^dence  provides  quite  a  ccmvincing  proof  for  a 
percolaticm-like  conduction,  the  question  arises  whether  the  effective  medium 
picture  with  the  intergrain  tunneling  transport  is  indeed  necessary  to  account  for  the 
resistivity  normalized  noise,  as  suggests  in  Ref.  9.  The  present  computer 
simulation  is  trying  to  answer  this  question  by  stuclying  the  transport  properties  of  a 
tystem  of  metallic  spheres  between  wfrich  a  simple  quantum  tunneling  conduction 
takes  place.  The  corresponding  model,  which  is  interesting  in  its  own  right,  is 
studiea  here  for  the  first  time.  The  rraults  appear  to  show  that  for  x>X(.  ,  both  , 
the  conductance  and  the  resistance  noise  are  dominated  by  the  metallic  network 
transport  Moreover,  we  propose  that,  at  metallic  contents  just  above  x>  ,  the 
coalescence  of  the  grains  forms  an  "inverted  random  void"  like  network  [1344],  i.e. 
a  tystem  of  intersecting  permeable  q>heres. 


THE  MODEL 

Our  mcxiel  consisti  of  hard  spheres  of  diameter  b  which  are  randomly 
dispersed  in  a  cube^  a  unit  volume.  Hence,  for  N  spheres  the  fractional  occupied 
volume  is  x«(ii/$)b^.  In  the  sunulations  the  N  tyiheres  are  placed  randomty  one 
after  the  other.  A  new  »here  is  accepted  to  the  tystem  only  u  it  does  not  overly 
a  previous  sphere.  For  the  precision  used  in  the  simulations  for  the  coordinates  of 
the  spheres*  centers  there  is  no  geometrical  contact  between  any  pair  of  spheres. 
Thus,  up  to  a  concentration  which  corresponds  to  the  randomty  close  packed 
structure  [8]  (x«0.64)  the  tystem  has  no  coimectivity.  If  we  assume  that  the  above 
tystem  simulates  metallic  partides  embedded  in  an  insulating  matrix,  all  the 
4>herical  partides  can  be  considered  electrically  connected  due  to  quantum 
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mechanical  tunneling.  For  simplicity  we  consider  simple  tunneling  or  a  constant 
tunneling  activation  energy  [1,2]  (Le.  we  ne^ect  the  interaction  energy  between  the 
spheres  or  we  assume  a  nigh  enoueh  temperature  [9])  which  is  a  reasonable 
approximation  for  a  system  >^ere  all  the  spheres  are  of  the  same  diameter  [1].  We 
assume  then  that  the  conductance  between  two  metallic  grains  can  be  written  as: 

g  •  goe*P(-*^o)  »  (2) 

where  go  is  a  constant  (taken  as  unity],  s  is  the  separation  (at  the  nearest  gap) 
between  two  spherical  particles  and  Tq  is  the  characteristic  tunneling  distance  (the 
decay  distance  of  the  electrons*  wave-function  (Q).  By  assigning  a  conductance  g 
to  each  pair  of  metal  spheres  we  have  oonstruetea  a  network  of>N(N-l)/2  resistors. 
This  model  Is  expeeteo  to  be  a  simple  simulation  of  a  granular  metal  when  x<x^ 
and,  for  the  ensemble  of  particles  which  do  not  belong  to  the  percolation  cluster 
when  x>Xo  .  Computationally,  however,  if  reasonable  statistics  are  desired,  the 
computer  memory  r^uirements  become  difficult  and  expensive.  The  difficulty  is 
reduced  [151  if  one  assumes  a  cutoff  distance  s^  such  that  for  5>S(.  one  takes  g»0. 
The  physical  justification  for  the  cutoff  is  the  fact  that  the  average  dist^ce  be^een 
the  nearest-neighbor  spheres’  centers  in  the  system  is  [16]  d=(&iTN)^=bAt^  and 
one  expects  then  that  resistors  for  which  s>Sp>>a  provide  a  negligible 
contribution  to  the  conductance  of  the  system.  We  set  s^  equal  to  6b.  For  typical 
granular  metals  [1],  r^  ^  b/S  and  d  b,  and  thus  settmg  S(.==6b  amounts  to 
neglecting  conductances  which  are  smaller  than  the  typical  largest  conductance  in- 
the  system  ^  a  factor  of  exp[-5(7-dA})]. 

Following  the  construction  of  the  system  we  have  computed  the  currents 
through  all  the  different  conductors,  using  our  conjugate-gradient  algorithm  [14,15], 
and  computed  the  resultant  conductance,  G,  and  the  resistivity  normalized  noise  So. 
For  simplicity  we  considered  the  case  where  the  resistance  fluctuations  are 
proportional  to  the  value  of  the  resistors.  Correspondingly  the  resistance  noise  is 
taken  here  as  [11,13]  ^  o 

Sr  =  A{Si„Vg„2y(Si„Vg„)2  ,  (3) 

where  A  is  a  constant  and  the  summation  is 
over  all  the  resistors  in  the  system.  This  ^  ® 
form  has  the  essentials  of  the  problem,  since 
it  contains  the  dtverginE  nature  of  die 
resistor  distribution.  Li  view  of  the 
exponential  dependence  of  g  [Eq.  (2)]  it  is  e 
not  eiqiecmd  that  intersrain  volume  effects  1 
[13,14]  will  alter  significanUy  any  of  the  ^ 
conclusions  of  the  present  work.  f 

O  0.4 

RESULTS 

The  typical  behavior  of  the  computed 
sam|>le  conductance  G,  as  a  function  of  the  o.2 
metallic  fractional  volume  x,  is  shown  in  Fig. 

2  for  two  values  of  rjb.  It  is  seen  that  the 
conductance  increases  monotonically  with  x, 
and  that  this  increase  is  faster  the  smaller  o.o 
the  value  of  r~  As  is  to  be  expected  from 
the  model,  tnere  is  no  cridou  x^  below  x 

sdiich  G«0  (as  in  well  defined  percolation  FIGURE  2.  The  computed 
qrstems).  On  the  other  hand,  the  qualitative  dependence  of  the  conductance  G 
features  of  the  data  in  He.  2  are  similar  to  of  a  system  of  metallic  spheres  in  a 
those  found  experimentally  in  the  various  unit  volume,  on  the  volume  fraction 
granular  metals  [3,9,12].  Considering  the  of  metal,  for  two  extreme  values  of 
results  of  Hg.  1  it  is  apparent  that  a  log-log  r^,  the  characteristic  tunneling 
presentation  of  the  data  given  in  Fig.  2,  for  x  c&tance  . 
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larger  than  some  jl.  ,  may  indicate  whether 
our  simple  tunneung  model  resembles  a 
percolation-like  behavior  [Eq.  (1)].  Indeed 
the  best  fit  of  these  data,  for  the  higher  x 
^ues,  to  a  power-law  behavior  yields  a 
"threshold”  vuue  of  x^^O.!  and  a 
corresponding  exponent  For  the  illustration 
of  the  nature  of  this  fit  we  show  in  3  a 
log-log  presentation  of  the  data  of  Fig.  2. 
The  rc^ts  show  that  for  the  higher  x  vuues 
the  data  can  be  well  fitted  by  a  power-law 
and  cannot  be  distinguished  from  typical 
percolation  simulation  data.  ]Bven  the 
deviation  from  a  power-law  behavior  at  low 
x-Xg  resembles  the  typical  behavior  due  to 
finite  size  scaling  [4,14]  and  thus,  for  the 
finite  samples  used  here,  this  deviation  does 
not  reveal  a  non-percolative  conduction. 
The  t  values  found  are,  however,  strongly  r^ 
dependent  Qualitatively  this  is  to  be 
expected  for  the  small  values  of  rjb  since 
for  these  values  onty  nearest  neighbors 
contribute  to  the  electrical  conduction  of  the 
system  and  thus  it  resembles  conduction  in  a 
percolating  system.  The  increase  of  t  with 
decreasing  r^  indicates  then  a  transition  from 
a  conductance  due  to  the  many  parallel 
diains  of  a  few  resistors  of  high  values,  to  a 
conductance  due  to  very  few  (percolation- 
like)  long  chains  of  many  resistors  of  smaller 
values.  Indeed,  for  an  experimental  system 
where  nearest  neighbor  (percolation -like) 
tuimeling  transport  tekes  place  between 
"metallic^  spheres  which  do  not  coalesce 
(Mogul-L  m  PVQ,  we  have  found  [16] 
previously  that  a  power  law  behavior  with  a 
nonuniversal  exponent  (t>tg)  is  obtained. 
We  may  conclude  then,  that  in  systems 
where  the  tunneling  is  of  short  range,  the 
conductance  behaves  according  to  the 
continuum  theory  of  percolation  [5,13,16].  In 
fact,  it  appMrs  that  we  may  define  a 
phenomenological  critical  value  for  r^,  which 
marks  the  transition  between  a  non- 
percolative  (t<tQ)  behavior  to  a 
monuniversal)  percolative  behavior  (t>tQ). 
This  value,  in  the  present  simulations  is 
found  to  be  rQQ«0.1b.  To  further  stress  the 
relationship  between  the  present  problem  and 
continuum  percolation  [5,13]  we  show  in  Fig. 
4  the  conductance  distributions,  P(g),  vdiia 
were  found  here  for  various  values  of  r^.  It 
is  dearly  seen  that  these  distributions  exhibit 
a  diverging  behavior  reminiscent  of  the 
conductance  distribution  fimetions  [5] 
P(g)*(i-o)g  ®  which  yields  the 
nonuniversal  critical  behavior  for  other 
continuum  systems  [13,14].  This  behavior  is 
known  [5]  to  be  manifested  by  the  relation 
t*tQ+a^l-a).  Hence,  when  the  system 
approadies  a  percolation  -  like  network 
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FIGURE  3.  The  data  of  Fig.  2 
presented  on  a  log(l/G)  vs.  logux- 
x^/xJ  scale.  The  apparent  cntical 
metallic  volume  fraction  is  X{.»0.1. 


10*  r 


10*  b 


10*  L- 


10* 


•I  '  ■  "  I - ' — ^ — '  I  '  ' 

Tunneling  Model 

s,»3b  N-7500 

X-0.25  b-0.04 

o 


•  U-0.05b 

•  r,-0.1b 
D  r,-0.15b 
o  r,«0.2b 

•  r»-b 


o 

o 


10-* 


■  ■  I  ■  ■  ■  ■  I 


0.1 


j__j — i_i 


FIGURE  4.  The  conductance 
distributions  P(gl  of  the  present 
simulation  for  a  tew  tunneli^  cutoff 
distances,  as  marked  in  the  f^ure. 
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(r^-^)  the  nonimiversality  is  due  to  the  singly-connected  (nearest-neighbor)  largest 
re^ors,  as  in  "conventional'’  continuum  percolation  [543>16].  In  summary,  the  data 

E resented  in  Figs.  2-4  show  that,  in  principle,  data  derived  from  a  ^em  of  grains 
etween  which  simple  tunneling  takes  place,  can  always  be  presented  in  a  form 
which  suggests  a  continuum  percolation- like  behavior.  This  presentation  is  both 
conceptu^  [5]  (see  above)  and  experimentally  {1^  justified  m  the  r^/b-^  limit. 
We  also  earned  out  computations  for  various  values  of  the  cut-ofi  distance  s^  in 
the  range  b  £  ^  £  6b  and  no  variations  of  significance  in  the  present  context  were 
ol^rved  [15].  Hence  we  do  not  discuss  these  variations  here. 

Turning  to  the  resistivity  normalized  noise  data,  we  found  again  (for  x>Xj^,  a 
power-law  oehavlor  of  the  form: 


Sr  «  (x  -  Xc) 


— K 


(4) 


This  is  shown  in  Fig.  S.  As  expected  from 
the  conductance  data,  the  vmues  of  k 
increase  with  decreasing  r^  but  they  are 
found  to  be  confined  to  a  relatively  narrow 
range  around  1.6  (k  ==1.6±0.5).  Noting  that 
the  universal  k  value  is  k^— 1.56,  the 
surprising  fact  is  that  the  deviation  k-Kq  (for 
K>Kg)  is  much  smaller  than  t-t^  for  the 
same  continuum  system.  This  is  just  the 
opposite  to  the  situation  in  ^ical  continuum 
s;^ems  [5,13].  Following  the  interpretation 
of  the  conductivity  data,  this  implies  that 
since  the  currents  through  the  many 
"parallel"  high  resistance  channels  are  small, 
tne  noise  power  is  dominated  by  the  singly 
connected  bonds  through  which  the  entire 
current  is  passing  as  in  a  genuine  percolation 
backbone.  Hence,  while  the  value  r^^. 
describes  the  transition  to  the  percolation - 
like  behavior  of  the  electrical  resistance,  the 
corresponding  value  for  the  resistance  noise 
is  prooabiy  smaller.  The  reason  for  this  may 
be  that  as  Tq  decreases  below  the  second 
moment  of  the  current  (i.e.  tne  resistance 
[5])  through  a  channel  made  of  nearest 
neighbors  dominates  the  multiplicity  of  the 
non -nearest  neighbors  current  paths  in  the 
system,  but  for  the  fourth  moment  (i.e.  the 
electrical  noise  [5,13])  this  may  not  yet  be 
the  case. 
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normalized  noise  So  as  a  function 
of  (x-Xp)/Xj.  for  X  >  0.1.  The 
values  of  k  (and  Xp)  are  marked  in 
the  figure. 


DISCUSSION 

The  first  conclusion  from  our  simulation  results  is  that  the  critical  behavior  of 
a  tystem  of  grains,  between  which  tunneling  transport  takes  place,  i.e.  the  value  of 
t,  depends  on  tg/b.  Relating  our  results  to  granular  metals  we  note  that  the  values 
of  b  vary  significantly  from  one  granular  material  to  another  [1].  We  should  expect 
then  that  smee  the  universal -like  behavior  (t’^O  will  appear  only  for  the  specific 
value  of  r^^,  different  materials  should  exhibit  different  t  values.  The  fact  tnat  in 
granular  metals  one  seems  to  get  a  universal  behavior  (Fig.  1  and  Ref.  3)  indicates 
that  a  different  mechanism  is  responsible  for  this  behavior,  namely  percolation.  We 
thus  reconfirm  the  original  suggestion  of  Ref.  3  that  in  the  x>Xp  range,  metallic 
network  percolation  appears  to  m  a  better  description  of  the  system  than  intergrain 
tunneling,  which  can  not  account  for  the  observed  universal  behavior. 

Let  us  now  compare  our  results  with  those  of  Mantese  et  al.  [9,10].  They 
concluded  that  even  for  x>>Xp  the  resistance -noise  is  associated  vrith  the 


intennw  tunneling  rather  than  with  the  continuous  metallic  network.  This 
concwsion  is  based  on  the  asumptions  that  the  relative  resistance  noise  associated 
with  intergrain  tunneling  is  independent  of  x  and  is  much  larger  than  the 
correspondmg  noise  assoaated  with  the  metallic  network.  In  contrast,  our  results 
show  that  the  intergrain  relative  noise  is  strongly  x  dependent  In  addition,  while 
the  metallic  grains  contribute  maybe  onW  a  low  relative  noise,  the  metallic  network 
may  yield  a  high  relative  noise,  due  to  the  metallic  necks  between  the  grains.  This 
noise,  in  principle,  may  be  even  larger  than  that  associated  with  intergrain 
tunneling.  Another  discrepancy  is  the  ^ue  of  the  resistivity -normalized -noise 
exponent  k.  The  results  of  our  model  show  (see  Fig.  5)  that  intergrain  tunneling  is 
reasonably  well  approximated  by  the  known  [13,1^  universal  percolation  behavior 
K/t:^l.  In  contrast,  the  experimental  findings  of  Mwtese  et  al.  [9,10]  show  a  veiy 
strong  nonuniversal  behavior  manifested  by  for  Pt-Al203  and  K/t=4.22  for 

M0-AI2O3.  Thus  we  must  conclude  that  these  experimental  results  may  be  better 
explained  ^  conduction  in  a  metallic  network  than  by  intergrain  conduction. 
Further  support  for  this  argument  can  be  found  in  their  experimentxd  observation  of 
the  change  in  the  temperature  dependence  of  the  relative  resistance  noise  at  the 
same  x.  as  the  transition  in  the  temperature  dependence  of  the  electrical  resistivity 
to  metmllc— like. 

The  question  now  arises  what  kind  of  a  metallic  network  can  cause  the 
observed  universal  value  [3]  for  t  and  nonuniversal  value  [10]  for  k.  So  far  the  only 
system  known  to  obey  these  two  conditions  simultaneously  is  the  inverted  random 
void  model  [13,14]  '(where  overlap  of  permeable  conducting  spheres  takes  place)  for 
which  in  three  dimensions  [14]  t=tQ  and  k=Kq+3  i.e.  K/t=2.36.  Noting  in  particular 
the  agreement  with  one  ot  the  experimental  values  [10]  it  appears  that  upon  the 
coalescence  of  the  metallic  particles,  for  x  just  above  ,  tne  "necks”  which  are 
formed  resemble  overlaps  ot  two  permeable  spheres.  Hence  we  believe  that  for 
x>Xg  both  the  conductance  and  tne  corresponding  noise  are  associated  with  ^e 
inverted-void-like  metallic  network.  It  is  only  for  xC^  that  the  electrical 
behavior  of  granular  metals  is  determined  by  intergrain  tunneling. 
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The  photocarrier  grating  (PCXl)  technique  is  used  for  the  first  time  in  the  high  electric  field 

regime.  The  results  are  shown  to  confirm  previously  unproven  theoretical  predictions 

for  this  regime.  It  is  demonstrated  that,  by  application  of  the  PCO  technique  in  the  high-field 

#  regime,  accurate  values  for  the  ratio  of  the  two  carriers  mobility-lifetime  products  can  be 
deduced.  This  is  in  contrast  with  the  fact  that  their  sum  cannot  be  derived  accurately  from  the 
measurement  of  photoconductivity  because  one  cannot  determine  accurately  the  carrier 
generation  rate.  Combining  the  ratio,  determined  by  the  high  field  PCG,  with  the  low-field 
yields  accurate  values  for  the  mobility-lifetime  products  of  both  carriers. 


I.  INTRODUCTION 

The  two  most  informative  parameters  which  charac¬ 
terize  the  photoelectronic  properties  of  a  semiconductor 
are  the  mobillty*iifetims,  fir,  product  of  the  two  charge 
carriers.'**  These  parameters  contain  the  information  re¬ 
garding  the  transport  mechanism  (/r)  and  the  recombina¬ 
tion  kinetic  mechanism  (r).  Traditionally  /ir  products 
have  been  determined  under  transient  conditions^  or  under 
steady-state  conditions.^  It  is  well  known  that  the  two 
types  of  measurements  yield  different  results, but  in  de¬ 
vices  it  is  usually  the  steady-state  ht  value  which  deter¬ 
mines  their  performance.'*^  For  both  basic  physics’  and 
device  physics*  it  is  important  in  many  cases  to  determine 
these  values  accurately.  This  paper  proposes  and  demon¬ 
strates  a  new  technique  for  accurate  determination  of  the 
/Ur’s  of  both  carriers. 

While  by  now  it  is  quite  established  that  in  low  mobil¬ 
ity  materials  the  photocairier  grating  technique  (PCG) 
yields  accurate  values  for  the  minority  carrier  ^he  carrier 
with  the  lower  drift  mobility"*)  steady-state  fir  product,’ 
there  is  no  accurate  technique  to  derive  the  same  product 
for  the  majority  carriers.  Thus,  the  latter  quantity  has  been 
derived  thus  far  from  the  measurement  of  the  photocon¬ 
ductivity  or.  Assuming  that  the  electron  is  the  majority 
carrier,  the  relation  between  the  above  quantities  is  ex¬ 
pressed  by  the  equation"* 


c^gilier!!  +  likTp)G 


(1) 


and  the  inequality  >  Ph'^'p-  Here  q  is  the  (positive) 
electronic  charge,  G  is  the  carriers  generation  rate, 
Heifip)  is  the  band  mobility  of  the  electron  (hole),  and 
T  n  (r  ^)  is  the  recombination  time  for  the  electron  (hole). 
The  experimental  difficulty  with  the  application  of  Eq.  ( 1 ) 
is  that  a  and  in  particular  G  cannot  be  determined  accu¬ 
rately.  The  difficulties  are  associated  with  inaccurate  deter¬ 
mination  of  the  sample  geometry,  the  quantum  efficiency  of 
the  carrier  generation,  the  absorbed  light  fraction  (even 
after  correcting  for  reflectivity  and  scattering),  and  the 
inhomogeneity  of  the  impinging  light  beam.  Usually,  when 


one  derives' a  /rr  value  from  photoconductivity  [by  using 
Eq.  (1)  or  similar  expressions'*^],  one  can  only  say  that 
this  value  is  within  a  factor  of  two  of  the  “correct”  value. 
This  problem  may  become  severe  if  one  wants  to  determine 
when  a  I"  this  ca.se,  one  determines 

from  another  measurement  (e.g.,  the  low  field 
PCG")  and  with  the  inaccurate  sum  [Eq.  (1)]  one  may 
obtain  a  “negative”  value  for  J.  Another  example  for 
the  need  to  overcome  the  inaccuracies  is  in  the  problem  of 
the  ratio  between  the  steady  state  /itr  value  and  the  tran¬ 
sient  fiT  value."*  It  may  be  that  this  ratio  is  twice  as  large 
or  twice  as  small  as  that  reported  in  the  literature.  **‘ 

Following  these  considerations  it  appears  that  there  is 
a  need  for  a  technique  in  which  the  uncertainties  in  the 
measurement  will  be  removed.  One  of  the  major  advan¬ 
tages  of  the  photocarrier  grating  (PCG)  technique""'*  is 
that,  one  measures  photoconductivity  ratios,  rather  than 
absolute  values  of  the  photoconductivity,  under  the  same 
G.  This  technique  has  been  used,  however,  thus  far  in  the 
zero-field  limit,’*"  where  it  yields  essentially  the  minority 
carriers  fir,  and  in  the  low-field  (or  "relaxation”)  regime, 
where  it  yields  an  ill-defined  fir  product.'’*'’  It  can  be 
shown'’  that  the  latter  is  a  result  of  the  fact  that  in  the 
low-field  regime  the  measured  quantity  (P,  see  below)  is  a 
complicated  combination  of  the  two  carriers  fir  products 
and  other  parameters  and  thus  it  cannot  provide  unique, 
single  carrier,  fir  values.  In  this  paper  we  examine  the 
previously  unstudied  high-field  regime  and  show  that  in 
this  regime  one  can  derive  a  unique  value  for 
;i*T*/(|i,T  J).  A  combination  of  this  result  with  the 
fS/,T^  result  derived  from  the  zero  field  PCX3  measurement 
then  yields  an  accurate  value  for 

II.  THEORETICAL  BACKGROUND 

Turning  to  the  utilization  of  the  PCG  method  for  the 
present  problem,  we  follow  the  concepts  of  the  zero-field 
PCG  theory'*  while  introducing  the  specific  variations  and 
definitions  associated  with  the  high-field  regime.  A  more 
detailed  account  of  the  analysis  and  the  physical  picture  of 
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the  photocarrier  grating  under  these  conditions  will  be 
given  elsewhere.”  In  the  PCX3  configuration,  two  light 
beams  which  can  interfere  generate  carriers  in  a  photocon¬ 
ductor."  When  they  are  made  to  interfere,  a  sinusoidal 
carrier  generation  along  the  direction  between  two  copla- 
nar  electrodes  (hereafter  the  x  direction)  takes  place. 
When  they  are  made  not  to  interfere,  a  uniform  illumina¬ 
tion  results.  The  ratio  between  the  additional  photocon¬ 
ductivities  which  result  due  to  the  application  of  one  of  the 
beams  (the  low-intensity  beam),  in  the  above  two  situa¬ 
tions  is  the  quantity  which  is  being  determined"  experi¬ 
mentally,  /9.  The  photocarrier  generation  function  is  given 
then  by  C  =  Go  -f  AG(x),  where:  AG(x)  =  AgCO&(kx). 
The  generated  electron  concentration  distribution  which 
results  is  n{x)  =  Hq  +  An(x),  where  Ar(x) 
=  Af  cos(Ajc  -I-  v),  and  the  corresponding  hole  concentra¬ 
tion  distribution  is  p(x)  =  po  +  Ap(x),  where  Ap(x) 
=  Ah  cos(kx  -f-  0).  Here  Go(Ro?Po)  >s  the  carrier  genera¬ 
tion  (electron  concentration,  hole  concentration)  when  no 
grating  is  formed,  k  =  2ir/A  where  A  is  the  grating  period, 
and  V  (d)  is  the  phase  shift  of  the  electron  (hole)  concen¬ 
tration  grating"  with  respect  to  the  carrier  generation  grat¬ 
ing  (^  =  V  =  0  with  no  field).  Ag  is  the  amplitude  of  the 
generation  grating  and  A,(Ah)  is  the  excess  electron  (hole) 
concentration  amplitude  in  the  presence  of  the  grating.  The 
theory  developed  thus  far"’'*'"  is  applicable  only  if  the 
linearization  of  the  continuity  equation  applies,  i.e.,  if 
Ag<Go,  and  if  the  material  is  assumed  “intrinsic**  in  the 
sense  that  rq  and  po  are  much  larger  than  the  carrier  con¬ 
centrations  in  the  dark.  The  task  of  the  theory  is  then  to 
present  A,,  Ah,  4>,  and  v  in  terms  of  the  pr  products  while 
the  task  of  the  experimental  procedure  is  to  correlate  the 
former  parameters  with  the  quantity  measured  by  the  PCG 
technique,  p.  For  the  presentation  of  the  experimental  re¬ 
sults  one  uses  the  normalized  grating  amplitude Yeir 
which  can  be  expressed  as'* 

lit[=  I  (M,  cos  V  +  Ph^h  cos  d)^ 

■f  (M^,sin  v  +  ^V<*sin  d)^l/[^*o/Gol  (2) 

and  is  related  to  ^  by"  '*’" 

yfff=(l -/3)/(2n^),  (3) 

where  Yo 's  an  experimental  quality  factor  (of  the  order  of 
unity")  and  y  is  the  light  intensity  exponent  of  the  pho- 
toco.nductivity.'* 

Under  the  small  signal  grating  conditions  the  linear¬ 
ization  of  the  transport-continuity  equations  is  possible  and 
for  the  electrons  the  corresponding  equation  takes  the 
form'^'^ 

/5^Ar\ 

-  (Pt9no/e){Ap/eg  -  Ar/0,)  -  D,  r^j 

fdAn\ 

\  +  (.An/t„  +  Ap/ti,)=AgCOs(kx).  (4) 

The  equation  for  the  holes  takes  a  similar  form.  The 
first  term  in  Eq.  (4)  is  the  space  charge  term,  e  is  the 
dielectric  constant,  and  the  d*s  are  the  free  to  total  electron 
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(or  hole)  ratios.'*  The  second  term  is  the  diffusion  term 
and.  we  assume  throughout  that  Z),  =  {kgT/g)p^  where 
kgT  is  the  thermal  energy.  The  third  term  is  the  drift  term 
where  £o  is  f^e  applied  dc  field,  and  the  last  term  on  the 
left-hand  side  of  (4)  is  the  small  signal  recombination 
term'^  (characterized  by  the  small  signal  recombination 
times  r,  and  tp).  For  aU  practical  purposes'*  t„  =  r  if  and 
=  rf.  Equation  (4)  and  its  coupled  hole  equation  can  be 
solvjd  analytically  by  substituting  the  above  suggested  so¬ 
lutions  in  these  equations,'*  but  the  expressions  for  A^ 
A„  and  v  are  so  cumbersome  that  they  are  of  no  prac¬ 
tical  use  for  deriving  unique  pr’s  from  the  experimental  p. 
The  alternative  approach'*  is  to  perform  the  experiments 
under  conditions  which  correspond  to  a  simplified  mathe¬ 
matical  version  of  Eq.  (4)  and  its  hole  counterpart,  such 
that  a  desirable  unique  and  simple  relation  between  the 
above  mentioned  quantities  will  result.  Such  a  simplifica¬ 
tion  occurs  if  the  space-charge  term  is  made  negligible 
compared  with  the  diffusion  term  and  the  diffusion  term  is 
made  small  but  non-negligible  compared  to  fli  •  dritt 
term.'*  If  we  assume  that  Eq.  (4)  has  the  sinusoidal  solu¬ 
tions  suggested  above,  the  condition  mentioned  can  be 
written  [see  Eq.  (4)]  as 

(qno/€e„)  [  {pfi^ph^p)  -  1  ]  <t^kgT/q  <  kE^.  ( 5 ) 

Under  condition  (S),  the  algebraic  equation  obtained 
by  using  the  suggested  solutions  of  Eq.  (4)  (and  its  hole 
counterpart)  can  be  easily  shown  to  yield  reasonably  sim¬ 
ple  relation  between  Ytir  And  the  pr  products. '*  In  partic¬ 
ular,  assuming  that:  Pft„  -  ppp,  >  kE^^Pnfp  (i  e--  fields 
higher  than  kkgT/q  =  16(X)  V/cm  bul  not  too  high,  see 
below)  the  solution  for  Ya  (see  Eq.  (2)]  becomes'* 

Ya=(lkkgT/qEo){p,tihtntf/[iy;in)^  -  (6) 

If  we  further  assume  that  pp„  >  IpPp,  this  equation  can  be 
further  reduced  to 

Ya=  ( 2kkgT/q )  [phT  p(p,T^„)]{\/Eo).  ( 7 ) 

Since,  in  the  experiment,  2kkgT/q  is  known,  a  presentation 
of  the  measurable  quantity  Ytn  (see  Uq.  (3)]  as  a  function 
of  \/Eo  is  expected  to  yield  the  steady-state  carriers  pr's 
ratio.  On  the  other  hand,  the  known  result  (in  the  Eq 
—  0  limit)  for  the  ambipolar  diffusion  length  L  is  that"  '^ 

Y,„=(\^k^L^)-\  (8) 

where 

L^=(2kgT/q)phT^.  (9) 

As  pointed  out  above,  the  principle  of  the  method  is  to 
substitute  the  experimental  results  for  Yeir Eqs.  (7)-(9) 
in  order  to  derive  the  two  pr  products  without  resorting  to 
the  inaccurate  estimates  of  o  and  G. 

Before  turning  the  experimental  results  let  us  consider 
the  consequences  of  condition  (5).  For  the  typical  grating 
of  A  =  10~*  cm  this  condition  implies  that  Eq  >  1600 
V/cm,  and  for  a-Si:H  if  (Pfd„}/(phdp)  at  10,  one  must  not 
exceed'*  a  light  intensity  of  1.5  mW/cm^  (for  higher  drift 
mobility  ratios  lower  intensities  should  be  used).  For  the 
demonstration  of  the  proposed  method  under  convenient 
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FIC.  I.  The  meetured  dependence  of  the  excee*  carrien  normalized- 
gratmg-amplitude  vz  the  applied  electric  6eld,  under  three  different  back¬ 
ground  He-Ne  illuminations.  The  ambipolar  diffusion  lengths  are  deter¬ 
mined  from  the  zero  field  limit. 

experimental  conditions,  this  constraint  led  us  to  choose  an 
a-Si:H  material  for  which  <  10  as  well  as 

the  lowest  light  intensity  which  enabled  a  convenient  mea¬ 
surement  (0.7  mW/cm^).  A  rough  preliminary  estimate  of 
this  ratio  was  determined  from  the  corresponding  photo¬ 
conductivity  and  the  diffusion  length  measurements  by  ap¬ 
plying  the  method  of  Ritter  et  al.  We  note  in  passing  that 
the  material  used  in  the  following  experiment  had  a  pho¬ 
toconductivity  of  1.7  X  10“’  (ftcm)“'  at  100 
mW/cm*  and  that  for  this  light  intensity  is 

larger  than  the  corresponding  value  for  the  weak  light  in* 
tensity  which  we  have  used  for  the  verification  of  the  linear 
dependence  of  Ytt  on  1/Eq.  This  is,  of  course,  a  result  of  the 
“faster”  displac<nnent  of  the  quasi-Fermi  level  for  elec¬ 
trons,  in  comparison  with  that  of  the  holes,  with  increasing 
illumination  intensity  in  a-Si:H. We  emphasize  that  in 
contrast  to  the  common  —  0  case,”''*  we  are  interested 
here  in  nonambipolar  conditions  for  the  excess  carriers. 

III.  EXPERIMENTAL  RESULTS 

The  experimental  details  of  the  present  work  were 
much  the  same  as  those  described  in  detail  previously.”''’ 
The  measured  dependence  of  Y«g  on  the  applied  electric 
field  £|)  is  shown  in  Fig.  1.  The  conspicuous  qualitative 
feature  of  the  data  is  the  inflection  point  in  the  vs  £o 
dependence,  and  in  particular  its  shift  towards  the  lower 
values  of  Eq  with  decreasing  illumination  intensity.  This 
result  is  consistent  with  the  predictions  from  computer'’ 
and  analytic'*  solutions  of  the  continuity  equations  [see  Eq. 
(4)],  and  with  the  expected  transition  from  ambipolar  to 
nonambipolar  transport  in  the  grating  [Eq.  (3)]  with  de¬ 
creasing  dielectric  response  time  and  inCTeasiiig  field.  The 
transition  (which  is  signaled  by  the  inflection  point)  is 
necessary  for  the  applicability  of  Eq.  (6).  As  discussed 
above,  and  as  can  be  gathered  from  the  data,  this  condition 
is  well  fulfilled  for  illumination  levels  below  2.S  mW.  Also 
given  in  the  figure  are  the  ambipolar  diffbsion  lengths  de¬ 
rived  by  using  the  y^  values  at  -  0  and  relation  (8) 
(with  A  »  9  X  10“’  cm),  llieae  values  have  been  con¬ 


FIQ.  2.  The  results  of  Fig.  I  for  the  He-Ne  illumination  of  0.7 
mW/cffl’,  plotted  as  a  (unction  of  the  inverse  applied  electric  field.  The 
slope  of  the  plot  and  the  zero  field  result  yield  the  /tr  products  for  the  two 
carriers. 

firmed  independently  by  a  more  detailed  PCG  measure¬ 
ment.'’  We  should  point  out  that  the  measurement  for 
Y^  <  O.IS  requires  an  extremely  high  sensitivity  of  the 
technique  which  is  beyond  our  capability  at  present. 

Following  the  observation  of  nonambipolar  transport 
at  2S(X)  V/cm  (the  inflection  point,  for  =  0.7  mW/cm, 
see  Fig.  1 )  we  expect,  for  higher  fields,  the  linear  y^  vs 
l/£^  dependence  predicted  by  our  “diffusion/drift  pertur¬ 
bation”  approach'*  [Eq.  (7)].  Indeed,  this  expectation  is 
fulfilled  by  the  y^  vs  I/Eq  plot,  shown  in  Fig.  2.  Compar¬ 
ing  the  slope  in  the  linear  regime  with  Eq.  (7)  we  get  a 
value  of  0.24  for  the  minority  to  majority  fir  ratio.  Com¬ 
bining  this  with  the  fi/,Tp  result  [derived  as  explained 
above  from  Eqs.  (8)  and  (9)],  we  find  that  the  sum  of  the 
two  fir’s  is  2  X  10  “*  em’/V. 

IV.  DISCUSSION 

In  order  to  check  the  qualitative  validity  of  the  above 
novel  theoretical  and  experimental  approaches,  we  have 
compared  the  above  sum  of  the  fir's  with  the  one  derived 
from  photoconductivity,  although,  as  is  apparent  from  the 
above  discussion,  these  may  differ  by  a  factor  as  large  as  2. 
Our  best  experimental  estimate  of  a  under  the  0.7-mW/ 
cm^  illumination  (electrode  separation  0.4  mm,  beam  di¬ 
ameter  3  mm,  and  sample  thickness  1  ^m)  was  1.3  X  10  “  ’ 
(ncm)“'.  The  corresponding  He-Ne  laser  photon  flux 
was  2.1  X  I  o'’  photons/cm’s.  Assuming  the  commonly 
accepted’*’  30%  reflection,  and  an  absorption  coefficient  of 
2  X  lO’  cm  “ '  we  estimated  a  carrier  bulk  generation  rate 
of  2.8  X  lO”  carrier-pairs/cm’ s.  Using  these  values  in  Eq. 
( 1 ),  we  find  that  ^*7* N  +  MsT"*"  2.9  X  10“*cm’/V. 
Considering  the  more  reliable  procedures  for  the  determi¬ 
nation  of  this  sum  firom  the  measured  photocurrent,''’'’' 
we  found  somewhat  different  estimates.  All  the  estimates 
can  be  summarized,  however,  by  (3  ^  1 )  X  10  “ '  em’/V, 
(which  agrees  with  the  uncertainty  range  suggested  in  the 
above  given  discussion).  The  important  point,  of  course,  is 
the  good  agreement  of  this  value  with  the  value  derived 
from  the  data  of  Fig.  2.  This  agreement  proves  the  validity 
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of  our  suggested  new  method  for  the  derivation  of  the  /xr*s. 
For  the  reasons  given  above,  the  difference  between  the  two 
results  is  attributed  to  the  uncertainties  associated  with  the 
extraction  of  the  ;tr  from  photoconductivity,  rather  than  to 
possible  shortcomings  of  the  present  method. 

In  conclusion,  we  have  presented  the  first  study  of  the 
PCG  in  the  high  electric  field  regime  and  have  shown 
quantitative  agreement  between  experimental  results  and 
the  simplest  analytical  predictions  for  this  regime.  The 
agreement  with  the  corresponding  photoconductivity  dau 
shows  that  the  suggested  technique  can  be  useful  for  the 
determination  of  accurate  values  for  the  two  carriers  /zr 
products. 
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Thus  far  the  many  reports  concerning  the  utilization  of  the  photocarrier  grating  (PCG) 
method  have  assumed  that  ambipolar  transport  takes  place  in  such  a  PCO  when  it  is  imposed 
on  hydrogenated  amorphous  silim  (a-Si:H).  This  assumption,  which  is  decisive  in  the 
interpretation  of  the  experimental  results  in  terms  of  the  ambipolar  diffusion  length,  has  not 
been  tested  thus  far.  In  this  letter  a  corresponding  testing  criterion  is  proposed,  anid'it 
is  demonstrated  that  whenever  ambipolarity  is  lost,  the  PCXJ-derived  diffusion  lengths  may  be 
wrong.  The  finding  that  ambipolarity  is  maintained  in  device  quality  o-Si:H  is  shown  to 
confirm  the  theoretical  suggestion  tlut,  whenever  observed,  the  ambipolarity  in  a>Si:H  is  due 
to  shallow  trapping  effects. 


By  now  the  photocarrier  grating  (PCG)  technique  for 
the  estimation  of  the  minority-carrier  mobility  lifetime 
(/xr)  product  in  amorphous  semiconductors  has  become 
the  most  accepted  and  widely  used  method  for  this 
purpose.''^  The  very  fundamental  assumption  associated 
with  the  interpretation  of  the  corresponding  experimental 
resultt  is  that  ambipolarity  (or  local  charge  neutrality)  is 
maintained  in  the  PCO.  Under  these  conditions  the  phys¬ 
ical  quantity  derived  is  the  important  photoelectronic  pa¬ 
rameter  which  is  known  as  the  ambipolar  difiiision  length 
L.  This  crucial  assumption  has  not  been  confirmed  thus 
far.  In  fact,  it  is  quite  surprising  that  in  spite  of  the  many 
reports''^  in  which  the  results  were  interpreted  in  terms  of 
L  no  one  has  reported  nonambipolar  transport.  In  view  of 
the  importance  of  the  knowledge  of  the  minority-carrier  pr 
(  kI.^)  product  for  the  physics’  and  application*  of  hydro¬ 
genated  amorphous  silicon  (a-Si:H)  and  in  view  of  the 
present  abundance  of  the  PCG  method,  a  critical  exami¬ 
nation  of  this  assumption  appears  necessary. 

The  PCG  technique  is  based  on  a  comparison  of  the  ac 
(due  to  light  chopping)  photocurrent  through  a  photocon¬ 
ductor  under  two  conditions;  when  it  is  subjected  to  spa¬ 
tially  modulated  illumination  and  when  under  the  same 
illumination  intensity  there  is  no  such  modulation.  The 
relation  between  the  corresponding  ratio  of  the  two  ac  pho¬ 
tocurrents,  ff,  and  the  microscopic  transport  and  kinetic 
parameters  of  the  material,  which  are  summarized  (see 
below)  by  the  parameter  Ytt>  <*  given  by'*’ 

ydr»l(l-m2n^)]''*.  (1) 

Here  yo  <*  *n  experimental  quality  factor  which  is  of  the 
order  of  unity'’’  and  y  is  the  well-known  light-intensity 
exponent  of  the  samples’  dc  photoconductivity  oq  (i.e., 
(Tg  oc  GS,  where  Gq  is  the  carrier’s  generation  rate  under  the 
relatively  large  uniform  dc  background  illumiiution). 

The  purpose  of  this  letter  is  to  suggest  a  test  by  which 
one  can  qualitively  identify  the  existence  or  nonexistence  of 
am'oipolar  transport  in  the  PCG  and  to  demonstrate  exper¬ 
imentally  that  one  can  find  ambipolar  as  well  u  nonambi¬ 


polar  transport  in  a-Si:H  materials.  This  test  is  possible 
now  due  to  the  development  of  an  analytic  theory  of  the 
PCG.’ 

Since  the  concept  of  ambipolarity  is  somewhat  differ¬ 
ent  in  the  present  context  than  in  the  classical  theory  of 
semiconductors'"  sve  redefine  it  here  as  the  set  of  condi¬ 
tions  under  which  local  charge  neutrality  prevails  (see  be¬ 
low).  Under  these  conditions  the  PCG  measurement  can 
be  readily  shown'*’  to  yield  the  relation 

ydr=(I+*^I*)''.  (2) 

where  k  =  2ir/A,  A  is  the  grating  period,  and  L  is  defined 

by’" 

L^=(2kgT/q)  (3) 

Here  kgT  is  the  thermal  energy,  q  is  the  (positive)  elec¬ 
tronic  charge,  p,{pi,)  is  the  electrons’  (holes’)  band  mo¬ 
bility,  6„{6p)  is  the  ratio  between  the  free-electron  (hole) 
concentration  no(po)>  the  total  (free  and  trapped)  con¬ 
centration  of  electrons  (holes),  N  +  nQ  [P  +  Pn),  and  tg  is 
the  “common”  two-carrier  recombination  time.’’'’*'’  Un¬ 
der  the  PCG  configuration,  in  the  zero  applied  field  limit, 
the  excess  free-electron  (hole)  concentration  that  forms 
under  the  generation  grating  is  described  by  An 
s  A,cos(kx)  [Ap  B  AACOs(kx)],  where  A,(Aa)  is  the  ex¬ 
cess  electron  (hole)  concentration  amplitude.  Hence  it  is 
obvious  that  local  charge  neutrality  will  be  maintained  if 
Ap/Op  An/0^  Physically  this  condition  means  that  all 
(free  and  trapped)  generated  electrons  are  not  spatially 
separated  from  all  generated  holes.  Since  the  carriers  dif¬ 
fuse  from  regions  of  excess  illumination  and  the  didectk 
response  time  Tg  {  —  e/og,  where  e  b  the  dielectric  con¬ 
stant)  b  finite,  some  separation  will  always  exist  In  the 
PCG  configuration,  whm  the  only  imposed  length  scale  b 
the  grating  period  A,  one  can  say  that  “local”  charge  neu¬ 
trality  prevaib  if  the  two-carrier  separation  in  the  PCG  b 
much  smaller  that  A.  Thb  happens  when  A  b  much  larger 
than  the  dielectric  diffusion  lengths  of  the  two  carriers.  The 
latter  lengths  are  the  dbtances  traveled  by  the  two  carriers 
due  to  diffusion  within  the  time  xg.  Thb  dbtance  for  all  the 
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(ftce  and  tnpped)  dectrons  is  and  for  all  the 

holes  it  is  Here  />,(/>*)  is  the  diAuion  coef- 

fleknt  of  free  ekctrons  (holes).  Hence  the  lequirements 
for  ambipoiarity  in  the  PCO  are  A  >  (i>^,r^)''^al.,^and 
A  > 

This  conclusion,  from  the  above  simple  physical  pic¬ 
ture,  has  also  been  derived  rigorously.’  The  exact  solutions 
found  (which  are  given  explicitly  in  Ref.  9)  yield  the  de¬ 
pendencies  of  A,  and  Ai,oii1^  and  /  where 

/*^/^**  lO(/(€0,)  -I-  ik^D*]/[ao/(6dJ  -I-  I^D^,  (4) 

and  where  Oo  ==  fCMA)  +  MaAi)-  Considering  the  definition 
(4)  it  is  apparent  that  /  is  a  quantity  that  measures  the 
degree  of  ambipoiarity.  Indeed,  when  the  above  ambipo¬ 
iarity  requirements  for  Lti  and  fulfilled,  /  is  re¬ 
duced  to  and  the  ambipoiarity  condition  A, 

(6^0f)Ai,  is  fiilfllled.  Correspondingly  A/,  is  given  by* 

AM=^G,e/ll/tt  +  2k^Dg].  (5) 

With  application  of  Eqs.  (3)  and  (S)  and  the  definition'* 
of  Ytf  as  the  ratio  of  the  excess  photoconductivities  in  the 
presence  and  absence  of  the  grating,  i.e.,  as’*'* 

ydr=9(/^e^e  +  M*4*)/(C|Oo/Co], 

one  can  immediately  see  that  under  the  ambipolar  condi¬ 
tions  Yff  is  reduced  to  the  well-known'*^  relation  given  by 
Eq.  (2). 

If  the  conditions  for  and  Lw  oof  fbifliled  si¬ 
multaneously  we  have  to  introduce /from  Eq.  (4)  into  the 
expressions  for  Af  and  A^  and  then  use  the  definition  of 
Ytt-  For  example,  in  the  extreme  nonambipolar  limit,  i.e., 
when  /  =  Di/D^  Ytx  hss  the  form 

Ya={l+BI^)/{Ck*  +  D/^  +  E),  (6) 

where  B,  C,  D,  E  are  quantities  determined  (through  cum¬ 
bersome  expressions)  by  the  various  material  parameters 
that  characterixe’  A^  and  A/^  In  the  intermediate  case,  i.e., 
when  the  condition  (or  is  reversed  but  the  condition  for 
Lu  is  fulfilled  [so  that  /  s  ~ ']  Yts  l^os  the 

form 

Ya=‘l<^/iak*  +  blc^  +  c)  +  l/{a’k^  +  b').  (7) 

where  again  a,  b,  c,  a',  b'  are  determined  'oy  the  various 
material  parameters.  The  latter  case  describes,  for  example, 
an  »-type  material  for  which  >  A  >  L^-  The  condi¬ 
tion  >  Lu  definitely  holds  in  undoped  and  phosphor- 
ous-doped  a*Si:H  where  (assuming  the  Einstein  relation'") 
one  would  expect  that  >  D/ftp  (since'* 

>  XQ^iifip).  Hence  in  discussing  n-type  a-Si:H  we  do  not 
have  to  use  the  most  general  expression  for  the  parameter 
/  [Eq.  (4)]  and  we  may  consider  the  three  simpler  cases 
that  yidd  the  results  given  by  Eqs.  (2),  (6),  and  (7). 
The  important  point  to  realize  is  that  unlike  the  ambipolar 
case  given  by  Eq.  (2),  the  cases  described  by  Eqs.  (6)  and 
(7)  yield  a  maximum  in  the  vs  dependence.  Hence  a 
presence  of  such  a  maximum  in  the  experimental  data  is  a 
clear  signature  of  nonambipolar  transport.  A  more  detailed 
discussion  of  the  various  coefficients  in  Eqs.  (6)  and  (7) 
will  be  given  elsewhere. 
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FIG.  I.  The  mcaMUcd  dependence  of  the  normalixed  grating  amplitude  . 
on  the  iquared  wave-vector  of  the  photocarrier  grating,  at  obtained  on 
undoped  device  quality  «-Si  Jl. 

Let  US  estimate  now  the  magnitude  of  the  parameter  / 
[Eq.  (4)]  and  the  corresponding  expected  Yttr  vs  depen¬ 
dence.  In  trying  to  quantify  /  one  notes  that  we  do  not 
know  the  values  of  D,  and  Di,  in  o-Si:H  and  even  for  fi,  and 
fih  we  have  only  rough  estimates.  However  the  order  of 
magnitude  of  these  quantities  is  quite  well  established'*  as 
ftf  2:  10  cm*/V  s  and  ~  1  cm^/V  s.  Assuming  the  Ein¬ 
stein  relations  one  can  then  estimate  that  D, 
=  {kgT/q)Ht  s  0.25  cmVs  and  that  Di,  =  0.025  cm^/s. 
Consideration  of  Fermi-Dirac  statistics'*  and  our  knowl¬ 
edge  of  the  band-tail  widths'*  are  not  expected  to  increase 
these  values  by  more  than  a  factor  of  3.  The  photoconduc¬ 
tivity  is  well  known  for  device  quality*  a-Si:H;  for  a  gen¬ 
eration  rate  of  10*'  cm  ~  *  s  " '  ( I -sun  intensities)  it  is  of 
the  order  of  ao=10”*  (Ocm)"'.  Since  for  0-Si:H, 
e  SB  10  " F  m  ~ '  and  since  in  the  typical  grating  geome¬ 
try  As  1  ^m,  we  find  that  under  the  above  conditions  /  of 
Eq.  (4)  is  of  the  order  of 

/=  ( lO*/0p  -I- 10*)/(  10*/e„  -I-  10’).  (8) 

For  obtaining  the  ambipolar  limit,  i.e.,  in  order  that  the 
Teff  (k*)  dependence  given  by  Eq.  (2)  will  be  maintaining, 
we  must  have  9p  <  I  and  6„  <  10  ~  These  two  conditions 
are  definitely  fulfilled  considering  the  above  values  of  fi, 
and  Hh  and  the  corresponding  values  of  the  drift 
mobilities.'* 

Turning  to  the  experimental  results  derived  from  the 
measurements  on  the  PCG  we  note  that  these  results  are 
usually  presented  by  ^(  A)  plots.'**  Here,  however,  in  order 
to  make  a  comparison  with  the  predictions  (2),  (6),  and 
(7),  we  use  the  dependence  of  Ya  Ube  values  of  which  are 
determined  from  the  measured  B  and  relation  ( 1 )]  on  k*. 
In  Fig.  I  we  show  the  predicted  [Eq.  (2)]  and  the  experi¬ 
mentally  observed  Ya  vs  k*  dependence  for  device  quality, 
undoped  a-Si:H.  In  the  present  case  we  found*  that 
Toy'^  SB  0.91  which  is  consistent  with  the  independent  de¬ 
termination  of  y  (  =0.85).  The  excellent  agreement  be¬ 
tween  these  results  and  the  predicted  behavior  of  Eq.  (2), 
as  well  as  the  above-mentioned  expectation  from  Eq.  (8), 
seem  to  confirm  our  suggestion’  that  the  small  values  of 
9^  and  6p  are  the  reason  for  the  ambipolar  behavior  in  the 
PCG  imposed  on  undoped  a-Si:H.  This  proof  is,  however, 
incomplete  since  it  could  be  that  the  experimentally  ob¬ 
served  Eq.  (2)-like  behavior  is  “universal"  and  has  noth¬ 
ing  to  do  with  our  /  parameter.  In  particular,  this  possi- 
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bility  cannot  be  ignored  since  no  other  experimental  Ytt  vs 
dependence  has  been  reported  in  the  literature*'^  thus 
far.  In  order  to  provide  a  decisive  proof  for  the  applicabil¬ 
ity  of  the  theory,’  and  the  meaning  of  /  as  outlined  above, 
one  has  to  show  that  being  guided  by  our  expression  for  / 
[Eq.  (4)]  it  is  possible  to  obtain  experimentally  nonambi- 
polar  behaviors  such  as  those  given  by  Eqs.  (6)  and  (7).  In 
other  words,  one  has  to  show  that  by  upsetting  the  fulfill¬ 
ment  of  the  conditions  for  or  for  one  can  obtain  a 
deviation  from  the  ambipolar  behavior  given  by  Eq.  (2). 

We  have  chosen  to  demonstrate  a  nonambipolar  be¬ 
havior  by  using  a  doped  material.  Our  approach  is  based 
on  the  finding'^  that  with  a  few  ppm  of  phosphorous  dop¬ 
ing  the  drift  mobility  of  the  electrons  increases  while  the 
drift  mobility  of  the  holes  does  not,  i.e.,  d„  increases  while 
dp  does  not.  On  the  other  hand,  it  is  known  that  such 
doping  causes  an  increase  in  the  photoconductivity.  In  or¬ 
der  to  offset  this  increase  we  have  used  thinner  films  (4200 
A  instead  of  2  /xm )  than  those  used  for  the  undoped  u-Si:M 
films. Hence,  for  the  doped  films,  the  first  term  in  the 
numerator  and  the  second  term  in  the  denominator  of  Eqs. 
(4)  and  (8)  become  dominant.  Correspondingly  we  expect 
a  nonambipolar  Ytt  vs  dependence  [as  given  by  Eq.  (7)1. 
i.e.,  an  apparent  peak  in  the  corresponding  plot.  The  re¬ 
sults  shown  in  Fig.  2  confirm  this  expectation.  In  turn,  this 
observation  proves  the  results  of  our  general  theory,’  the 
meaning  of  the  “ambipolarity”  coefficient  /,  and  our  con¬ 
jecture  regarding  the  importance  of  d„  and  dp. 

Another  extreme  nonambipolar  case  may  be  expected 
in  very  thin  films  of  a-Si:H  where  the  high  defect  density" 
(or  the  strong  effect  of  the  surface  recombination  velocity) 
is  expected"  to  yield  low  values  for  Oq  L.  Correspond¬ 
ingly  we  have  chosen  an  undoped  2S0-A  film  in  which  fir 
is  three  orders  of  magnitude  smaller  than  in  our  2-fxm 
undoped  film.  For  this  case  one  would  expect,  on  the  basis 
of  the  above  large  difference  in  the  majority-carrier  fir's, 
the  behavior  characterized  by  Eq.  (6).  Carrying  out  the 
Ytg  vs  measurement  on  these  films  has  confirmed  the 
expectation  for  a  nonambipolar  transport.  This  can  be  seen 
by  the  predicted  presence  of  a  corresponding  peak  in  Fig. 


on  the  iquand-wave  vector  of  the  photoearrier  grating  at  obtained  for  a 
doped  lample  (PHj/SiH,  •  }  x  10  ~  *)  and  a  very  thin  undoped  lample  of 
ii-Si:H.  Note  that  yyo  ■*  ■  constant  (of  order  I )  and  thus  the  qualitative 
contrast  between  the  results  shown  in  this  figure  and  those  shown  in  Fig. 
1  is  apparent. 


2.  Hence  the  guidance  provided  by  our  ambipolarity 
coefficient’  /  is  shown  again  to  have  a  predictive  power  for 
the  ambipolarity  test.  ^ 

In  conclusion,  we  have  shown  theoretically  and  exper¬ 
imentally  that  a  nonmonotonic  decrease  of  the  normalized 
grating  amplitude  with  the  squared  grating-wave  vector  is 
a  signature  of  a  nonambipolar  transport  in  the  PCG.  Our 
procedure  provides,  then,  an  immediate  qualitative  test  for 
the  interpretation  of  the  results  in  terms  of  the  ambipolar 
diffusion  length.  A  quantitative  application  of  the  ambipo¬ 
larity  factor  /  further  confirms  that  the  claim  made  in  the 
many  reports,  in  which  the  value  of  L  was  derived  from 
PCG  measurements  on  undoped  a-Si:H,  is  justified,  and 
that  the  reason  for  this  is  the  high  trapping  rate  in  this 
material. 
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We  hare  reccatiy  ibewa  that  with  the  eaaawaaly  aeaiiahle  laaahitiea  m  Aa|cr  Hcclraa  ■yactreecoyy  meaeate- 
mcBB,  niBMcaat  emta  caa  he  taemad  ia  the  ateaearad  Aafcr  Uae  iateaeitiia.  To  arewaaie  thie  dilTiealty,  a 
aaimaal  rciatioa  has  hoea  darirad  thcorctieally  whcrehy  the  axyeriaNotally  — a  yae  iaiiaaiiiri  caa  anily  be 
eorrectad  so  aa  to  yield  good  oaiaHMi  far  the  irae  iaiiaaltiw.  La.  those  that  woaM  hare  baea  aieawaed  were  tbe 
isoahaioa  MaMy  gaoA  The  ralMity  of  oar  canaeiioa  praeadvc  was  laeeatly  dimeaaBated  for  the  Ugb-eaergy 
KLL  aad  law  laerfy  LMM  Raao  of  Si,  Al  aad  Mg  as  watt  as  for  the  highaimgy  LMM  sad  iow-eacrgy  MNN 
liaas  of  Ca.  Ia  thia  gaper  wa  axtead  these  atadiaa  to  the  highHiaugy  LMM  Kaaa  of  Za,  Gc,  Fc,  Co  aad  NL  We 
ptaaaat  hate  the  iatiiaaic  liaaahapaa  af  thcaa  Vaaa  as  wdl  as  the  Aagar  seaaitiritiee  iclaiirc  to  ailrer,  aieaaared  with 
dttfcrcat  lasaiatioas.  The  eanaetioa  pracadan  apptted  to  the  data  yieMa  the  trae  seasiihiiiea  to  a  good  apptaa* 
haatioa  aad  is  thcrefsia  hiipartaat  fsr  gaaatifleatioa  aad  for  theoretical  calcalaiioas  of  Aagw  yields. 


INTRODUCTION 


When  the  resolution  width  used  in  measurements  of 
Auger  electron  spectroscopy  (AES)  becomes  compara> 
ble  to  or  larger  than  the  intrinsic  Auger  linewidt^  the 
signal  intensity  is  no  longer  a  linear  function  trf 
resolution.'*’  As  a  result,  atomic  conceotratioos  derived 
from  line  intensity  measurements  usually  depend  on  the 
resolution  used,  a  situation  that  is  detrimental  for  quan¬ 
tification.  Since  the  resolutioR  width  is  usually  pro¬ 
portional  to  the  Auger  energy,  this  effect  is  more 
pronounced  in  measuicments  of  high-energy  Auger 
lines.  Contrary  to  the  common  belief  (see  e.g.  Ref.  3) 
that  the  width  of  tbe  Auger  lines  is  of  the  order  of  3-4 
eV,  the  widths  of  the  higllenergy  lines  of  many  elements 
are  below  1  eV,*  a  fact  that  enhances  the  dependence  of 
tbe  measured  conceotratioos  on  the  resolution  used. 

In  order  to  overcome  these  difficulties,  we  have 
derived'*’  a  universal  relation  whereby  the  experimen¬ 
tally  measured  peaic-to-peak  amplitudes  of  the  differen¬ 
tiate  Auger  sign^  referred  to  as  the  ‘derivative 
amplitude',  can  easily  be  corrected  so  as  to  represent 
quite  accurately  in  most  cases  tbe  true  value  of  the 
4ierivative  amplitude.  All  that  is  required  is  to  multifdy 
the  as-measured  amplitude  for  .the  line  of  each  element 
‘ET  by  1  +  Here  HXEl)  is  the  intrinsic 

width  of  the  line  and  is  the  resolution  width  at 
the  line's  energy'*’  (see  bdkrn).  Thus  the  eoneatd  sensi- 
^«n  tvhitive  to  a  standard  line  ‘St*  (such  as  silver) 
is  related  to  the  corresponding  as-messured  sensitivity 
S^by 

(I) 

where  the  correction  factor  F  is  given  by 

1 -I- rHyaywyEi)]’ 

1  +  CH«SlVK«St)]’ 


*  Anther  lo  whom  eermpoedMcn  tbould  bt  addrawnd. 


Strictly  speaking,  our  analysis  is  valid  for  Gaussian  line- 
shapes,  but  it  is  usually  quite  adequate,  with  somewhat 
reduced  accuracy,  for  more  complex  lineshapes  as  well. 

lo  this  paper  we  display  the  intrinsic  lineshapes  of  tbe 
high-energy  LMM  lines  of  Zn,  Ge,  Fe,  Co  and  Ni  as 
obtained  from  extremely  high-resolution  measurements. 
In  addition,  we  present  the  atomic  sensitivities  at  differ¬ 
ent  resolutions.  The  combined  dau  lead  to  the  cor¬ 
rected,  true  sensitivities  of  these  elements. 


EXPERIMENTAL 


Tbe  samples  studied,  together  with  the  silver  standard, 
were  mounted  simultaneously  in  the  vacuum  chamber 
of  a  Physical  Electronics  model  S60  Auger  microprobe. 
Their  surfaces  were  cleaned  by  argon  ion  bombardment 
Two  modes  of  operation  were  used  iii'the  AES  measure¬ 
ments,  the  'ESCA  (electron  spectroscopy  for  chemical 
analysis)  mode'  and  the  ‘Auger  mode'.  In  both  modes  aa 
electron  beam  was  used  to  excite  the  Auger  transitions. 
In  tbe  ESCA  mode  a  reurding  field  is  applied  to  the 
Auger  electrons  so  that  they  enter  the  analyser  at  a  low 
(constant)  energy.  The  resolution  attained  in  this  mode 
was  extremely  high,  amply  sufficient  to  measure  accu¬ 
rately  the  natural  shapes  and  intrinsic  widths  of  all 
lines.  Tbe  intrinsic  wid^  W,  of  each  line  is  given  at  this 
high  resolution  by  the  energy  separation  between  the 
positive  and  negative  peaks  of  the  Auger  derivative. 

Unfortunately,  our  ESCA  mode  is  unsuitable  for 
quantitttive  intensity  measufements  because  of  losses  in 
beam  intensity  intr^uced  by  (he  retarding  field.  Such 
losses  have  a  complicated  dependence  on  energy,  so 
that  a  quantitative  comparison  of  derivative  amfditudes 
for  different  Auger  lines  is  unfeasible.  Tbe  derivative 
amplitudes  were  therefore  measured  in  the  usual  Auger 
mode,  Le.  aritbout  tbe  application  of  the  reurding  field. 
In  addition,  this  is  the  mode  commonly  used  in  quanti- 
utive  analyses  and  is  hence  of  more  interest  The 
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with  1  primary  beam  enerc'  of  5  keV.  The  derisive 
spectra  were  obtained  by  nine*point  computer  differen¬ 
tiation  of  the  measured  Af(£Hpectra.*  The  measure¬ 
ment  steps  were  sufficiently  fine  (either  0.05  or  O.I  eV) 
that  the  ‘differentiation  potential**  was  always  less  th«n 
0.7  of  the  measured  widths  of  all  the  Auger  lines 
studied.  Hence  this  procedure  yielded  accurate  values 
for  the  derivative  spectra*  in  as  far  as  the  differentiation 
procedure  is  concerned. 

T^e  resolution  width  W,  was  determined^  by  differen¬ 
tiating  the  signal  onginating  from  the  elastically  scat¬ 
tered  electrons  at  a  sufficiently  high  primary  beam 
energy  (~2  keV)  and  measuring  the  separation  between 
the  positive  and  negative  peaks.  In  some  measurements 
W,  was  determined  by  measuring  the  actual  width  of  an 
Auger  line  of  narrow  and  known  intrinsic  width.  The 
measured  width  IV.  is  related  to  W;  and  W,  by  ^ 


(3) 


Obviously  in  the  limit  of  infinitely  good 

resolution  (H'  <•!♦(),  while  for  an  infinitesimally  narrow 
line  <  W,)  =  W,.  Since  W,  is  proportional  to 

energy,  its  value  at  any  Auger  line  can  readily  be  evalu¬ 
ated. 


RESULTS 


Figwr*  1.  Derivative  eoactre  of  the  Zn  LMM  line  in  the  ESCA 
mode  (bonom  curve)  ertd  in  the  Auffw  mode  at  three  different 
resotutione  (three  upper  curvet). 


resolution  was  varied  by  adjusting  the  analyser  aper¬ 
tures.  In  this  manner  three  resolutions  were  obtained, 
around  0.3%,  0.6%  and  I.3*/4  of  the  Auger  energy  E. 
Care  was  taken  at  each  setting  to  perform  the  measure¬ 
ments  on  the  different  samples  under  identical  experi¬ 
mental  conditions.  The  measurements  were  carried  out 


The  silver  line,  used  as  the  standard,  and  ail  the  high- 
energy  LMM  lines  studied  were  measured  both  in  the 
ESCA  mode  (resolution  width  W,  either  60  or  120  meV) 
and  in  the  Auger  mode  at  three  different  resolutions 
(0.33%,  0.55%  and  1.3%  of  the  line  energy).  The 
detailed  linesbape  of  the  silver  line  has  been  reported 
previously.*  Figure  1  displays  the  evolution  of  the  Zn 
LMM  lineshape  as  the  resolution  is  degraded  from  that 
in  the  ESCA  mode  (bottom  curve,  IV,  *  60  meV)  to 
those  of  the  Auger  mode  (three  upper  curves).  The 
curves  in  the  figure  were  shifted  along  the  >>-axis  and  the 
line  intensities  were  normalized  to  10.  As  seen  in  the 
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Figwra  4.  Th«  Co  LMM  lino  mooturod  in  the  ESCA  modo. 
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bottom  curve  of  Fig.  1,  the  zinc  line  is  intrinsically  very 
narrow  (0.6S  eV)  and  has  a  satellite  about  3  eV  above 
the  main  line.  Tlie  utellite  is  still  resolved  (with  lower 
amplitude)  in  the  first  upper  curve  {W,  ~  3.3  eV),  our 
best  resolution  in  the'  Auger  mode,  but  is  completely 
washed  out  in  the  much  poorer  resolutions  of  the  two 
uppermost  curves. 

Figures  2-S  display  the  intrinsic  lineshapes  of  the 
high-energy  LMM  transitions  of  the  rest  of  the  elements 
studied,  as  measured  in  the  high-resolution  ESCA 
mode.  The  line  intensities  in  these  figures  were  normal¬ 
ized  to  full  scale.  We  note  that  the  Ge  LMM  line  (Fig. 
2),  just  like  the  Zn  LMM  line  (Fig.  1),  is  intrinsically 
very  narrow  (1  eV).  Such  linewidths  are  characteristic 
also  of  the  KLL  lines  reported  previously.^  In  contrast, 
the  LMM  lines  of  the  three  transitidn  metals  are  appre¬ 
ciably  wider,  as  can  be  seen  in  Figs  3-S. 

Measurements  in  the  Auger  mode,  such  as  those  dis¬ 
played  in  Fig.  1,  together  with  those  of  the  silver  line 
(not  shown  here)  were  used  to  determine  the  derivative 
amplitudes  of  the  different  lines  and,  through  their 
ratios,  the  atomic  sensitivities.  The  as-measured  sensiti¬ 
vities  of  the  Zn  and  Ge  lines  relative  to  Ag  at  three 
resolutions  are  depicted  by  the  open  symbols  4n  the 
semilogarithmtc  plots  of  Fig.  6.  As  expected,  the  as- 
measured  sensitivities  decrease  as  the  resolution  is 
degraded.  The  full  symbols  represent  the  corrected  sen¬ 
sitivities  Smt,  obtained  by  the  use  of  Eqns  (1)  and  (2) 
with  M|(Ag) «  1.3  eV.*  These  are  seen  to  be  constant  to 
within  10%-20%,  independently  of  resolution.  The 
horizontal  lines  depict  the  averages  of  the  corrected 
sesitivities  for  the  two  elements.  The  curves  below  rep¬ 
resent  the  calculated  dependences  of  on  resolution, 
as  evaluated  from  Eqns  (1)  and  (2),  and  are  seen  to  fit 
well  the  measured  points.  These  curves  immediately 
indicate  what  values  of  s'*  to  be  expected  at  any 
resolution.  For  the  case  at  hand  one  would  need  a 
resolution  of  0.01%  or  so  in  order  to  measure  the  true 
sensitivities.  Obviously  such  resolutions  are  unat¬ 
tainable  by  far  in  AES  measurements.  Hence  one  has  no 
recourse  but  to  use  the  correction  method  described 
here  if  one  wishes  to  arrive  at  reasonable  estimates  for 
the  true  sensitivities. 

Similar  plots  are  displayed  in  Fig.  7  for  the  Ni  and  Fe 
lines.  The  correction  facton  for  Ni  are  smaller  than 


Resolution  {%) 


Figure  7.  As-nwaturad  and  eonactad  tanaitivities  of  tha  LMM 
linat  of  Pa  and  Ni  ralativa  to  the  MNN  lina  of  Ag  at  funetioni  of 
inttrumantal  ratolution. 


those  for  Zn  and  Ge  owing  to  the  larger  width  of  the  Ni 
line  (1.7S  eV).  The  agreement  with  theory  in  this  case  is 
not  very  good.  This  is  not  surprising  since  the  Ni  line  is 
far  from  Gaussian  (see  Fig.  S).  Hence  the  correction 
procedure  yields  only  a  rough  estimate  of  the  true  sensi¬ 
tivity,  yet  the  divergence  in  the  corrected  sensitivities 
does  not  exceed  20%-30%.  The  Co  line  is  also  non- 
Gaussian  (see  Fig.  4)  and  the  agreement  of  the  sensi¬ 
tivity  values  with  theory  is  similarly  poor.  As  for  the  Fe 
line,  the  agreement  with  theory  is  quite  good,  as 
expected  in  view  of  its  Gaussian  iineshape  (see  Fig.  3).  It 
should  be  noted  that  here  the  corrected  sensitivities  are 
lower  than  the  as-measured  sensitivities.  This  is  because 
the  width  of  the  Fe  line  (3.6  eV)  is  much  larger  than  the 
width  of  the  Ag  line  (1.3  eV),  so  that  the  correction  of 
the  silver  intensity  is  larger  than  that  of  the  Fe  line 
intensity  (see  Eqns  (1)  and  (2)). 

The  main  results  for  Zn.  Ge.  Fe.  Co  and  Ni  are  sum¬ 
marized  in  Table  1.  The  table  lists  the  intrinsic  widths 
W^  of  the  lines,  as  obtained  from  Figs  1-S,  the  as- 
measured  sensitivities  relative  to  Ag  at  three  resolutions 
and  the  average  corrected  sensitivities  <5cw,>. 


Figured.  As-measurad  and  cowaetad  Mmitivitiat  of  tha  LMM 
Hnaa  of  Zn  and  Oa  tafaiMa  to  tha  MNN  lina  of  Ag  at  functions  of 
inatrumanial  raaotution. 


TaMc  1.  Aa-tBcaaared  acaaitivitMS  at  three  readatiaaa 
aad  arerage  corrected  acashiritica  <5^,)  ^  Ihv 
LMM  Haaa  af  Za,  Gc,  Ft;  Co  aad  Ni  rciaiivc  la  the 
MNN  Hue  of  Ag.  The  data  were  takca  with  a  gihaary 
haam  caergy  £,  of  5  keV.  Abo  bcladid  are  the  Hat 
eaargita  £  aad  the  iairiaaic  liaewidtha  IF, 
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0.14 
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0.13 

Co 
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4.4* 

0.28 

0.46 

0.46 

OJ* 

Ni 

848 

1.78 

0.41 

0.42 

0.39 

1» 

*Thit  valua  taptaatntt  tha  anorgy  aaoaration  baiwoan  tha  pothiva 
peak  and  tha  moan  of  tha  two  nagativo  paakt  of  tha  dlffarantiatad 
Augar  signal  (too  Fig.  4). 

*Thota  valuot  ara  rough  ottimaiot  (too  last). 


HIGH-ENERGY  LMM  AUGER  TRANSITIONS 


DISCUSSION 


The  use  of  the  high-resolution  ESCA  mode  with  fine 
steps  in  the  dau  accumulation  provides  a  detailed 
P’cture  of  the  natural  lineshapes  of  all  the  Auger  lines 
studied  All  the  lines,  except  the  Fe  line,'  exhibit  double 
peaks  or  prominent  satellites  in  their  close  vicinity. 
Such  fine  structure  should  be  useful  in  studying  the  dif¬ 
ferent  initial  and  final  states  involved  in  the  Auger  tran¬ 
sitions  (see  e.g.  Ref.  7). 

The  narrow  widths  of  the  high-energy  LMM  lines  of 
Zn  and  Ge  (1  eV  or  less),  just  like  those  of  the  high- 
energy  KLL  lines  reported  previously,*  arc  not  sur¬ 
prising  since  these  lines  originate  from  deep  level 
I  transitions.  As  a  result  of  the  narrow  widths,  the  errors 
in  the  as/tneasured  sensitivities  are  quite  large,  requiring 
correspondingly  large  correction  factors.  The  consider- 
ably  larger  widths  of  the  high-energy  lines  of  Fe.  Co 
and  Ni  are  probably  due  to  the  unfUl^  d-shells  in  these 
elements. 


I  PUASE  ^ 

The  observation  that^he  as-measured  sensitivities  at 
different  resolutions  can  vary  whereas  the  corrected  sen¬ 
sitivities  are  the  same  to  within  10%-20%  demonstrates 
once  again  the  validity  of  our  simple  correction 
procedure — this  much  in  spite  of  the  fact  that  the  Auger 
lines  studied  were  not  stnctly  Gaussian,  some  deviating 
quite  significantly  from  a  Gaussian  lineshape. 

The  corrected  sensitivities  provide  good  estimates  for 
the  true  sensitivities.  They  can  differ  by  one  to  two 
orders  of  magnitude  from  the  as-measured  values.  It  is 
these  true  sensitivities  that  the  results  of  any  theoretical 
calculations*'*  of  Auger  yields  should  be  compared 
with.  This  is  particularly  important  for  Auger  lines  orig¬ 
inating  from  different  transitions.* 
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Until  the  recent  development  of  the  photocarrier  grating  (PCG) 
method  it  was  impossible  to  follow  the  phototransport  properties  associated 
with  minority  carriers  in  materials  in  which  the  diffusion  length  of  these 
carriers  is  in  the  submicron  range.  The  two  systems  in  which  significant 
knowledge  of  the  minority  carrier  properties  have  already  been  extracted  by 
the  application  of  the  PcG  are  amorphous  silicon  and  polycrystalline  Cu— 
ternary— chalcopyrites.  While  the  basic  physics  of  the  method  is  well 
understcxxl  the  interpretation  of  its  experimental  results,  in  terms  of  the 
microscopic  transport  and  recombination  mechanisms,  is  not  settled.  In  this 
presentation  we  give  a  short  review  of  the  present  understanding  of  the 
PCG  methcxl  and  its  application  to  the  above  materials. 


I.  INTRODUCTION 

The  special  feature  which  singles  out  semiconductors  from  other 
materials  is  that  the  electrical  conduction  in  those  materials  is  carried  out 
by  two  types  of  carriers.  The  carrier  which  dominates  this  conduction  is 
known  as  the  majority  carrier.  On  the  other  hand  the  other,  minority, 
carrier  is  the  carrier  wnich  controls  the  electrical  behavior  of  inhomogeneous 
semiconductors  in  general,  and  in  semiconductor  devices  in  particular. 
Hence,  in  addition  to  the  interest  in  the  basic  physics  of  the  electronic 
structure  and  the  transport  mechanism  associated  with  the  minority  earner, 
there  is  a  strong  motivation  to  study  the  properties  of  this  carrier  from  the 
device  physics  point  of  view.  Correspondingly  an  intensive  effort  was  madp 
in  order  to  characterize  the  minority  carrier  properties  in  semiconductors  . 
Of  these  properties  the  minority  carrier’s  diffusion  length,  Lj^,  was  studi^ 
more  than  any  other  parameter.  yields  information  on  both,  the 

transport  and  the  kinetics  of  the  minority  carrier,  and  is  the  dominant 
parameter  in  determining  the  performance  of  the  well  understood  and 
widely  used  bipolar  devices. 

It  turns  out  that  the  many  methods  which  were  applied  to  crystalline 
semiconductors  have  used  in  one  way  or  another  the  fact  that  Ljjj  ^  10 
jun.  The  development  of  new,  in  particular  thin  film,  semiconductors  where 
Lj_  <  1  pm  has  challenged  the  researchers  of  those  materials  to  develop 
suiUble  method  for  the  measurement  of  their  L^j.  TTie  first  attempt^ 
which  were  based  on  methods  used  for  crystalline  semiconductors,  failed 
mainly  because  in  the  corresponding  inhomogeneous  configvrations  the  short 
Lju  could  not  be  separated  ftonL  the  drift  length  contribution  due  to  the 
presence  of  internal  electric  fields^. 
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A  breakthrough  in  the  possibility  of  the  measurement  of  came 
when  Ritter,  Zhidov  and  Weiser*  suggested  the  use  of  a  photocarrier 
sinuosidal  grating  which  is  based  on  the  possibiliw  of  comparine  L_  with  a 
variable  carrier— concentration  grating  period  A.  It  is  obvious  then  that  for 
Ljjj  >>  A  the  existance  of  the  grating  does  not  effect  the  measured 
photoconductivity,  or-,  while  otherwise  a-  is  effected  by  the  presence  of  the 
grating.  The  relatitm  between  cTg  and  A  yields*^  then  the  value  of  Lj^.  It 
took  however  three  years,  aftCT  the  work  of  Ritter  ct  to  prove 

experimentally  that  indeed  the  length  parameter  determmed  from  the 
measurement,  L,  is  associated  with  the  minority  carriers'*.  The  exact 
relation  between  L  and  is  still  however  a  matter  of  intensive  discussion. 
On  the  other  hand  the  latter  experimental  proof  enabled  the  utilization  of 
the  PCG  technique  to  resolve  many  material -physics  and  device -physics 
issues  in  Iwdrogenated  amorphous  silicon**  (a— Si:H)  and  Cu— ternary 
chalcopyritcs^(CuDiSe2  and  CuGaSe2). 


IL  THEORfrriCAL  DEVELOPMENTS 

The  basic  physical  picture  of  the  photocarrier  grating  formed  by  a 
small  sinuosoidal  carrier  generation  grating,  which  is  superimposed  oh  a 
much  larger  uniform  generation  background,  is  now  well  understood^. 
Quantitatively,  it  is  generally  agreed  that  the  length  L  derived  from  the 
PCG  measurement  can  be  written  as  . 

L  =  V2Dt  (1) 

where  D  is  some  "effective"  diffusion  coefficient  and  t  is  some  "effective" 
small  signaL  recombination  time  that  can  be  written  (for  high  enough  ct- 
values)  as^’®: 


lA  =(8U/an)p  +  (8U/ap)n  ,  (2) 

where  U  is  the  common,  two  carrier,  recombination  rate  in  the  uniform 
background  of  the  electron  (n)  and  hole  (p)  concentrations.  Until  recently 
even  Eq.  (1)  was  in  doubt  since  it  is  not  obvious  that  a  small  superimposed 
sinuosidal  carrier  generation  will  yield  a  small  superimposed  sinusoidal 
carrier  concentration  .  The  questions  of  present  debate are  the  meanings 
of  an  "effective"  D  and  an  "effective"  t.  These  concepts  depend  of  course 
on  the  set  of  assumptions  made  in  the  theoretical  calculations,  and  their 
justification,  when  applied  to  a  given  material  and  under  given  conditions. 
We  list  here  the  set  of  assumptions  made  in  order  to  get  the  simplest 
expression  for  L  in  terms  of  the  microscopic  parameters  of  the  material,  and 
we  refer  to  objections”''^  made  regarding  this  set 

In  the  value  of  D  the  question  is  whether  one  should  considu  only 
shallow  trappi^  or  whether  deep  trapping  should  also  be  included  \  *^0 
other  Quesuon^is  whether  the  quasi  Fermi  level  (or  the  carrier  generation 
rate)  enects  the  value  of  D  and  whether  one  can  assume  the  ful^ment  of 
the  Einstein  relation  in  systems  studied  thus  far  by  the  PCG*^.  The 
expression  for  t  is  even  more  problematic  since  one  has  to  assume  a  priori 
a  recombination  mechanism  and  thus  an  explicit  expression  for  U.  Also,  as 
with  the  quantity  D,  the  question  of  which  trapping  should  be  considered  is 
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still  opcn^  Correspondin^y  the  interpretation  of  the  value  of  L  and  its 
light  intensity  and  temperature  dependences,  in  terms  of  electronic  structure 
and  transport  mechanism,  is  determined  by  the  assumptions  made  regarding 
the  material  under  study.  The  simplest  assumptions  are  the  fulfillment  of 
the  Einstein  relation,  constant  carrier  mobilities,  only,  shallow  (or  no) 
trapping,  and  a  Shockley  Read  recombination  mechanism-^'®  (single  level,  n, 
p  <<  Pp  Np  where  Pj.  and  Nj.  are  the  correspondine  concentrations  of 
available  recombination  centers).  These  assumptions  yield  the  expression^; 

=  (2kT/q)[M.gM.ji'rnTp]  /  +  p.jjTpJ,  (3) 

where,  kT  is  the  thermal  enerw,  q  is  the  electronic  charge,  p,g(p.j^)  is  the 
electron  (hole)  mobility  and  T„(Tn)'  is  the  electron  (hole)  recombination  time 
in  the  uni  form  background  or  the  two  carriers.  One  notes  of  course  that  if 
the  minority  .carriers  are  the  holes  (p«h"n  l^-e'n)  expression  is 

reduced  to  =  (2kT/q)p.jjTp  and  thus  tne^assumption: 

L  =  Ito  (4) 

is  fuJl^ed.  As  stated  out  above  there  is  cuirentlv  sufficient  experimental 
data^’^  to  assume  that  Eq.  (4)  is  a  fairly  good  approximation  for  the 
materials  studied  thus  far.  However  the  acrnracy  of  this  relation  is  definitely 
an  open  question  and  so  is  the  phyrical  information  that  can  be  extracted 
regarding  the  recombination  process 


III.  EXPERIMENTAL  DEVELOPMENTS 

At  present  there  are  some  twenty  PCG  experimental  systems  in  the 
world,  which  except  for  some  minor  details  arc  similar  to  the  ^tem  used 
initially^.  The  photocarricr  generation  grating  is  created  by  the  interference 
of  two  coherent  laser  beams,  which  are  of  very  different  li^t  intensities. 
Hence  the  generation  of  a  large  uniform  carrier  background  and  a  sraa^l 
superimposed,  sinusoidal,  carrier  grating  which  has  a  submicron  period'^’  . 
The  initial  experiment^  results  obtained  by  the  •  PCG  method  yielded 
observations  which  appear*  to  be  contrary  to  expectations ^  The  variation  of 
Lr  with  the  variation  of  materials  was  found  to  be  over  a  narrower  range 
than  the  corresponding  change  in  the  majority  carrier  p.T  product,  the  r.atio 
p,uT_/p.gTn  was  found  to  be  much  larger  than  the  one  observed  in  time— 
or— fiignt  measincments,  and  the  light  intensity  exponents  ^cre  found  to  be 
different  for  L^and  p-.T-  in  contrast  with  the  eroectations"  from  the  simple 
Shockley  Read  model,  it  turns  out"*^  that  the  mst  observation  is  simply  a 
result  of  the  A  range  available  at  present  in  PCG  experimental  systems, 
while  the  second  observation  relates  only  to  the  oldest  sets  of  data,  which 
were  taken  in  measurements  carried  out  on  nearly  compensated  materials. 
The  third  observation  is  obviously  a  result  of  the  application  of  the  above 
Shockley  Read  approach.  One  can  show  that  if  charge  neutrality  is 
considered  a  one  level  Shockley  Read— like  model  of  recombination  can  still 
explain  most  of  the  data  reported.  This  is  since  under  this  consideration 
(2)  takes  the  two  term  form: 

lA  =  Kjn  +  K2(Nj.+Pr) 


(5) 
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where  Kj[  ana  K2  are  constants.  Hence  the  range  of  possibilities  covered  by 
(5)  IS  enou^  to  account  for  the  data.  In  fact  for  materials  where  this 
is  not  the  case  it  is  enough  to  introduce  just  one  additional  (e.g.  correlated) 
recombination  center  in  order  to  explain  the  data  .  At  present  the  proper 
approach  seems  to  be  the  one  in  which  the  minimum  assumptions,  necessary 
to  explain  the  deviations  from  the  behavior  given  by  Eq.  (3),  are  made. 

Among  the  achievements  of  the  applications  of  the  PCG  technique  to 
a-Si;H  materials  we  should  mention^;  optimization  of  materials  by  nnding 
deposition  conditions  under  which  L  is  maximized,  determination  of  the 
majority  carrier  to  minority  carrier  transition  (as  a  function  of  boron 
doping),  evidence  for  the  existence  of  safe  hole  traps,  the  determination  of 
the  dependence  of  the  minority  carrier  quasi  Fermi  level  on  temperature, 
the  surface  recombination  velocity,  the  relation  between  and  the 

corresponding  quantity  in  time— of— flight  measurements,  and  the 
experimental  proof  that  it  is  the  minority  carrier  which  is  the  limitiM 
carrier  in  the  p— i— n  solar  cells.  In  the  Cu- ternary- chalcopyrites  the  PCG 
yielded  the  value  of  and  its  relation  to  the  variable  crystallite  size  in  the 
material.  Here  again  conditions  were  found  for  maximizing  as  a 
function  of  composition  and  deposition  parameters^. 

In  conclusion,  the  PCG  is  presently  the  only  available  method  for  the 
experimental  study  of  the  phototransport  properties  of  the  minority  carriers 
in  materials  where  the  corresponding  diffusion  length  is  in  the  submicron 
regime.  The  present  understanding  of  the  method  and  its  results  is 
definitely  sufficient  for  the  understanding  of  the  device  physics  of  the 
corresponding  materials.  From  the  progress  made  thus  far  it  appears  also 
that  additional  data  and  proper  theories  will  enable  to  make  this  method  an 
effective  tool  for  the  derivation  of  more  detailed  information  regarding  the 
basic  physics  of  these  materials. 
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Lineshapes  and  yields  of  high-energy  Auger  lines 
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Abstract  We  present  high-resolution  measurements  of  the  Auger  transitions  of  Au, 

Pt  and  W.  AH  the  iintt  studied  are  narrower  than  2  eV.  We  also  measured  the  atomic  sensitivities 
for  the  above  Auger  transitions  with  a  resolution  of  0.4%  of  the  Auger  energy  (somewhat  better 
than  commonly  available).  With  this  resolution  it  is  impossible  to  resolve  fine  structure  and 
gross  errors  are  incurred  in  the  measurements  of  Auger  signals.  We  corrected  the  measured 
sensitivities  with  the  method  developed  previously;  the  correction  factors  for  the  sensitivities  are 
around  ten.  These  large  factors  represent  also  the  errors  incurred  in  the  as-measured  sensitivities 
if  no  correction  is  applied. 


With  the  commonly  available  analyser  resolution  in  Auger  electron  spectroscopy  (AES),  of 
about  0.5%  of  the  Auger  energy,  the  lineshs^s  can  be  greatly  distorted  and  most  of  the  fine 
structure  obliterated.  Large  errors  can  also  be  incurred  in  the  measured  signal  intensities 
[1].  This  effect  leads  to  gross  errors  in  the  measured  atomic  sensitivities  and  hence  in  the 
evaluation  of  cross  sections  for  Auger  excitation.  In  order  to  overcome  these  difficulties, 
we  have  derived  [1]  theoretically  a  universal  relation  whereby  the  experimentally  measured 
peak-to-peak  amplitudes  of  the  differentiated  Auger  signal  can  be  corrected  so  as  to  represent 
quite  accurately  the  true  atomic  sensitivities. 

The  samples  studied  were  mounted  together  in  the  vacuum  chamber  of  a  Physical 
Electronics  model  560  Auger  microprobe.  T\vo  modes  of  operation  were  used  in  the  AES 
measurements  [1],  the  *ESCA  mode’  and  the  ’Auger  mode’.  In  both  modes  an  electron  beam 
was  used  to  excite  the  Auger  transitions.  In  the  E5CA  mode,  a  retarding  field  is  applied  to 
the  Auger  electrons  so  that  they  enter  the  analyser  at  a  low  (constant)  energy.  The  resolution 
attained  in  this  mode,  200  meV,  was  very  good.  Such  a  resolution,  combined  with  a  fine 
mesh  (0.05  eV),  was  amply  sufficient  to  measure  accurately  the  natural  shapes  and  intrinsic 
widths  of  all  lines.  Unfortunately,  our  ESCA  mode  is  unsuitable  for  quantitative  intensity 
measurements  [1].  The  Auger  signals  were  therefore  also  measured  in  the  usual  Auger 
mode,  with  a  resolution  of  0.4%  of  the  Auger  energy.  We  present  results  on  the  highest- 
intensity  MNN  lines  of  Au,  Pt  and  W,  corresponding  to  the  MsN67N67  transitions.  In  the 
ESCA  mode  the  lines  are  well  resolved  and,  in  addition  to  the  main  peaks,  several  satellites 
were  observed.  The  widths  of  the  main  peaks  and  the  satellites  are  less  than  2  eV  (see 
table  1).  In  order  to  obtain  the  atomic  sensitivities  (with  respect  to  Ag),  we  measured  the 
Auger  signals  of  Au,  Pt,  W  and  Ag  in  the  Auger  mode.  The  lines  appeared  now  as  single 
peaks  with  no  discernible  fine  structure.  Their  as-measured  atomic  sensitivities  (relative  to 
Ag),  Sfoem*  are  listed  in  table  1. 

The  high-energy  MNN  lines  of  Au,  Pt  and  W  are  quite  narrow,  less  than  2  eV.  Because 
of  their  high  energies  ('^  1700-2000  eV),  the  commonly  available  analyser  resolution  (about 
0.5%)  in  their  vicinity  would  be  8-10  eV.  Obviously,  it  is  impossible  to  resolve  such 
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IkUe  1.  The  intrinsic  linewidths  of  the  MsNc7Ne7  transitions  of  An;  Pt.  and  W,  and  their  as- 
measured  and  corrected  atomic  sensitivities  relative  to  Ag.  The  piimaiy  electron  beam  energy 
was  S  keV  and  the  analyser  resolution  0.4%.  The  effective  width  of  the  Ag  line  used  in  the 
H  calculations  of  Scon  was  1.3  eV  [I].  Also  included  are  the  theoretical  values  of  Mrockowsld 

and  Lichtman  {2]  and  those  listed  in  [3]. 


Element 

Ei  (eV) 

Wi  (eV) 

Scow 

Theory  [2] 

Handbook  [3] 

Au 

2015 

1.95 

0.033 

0.27 

0.0V 

0.036 

Pt 

I960 

1.90 

0.036 

0.30 

0.092 

0.042 

W 

1729 

1.45 

0.076 

0.83 

0.130 

0.080 

narrow  lines  with  that  kind  of  resolution.  Thus,  the  only  way  to  obtain  the  intrinsic  lineshape 

•  of  the  lines  is  by  operating  the  analyser  in  a  very-high-resolution  model,  such  as  the  ESCA 
mode  used  here.  The  narrowness  of  the  high-energy  lines  and  the  poor  analyser  resolution 
in  the  Auger  mode  pose  another  problem,  that  of  how  to  obtain  the  true  atomic  sensitivities 
associated  with  the  transitions.  Our  analysis  makes  it  possible  to  overcome  this  difficulty 
even  when  the  resolution  is  much  poorer.  All  that  is  needed  for  a  good  approximation  of 

^  the  true  sensitivity  is  to  multiply  the  as-measured  line  intensities  by  the  correction  factor 
[1]  [1  +  where  Wt  is  the  analyser  resolution  width  at  the  Auger-line  energy  and 

Wi  is  the  natural  Auger  linewidth.  This  procedure  has  to  be  applied  for  the  transition  in 
question  as  well  as  for  the  standard  transition  (Ag  in  this  case).  The  atomic  sensitivities 

^  are  then  given  by  the  corrected  line  intensities  relative  to  that  of  Ag.  The  fifth  column  in 
table  1,  5coit«  lists  the  values  of  the  corrected  sensitivites.  We  note  that  they  are  about  an 
order  of  magnitude  higher  than  the  as-measured  sensitivities  (Smas)-  They  are  also  higher 
than  the  corresponding  theoretical  values  derived  hx)m  Mrockowski  and  Lichtman’ s  [2] 
calculations  and  listed  in  the  sixth  column  of  the  table  (theory).  For  comparison  purposes, 

•  we  include  in  the  last  column  the  sensitivities  taken  from  [3].  They  are  fairly  close  to  our 
as-measured  sensitivities  obtained  with  comparable  resolution.  We  wish  to  stress  that  the 
corrected  values  represent  true  sensitivities.  It  is  these  values  that  should  be  compared  with 
theoretical  calculations  on  the  one  hand,  and  be  used  for  more  accurate  quantification  of 

^  Auger  measurements  on  the  other. 
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ABSTRACT 

Pulsed  measurements  on  the  solid-electrolyte  system,  which  proved  vety 
useful  in  the  study  of  ctystalline  semiconductors,  have  been  found  to  be  equally 
effective  when  applied  to  hydrogenated  amorphous  Si  films.  Here,  as  well,  the 
aSi:H/electrolyte  interface  is  essentially  blocking  to  current  flow  and,  as  a  result, 
surface  space-charge  layers,  ranging  from  large  depletion  to  very  strong 
accumulation  conditions,  can  be  induced  and  studied.  In  particular,  valuable 
information  can  be  gained  on  the  density  of  the  localized  bulk  states. 
Measurements  in  the  depletion  range  under  illumination  yield  directly  the  total 
density  of  occupied  states  in  the  entire  energy  gap.  Tliis  is  very  useful  in 
obtaining  a  quick  and  reliable  assessment  of  the  quality  of  the  amorphous  films. 
In  high-erade  films  we  find  that  the  total  density  of  occupied  states  is  around 
10**  cm'^.  The  data  in  the  accumulation  range,  on  the  other  hand,  provide 
useful  information  on  unoccupied  states  near  the  conduction  band  edge.  The 
blocking  nature  of  the  amorphous  Si/electroiyte  interface  is  utilized  also  to  apply 
a  sweep-out  technique  for  an  accurate  determination  of  ^.t,  the  product  of  the 
electron  mobility  and  lifetime,  even  when  this  value  is  verv  low.  In  a  rather 
poor-quality  film,  for  example,  we  find  p,T  to  be  5  x  10**  cmW. 

INTRODUCTION 

Considerable  effort,  both  theoretical  and  experimental,  has  been  devoted  to 
derive  the  density  of  states  spectra  in  a-Si;H  films  [1],  In  this  paper  use  is 
made  for  this  purpose  of  pulse  measurements  on  the  semiconductor/electrolyte 
(S/E)  system.  Such  measurements,  which  proved  to  be  very  useful  in  the  study 
of  crystalline  semiconductors  [2-4],  have  been  found  to  be  equally  effective  when 
applied  to  a-Si:H  films.  The  essentially  blocking  nature  of  the  S/E  interface 
allows  one  to  induce  by  an  applied  bias  space-charge  layers  at  the  a-Si:H 
surface,  ranging  from  large  depletion  to  very  strong  accumulation  conditions.  In 
this  manner,  the  entire  energy  gap  in  the  space  charge  region,  together  with  its 
localized  states,  can  be  swung  below  and  above  the  Fermi  level.  At  the  same 
time,  one  can  measure  the  surface  space-charge  density  as  a  function  of  the 
barrier  height  V^.  In  a-Si:H,  resides  predominantly  in  the  localized  states 
(except  in  strong  accumulation  conditions),  so  that  the  measurements  yield,  at 
least  in  principle,  the  energy  distribution  of  the  density  of  states. 

High-grade,  device-quality  films  are  usually  close  to  intrinsic  (resistivity 
"-10**  ohm-cm).  Such  high-resistivity  films  cannot  be  handled  by  our 
measurement  technique  because  the  surface  space-charge  capacitance  cannot 
be  charged  within  the  short-duration  applied  pulse  biases.  Accordingly,  most  of 
the  results  to  be  reported  here  have  been  obtained  under  illumination,  for 
which  the  photo  resistivity  is  t3rpically  10*  ohm-cm  or  less.  Our  measurements 
in  the  depletion  range  are  in  remarkably  good  agreement  with  the  theoretical 
curve  that  takes  into  account  the  presence  of  localized  state.s  [5].  The  data  thus 
yield  directly  the  total  density  of  occupied  states  in  the  entire  energy  gap.  This 
is  very  useful  in  obtaining  a  quick  and  reliable  assessment  of  the  quality  of  the 
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amorphous  films.  In  high-erade  materials  we  find  that  the  total  density  of 
occupied  states  is  around  10''  cm*'.  The  results  in  the  accumulation  range,  on 
the  other  hand,  provide  useful  information  on  the  unoccupied  states  near  the 
conduction-band  edge.  Analysis  of  the  data  in  this  case,  however,  is  more 
difficult  because  surface  states  are  apparently  also  involved. 

The  blocking  nature  of  the  a-Si:H/eiectrolyte  interface  is  utilized  also  to 
apply  a  sweep-out  technique  for  an  accurate  determination  of  fvr,  the  product  of 
the  electron  mobility  and  lifetime.  The  technique  is  effective  over  a  very  wide 
range  of  m.t.  In  a  rather  poor-quality  film,  for  example,  we  have  measured  a  ilt 
value  as  low  as  5x10*'  crer/V. 


EXPERIMENTAL 

Device  quality  a-Si:H  films  were  prepared  by  rf  glow-discharge 

decomposition  of  silane.  First,  a  thin  n+  layer  was  deposited  on  a  conducting 
glass  substrate,  followed  by  a  1  p.m  thick  intrinsic  film.  The  n+  layer  provides 
an  ohmic  contact  between  the  conducting  glass  and  the  intrinsic  film.  The 
conducting  glass  was  cut  into  squares  of  about  O.S  cm'  in  area  and  a  contact 
attached  to  the  conducting  glas.s.  I'hc  wire  lead,  eontact  area  and  the  entire 
sample  were  masked  by  epoxy  cement,  except  for  a  small  area  (~  2  mm')  of 
the  film’s  surface  to  be  exposed  to  the  electrolyte.  The  sample  and  a  platinum  . 
electrode  were  immersed  in  an  indifferent  electrolyte  such  as  Ca(N03)2  or 
(NHJjSO..  The  sample  was  illuminated  through  the  (transparent)  electrolyte  by 
n  ~2  mw  He-Ne  laser.  The  measurement  technique  has  been  dc.scribcd 
elsewhere  [3],  and  will  be  reviewed  only  briefly  here.  A  short  (0.1-40  p»sec) 
voltage  pulse  applied  between  the  Pt  and  the  sample  is  used  to  charge  up  the 
semiconductor  space-charge  capacitance.  The  voltage  drop  across  these 

electrodes,  measured  just  after  the  termination  of  the  pulse,  represents  to  a  very 
good  approximation  the  change  6V-  in  barrier  height  across  the  film’s  space- 
charge  layer  induced  by  the  applied  pulse.  The  change  in  space-charge 
density  is  obtained  from  the  voltage  developed  across  a  series  capacitor,  again  at 
the  termination  of  the  pulse.  In  general,  the  induced  charge  may  include  also 
charge  leaked  across  the  S/E  interface  arising  from  imperfect  blocking.  Methods 
have  been  developed  to  determine  directly  the  leakage  charge  and  subtract  it,  if 
necessary,  from  the  measured  charge  so  as  to  yield  [3,4].  In  the  ca.se  of  a- 
Si:ll,  however,  the  interface  is  cs.scntiaily  perfectly  blocking  over  the  entire 
range  of  barrier  heights  studied.  Pulses  of  varying  amplitude  are  applied  singly, 
one  per  data  point  taken.  In  this  manner  electrochemical  reactions  that  may 
take  place  at  the  a-Si:H  interface  are  practically  eliminated  and,  what  is  more 
important,  damage  to  the  amorphous  film  is  minimized.  We  found  that 
applying  a  large  number  of  pulses  degrades  the  material,  drastically  reducing  its 
resistivity  and  introducing  large  trap  densities. 

In  order  to  derive  the  vs.  V.  curve  from  the  measured  vs.  8Vc 
data,  one  has  to  determine  ana  quiescent  barrier  height  and 

space  charge  density,  respectively.  is  determined  quite  accurately  (to  within 
±  20  meV)  from  measurements  in  the  depletion  range.  The  entire  \s.  Vj 
curve  can  then  be  constructed,  using  the  relations  =  V.q  +  8V-,  = 

Qsco  SQsc-  follows  the  surface  electron  density  N.  =  Qgjq,  where 

q  is  the  electronic  charge,  rather  than  Q^,  is  plotted  against 

The  dark  and  photo  resistance  of  each  sample  (Mtween  the  a- 
Si:H/elcctrolyte  interface  and  the  n+  contact)  was  derived  from  current-voltage 
characteristics  measured  at  the  onset  of  the  applied  pulse,  before  the  surface 
space-charge  capacitance  can  be  charged  through  the  sample’s  resistance  (see 
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below).  In  this  manner,  the  blocking  interface  is  effectively  shorted. 

RESULTS  AND  DISCUSSION 

Typical  results  of  the  surffice  electron  density  Ng  vs.  the  barrier  height  Vg, 
obtain^  for  an  illuminated  a«Si:H  sample  from  measurements  with  a  pulse  of 
3-(isec  (squares)  and  30-psec  (circles)  duration,  are  displayed  in  the  semilog 
plot  of  Fig.  1.  In  the  depletion  range  (Vg<0),  the  Ng  values  are  negative  and, 
because  of  the  logarithmic  scale  used,  the  plots  are  those  of  their  absolute 
magnitudes.  The  free  electron  concentration  under  illumination  n^^j  in  the 
amorphous  film  has  been  derived  from  the  measured  photoconductivity  on  the 
assumption  that  the  electron  mobility  p.  is  10  cmWsec  [6,7].  The  corresponding 
theoretical  dependence  of  the  surface  electron  density  N.  on  barrier  height  Vg, 
as  obtained  from  a  solution  of  Poisson’s  equation  [S]  in  the  absence  of  localized 
states,  is  shown  by  the  dashed  curve.  The  experimental  points  are  seen  to  lie 
well  above  this  curve  indicating,  as  expected  for  amorphous  films,  that  the 
space-charge  layer  in  both  the  depletion  and  accumulation  ranges  is  dominated 

localized  states.  ^ 

When  a  depletion  layer  7 
is  formed  (by  an  applied  I 
negative  pulse),  free  and 
trapped  electrons  are  ^ 
expelled  from  the  surface  ^ 
region,  leaving  behind  the  ^ 
positively-charged  localized  g 
states  that  make  up  the  2 
immobile  space-charge  in  the  g 
depletion  layer.  Consider  ^3 
first  the  case  of  a  non-  j 
illuminated  sample.  The  first 
process  that  occurs  following  u 
the  application  of  the  ^ 
negative  pulse  is  the  ^ 
expulsion  of  free  electrons 
from  the  surface  region.  As 
a  result,  there  is  a  net  Figure  1.  Surface  electron  density  N.  vr. 

thermal  emission  of  trapped  barrier  height  Vg  in  an  illuminated  sample  of 

electrons  from  localized  a-Si:H. 

states  into  the  conduction 

band,  and  the  emitted  electrons  are,  in  turn,  swept  away  from  the  space-charge 
region.  Only  shallow  states,  down  to  0.3-0.4  eV  below  the  conduction  band 
edge,  are  able  to  do  so  within  the  measurement  time  (the  pulse  duration)  [Sj. 
Under  illumination,  on  the  other  hand,  the  situation  is  quite  different.  Hole- 
electron  pairs  are  continuously  generated  by  the  light.  The  electrons  are 
expelled  from  the  surface  region,  while  the  holes  are  attracted  to  the  surface 
where  they  can  recombine  with  the  trapped  electrons.  In  this  manner  electrons 
in  occupi^  states  throughout  the  energy  gap  can  be  discharged  and  expelled 
from  the  space-charge  layer,  irrespective  of  Ae  depth  of  the  states.  The  only 
factors  that  determine  whi(^  and  how  many  states  can  release  their  trapped 
electrons  within  the  measurement  time  are  the  hole  capture  cross  sections  of 
the  states  and  the  hole  trapping  kinetics.  Referring  to  Fig.  1,  we  see  that  for  a 
measurement  time  of  3  iisec  (squares)  only  a  fraction  of  the  trapped  electrons 
are  expelled,  but  when  the  pulse  duration  is  extended  to  30  jisec  (circles). 


practically  all  occupied  states  in  the  depletion  layer  are  discharged  and  expelled 
(see  below).  The  solid  curves,  passing  in  each  case  through  the  points,  were 
calculated  from  Poisson’s  equation  for  the  case  in  which  localized  states  of 
densities  (as  marked)  are  present  [5].  The  agreement  between  theory  and 
experiment  is  seen  to  be  remarkably  good,  adding  considerable  confidence  to 
our  analysis.  The  total  density  of  occupied  states  in  the  entire  energy  g^,  as 
derived  from  the  data  for  the  30  (isec  measurement  time  is  about  lO'^cm*^, 
which  is  the  expected  value  for  the  films  studied. 

When  accumulation 
layers  are  formed, 
unoccupied  states  are  filled  f'' 
up.  Here  again,  .the  E 
experimental  points  are  seen 
to  be  initially  well  above  the  ° 
dashed  curve  calculated  for 
the  case  of  no  localized  ^ 
states  present,  indicating  that  <£ 
localized  states  dominate  the  ^ 
space  charge  layer  in  this 
range  as  well.  Surface  states  < 
are  probably  also  involved.  ^ 

At  stronger  accumulation  o 
conditions,  however,  the  ^ 
calculated  curve  crosses  the  iii 
experimental  points  and  x 
climbs  above  them.  This  is  ^ 
not  understood  at  present,  «  ^  ^ 

and  more  work  is  in  progress  Figure  2.  Expelled  trap  density  vs.  pulse 

to  account  for  such  a  duration  (measurement  time)  for  three  a- 

behavior.  Further  studies  Si:H  samples  under  illumination, 

are  needed  also  to 

distinguish  between  surface  ^ 
and  bulk  states.  7  .|  q 

It  should  be  noted  that  | 
the  quiescent  barrier  height 
is  nearly  zero.  This  is  to  be  z"  i  o 
expected  since  illumination  ^ 
tends  to  flatten  the  bands  at  iTi 
the  surface  (5).  5  1 0 

The  expelled  trap  density  ^ 
under  depletion  conditions,  § 
as  derived  from  ^  1 0 

measurements  such  as  those 
shown  in  Fig.  1,  is  plotted  in  d  _ 

Fig.  2  against  pulse  duration  u  '  ^ 

(measurement  time).  The  ^ 
lower  two  curves  were  k  ig 
obtained  for  high-grade  a-  ^ 

Si:H  films,  while  the  upper 
curve  was  obtained  after  a  Figure  3.  Surface  electron  density  vs. 

film  has  been  degraded  by  barrier  height  V.  in  an  unilluminated, 

the  application  of  many  degraded  sample  of  a-Si:H. 
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voltage  pulses.  In  all  cases,  the  expelled  trap  density  increases  with  pulse 
duration,  but  tends  to  saturate  at  a  pulse  duration  of  30-40  usee.  We 
interpret  the  saturation  level  as  representing  the  total  density  of  occupied  states 
in  the  energy  gap.  We  are  in  the  process  of  developing  a  model  involving  the 
hole  trapping  kinetics  in  order  to  account  for  the  shape  of  the  experimental 
curve.  Such  a  model  might  yield  an  estimate  of  the  average  cross  section  for 
hole  capture  the  occupied  states. 

The  results  depicted  in  Fig.  3  are  similar  to  those  in  Hg.  1,  except  that  jthey 
were  obtained  for  a  non-illuminated  sample.  The  low  dark  resistivity  necessary 
for  the  application  of  our  pulse  technique  has  been  achieved  in  this  ca.se  by 
degrading  the  film  (see  above).  Here  we  And  a  rather  high  density  of  occupied 
states,  as  expected  for  poor-quality  material.  Since  the  sample  is  in  the  dark, 
the  trapped  electrons  must  now  be  discharged  through  thermal  emission  into 
the  conduction  band.  As  such  these  states  must  be  not  deeper  than  0.3 -0.4  eV 
below  the  conduction-band  edge  [S]. 

The  dark  and 

photoconductivity  of  a 
degraded  a-Si:H  sample  vs. 
applied  pulse  voltage  is 
displayed  in  Fig.  4  on  a 
semilog  plot.  The 

measurements  were  taken  at 
the  onset  of  the  pulse, 
before  the  space-charge 
capacitance  at  the  a-Si:H 
interface  can  charge  up.  For 
positive  pulses  (electrolyte 
positive  with  resp^  to  the 
n+  contact),  the  conductivity 
in  both  the  dark  and  under  8 
illumination  is  seen  to 
increase  with  applied  voltage. 

This  is  due  to  electron 
injection  from  the  n+  layer. 

In  the  negative  polarity,  the 
dark  conductivity  is 
independent  of  voltage,  as  it 
should.  The 

photoconductivity,  on  the 

other  hand,  is  seen  to  decrease  with  pulse  amplitude,  approaching  the  dark 
conductivity  at  large  negative  voltages.  We  attribute  this  behavior  to  electron 
sweep-out  by  the  applied  Aelds,  as  shown  by  the  following  considerations.  The 
continuity  equation  [5]  for  the  case  of  a  blocking  contact  can  be  written  as: 
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Figure  4.  Dark  and  photo  conductivity  of  a 
degraded  a-Si:H  sample  vs.  applied  voltage 
pulse.  Solid  curve  in  the  negative  pulse 
range  represents  eq.  (3)  for  jat  =  5  x  10*® 
cm^. 


dN/dt  -  L  -  -  jiEN/d  ,  (1) 

where  N  is  the  total  density  of  photo  electrons  in  the  sample  per  cm^  L  is  the 
number  of  hole-electron  pairs  generated  by  the  light  per  cm^c,  E  is  the 
applied  Aeld  and  d  is  the  film’s  thickness.  Strictly  speaking,  the  last  term  in  cq. 
(1)  should  be  |iEn|,i(d),  where  nij|(d)  is  the  electron  concentration  at  the  n+ 
contact.  For  simplicity,  we  have  replaced  nM(d)  by  the  average  electron 
concentration  N/d,  which  is  a  fairly  good  approximation.  Space-charge  effects 


are  negligible  in  this  polarity,  even  if  all  photo  electrons  are  swept  away  from 
the  sample  (leaving  behind  the  photogenerated  holes).  Hence  the  field  is  very 
nearly  uniform  at  the  value  V/d.  Solution  of  eq.  (1),  with  the  boundary 
condition  that  N  assumes  the  zero*field  density  =  Lt  at  the  onset  of  the 
pulse  (t  =  0),  yields 

N  =  No/(l  +  M-xE/d)  +  Noexp[-(l/T  +  »tE/d)t]/(l  +  d/jtxE).  (2) 

Now  d/uE  is  the  electron  transit  time  through  the  film  and  is  typically  lU*^  sec 
for  an  applied  voltage  of  1  V.  The  measurement  is  taken  about  1  M.sec 
following  the  pulse  onset,  so  that  the  second  term  in  eq.  (2)  can  be  neglected 
to  a  very  good  approximation.  Re-writing  eq.  (2)  in  terms  of  the 
photoconductivity  a  =  qp.N/d  we  then  have 

a  =  o'q/(1  +  ftxV/d^)  (3) 

The  solid  curve  fitted  to  the  photoconductivity  data  at  negative  voltages  in  Fig. 
4  is  a  plot  of  eq.  (3)  for  jtx  =  5x10**  cmW.  The  fit  with  the  experimental 
points  is  seen  to  be  very  good.  The  rather  low  value  of  p,T  is  reasonable  for 
the  degraded  sample  u.sed  in  the  measurement.  The  method  can  be  applied 
effectively  even  if  \lt  is  an  order  of  magnitude  lower.  It  should  be  noted  that  ix 
in  this  context  is  the  trap-controlled  mobility  [7,8]. 
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ABSTRACT 

In  this  paper  we  show  that  the  degree  of  order  of  the  Si  network  in  a-Si:H  is 
increasing  with  two  length  scales  from  the  surface  into  the  bulk.  The  major  manifestation 
of  the  disorder  is  the  variation  in  the  Si>Si  bond>stretching  rather  than  the  variation  in  the 
width  of  the  dihedral  angle  distribution.  The  results  are  interpreted  in  terms  of  the 
decrease  of  the  hydrogen  concentration  from  the  free  surface  into  the  bulk. 

INTRODUCTION 

A  considerable  understanding  of  the  lattice  dynamics  of  covalent  amorphous 
semiconductors  in  general,  and  of  hydrogenated  amorphous  silicon  in  particular,  has  been 
gained'*’  between  the  mid  70’s  and  the  mid  80's.  Tbis  undemanding  is  bas^  on  the 
identification  of  the  phonon  spectra  in  these  materials  as  associated  with  a  coupling 
parameter-weighed  phonon  density  of  states  of  the  crystalline  one  phonon  modes^  . 
Following  this  development  the  Raman  spectra  of  a-Si:H  became  a  tool  in  evaluating  the 
"order"  of  the  amorphous  silicon  network.  Both  the  tool  and  the  "order"  are  simply 
defined  by  the  similarity  to  the  spectrum'  of  crystalline  silicon,  c-Si,  while  the  amorphicity 
of  the  structure  (i.e.,  the  corresponding  lack  of  a  long  ranee  order  or  the  presence  of  the 
corresponding  amorphous  radial  distribution  functioiv)  is  maintained.  The  most 
conspicuous  feature  of  these  Raman  spectra  is  the  phonon  density  band  which  originates 
from  the  transverse  optical  (TO)  Si-Si  stretching  mode  in  the  tetrahedral  structure  of  c-Si. 
The  closer  the  peak  position  of  the  corresponding  band  in  the  spectrum,  to  that  of 
the  c-Si  line  at  52U  cm  ',  and  the  narrower  the  peak  width,  A-ro,  the  more  "ordered"  the 
material.  The  microscopic  meanings  of  these  features  are  that  the  Si-Si  bond  length  is 
closer  to  that  of  c-Si,  and  that  the  width  of  the  bond’s  dihedral  angle  distribution  function, 
AO,  is  narrower,  in  particular,  models  were  presented’*’  and  predictions  were  given  for 
an  almost  linear  relationship  between  Ato  and  A6.  Once  this  model  became  well 
established,  attempts  were  made  to  relate  the  degree  of  order  to  other  physical 
properties'"’,  and  to  evaluate  the  clTect.s  of  various  film  de|M)siiion  piiranieters  on  the 
degree  of  disorder.  The  pur|XKic  of  the  pre.sent  study  is  to  follow  the  ty|x:  and  degree  of 
network  order  from  the  free  surface  of  a  deposited  a-Si:H  film  toward  its  bulk. 

Our  approach  in  the  present  work  is  to  study  both,  the  effect  of  the  laser  excitation 
frequency  and  the  effect  of  the  film  thickness  on  the  Raman  spectra.  The  need  for  such 
a  simultaneous  application  of  both  approaches  is  that  application  of  the  first  approach 
only,  does  not  enable  the  association  of  the  observations  with  other  physical  parameters 
(such  as  the  optical  band  gap'')  that  were  studied  in  conjunction  with  the  Raman  spectra. 
On  the  other  hand,  the  results  of  the  second  approach  do  not  necessarily  represent 
inhomogeneity  within  the  films.  For  example,  it  is  known  that  for  thick  films  there  is 
hardly  any  thickness  dependence  of  the  Raman  spectra”**’,  does  this  mean  that  these 
films  are  uniform?  Hence,  mutually  consistent  results  of  the  above  two  approaches  may 
reveal  whether  the  previously  reported  film  thickness  dependencies  of  the  Raman  spectra 
have  to  do  with  inhomogeneity  in  the  network  structure  along  the  axis  of  growth,  and 
indicate  the  reason  for  these  dependencies.  We  also  note  in  passing  that  while  film 
thickness  dependencies  of  the  Raman  spectra  have  been  studied  previously,  we  do  nut 
know  of  a  simultaneous  report  of  the  thickness  dependencies  of  the  TA-TO  peak 
intensities  ratio,  ItaAto>  the  ^ak  width  Ato>  and  the  peak  position  ta-ro  TO  mode. 
Consequently,  no  comprehensive  picture  regarding  the  film  thickness  dependence  of  the 
disorder  has  been  derived  thus  far. 
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Following  the  above  considerations  it  appears  that  only  a  ramprehensive  study  of 
the  various  features,  and  the  simultaneous  approach  adopted  in  the  present  study  may 
provide  maps  for  the  inhomogeneity  of  the  various  degrees  of  network  order  in  a-Sj:H 
films.  Following  the  derivation  of  the  maps  we  will  discuss  the  possible  mechanisms  that 
determine  the  spatial  distribution  of  the  disorder. 


EXPERIMENTAL  DETAILS 

The  samples  used  in  the  present  study  were  device  quality  a-Si:H  films  that  were 
deposited  by  rf  glow  discharge  decomposition  of  silane  under  standard  conditions'^'^  '* 
(e.g.,  decomposition  temperature  Tj  =  270°C  and  growth  rate  of  4A/s).  These  conditions 
were  kept  the  same  for  all  the  samples  used,  except  the  deposition  time  which  was  varied 
in  order  to  yield  films  of  different  thicknesses,  l^e  substrates  used  were  Coming  7059 
glass  slides.  The  most  relevant  characterization  of  the  samples  is  their  optical  band  gap. 
Consequently  we  show  in  Figure  1  the  Tauc-band  gap  of  the  films  (the  first  set  was  used 
for  our  Raman  spectra  study).  These  data  were  derived  from  the  measured  optical 
transmission  and  reflection  of  the  films. 

The  Raman  spectra  were  recorded  using  a  Raman  microprobe'*  (from  Instruments 
S.A.)  and  the  applied  laser  excitation  sources  were  S14.S  nm  radiation  from  a  6W  model 
Innova  906  Ar*  laser,  and  620nm  radiation  from  model  CR-S90  dye  laser,  both  from 
Coherent,  Inc.  The  measurements  were  carried  out  in  back  scattering  geometry,  using  a 
lOOX  microscope  objective. 

In  order  to  avoid  any 
heating  (or  degradation)  of 
the  sample,  the  laser  power 
at  the  sample  was  kept 
below  6  mW  and  the 
incident  beam  has  been 
slightly  defocussed  on  the 
samples.  The  samples’ 
condition  was  displayed  on  a 
video  monitor  interfaced  to 
the  microscope.  The 
detector  system  was  a  water 
cooled  1024  channel  diode 
array  from  Princeton 
Instruments.  The  raw  data 
was  smoothed  using  a  fast 
Fourier  filtering  program 
(Peakfit  V3.0)  developed  by 
AISN  Software  Inc. 


Figure  1.  The  Tauc-optical  gap  as  a  function  of  sample 
thickness  used  in  the  present  Raman  study  (first  set). 


EXPERIMENTAL  RESULTS 

In  Figure  2  we  show  typical  Raman  spectra  taken  by  the  application  of  the  above 
described  two  laser  excitations.  It  is  seen  that  for  this  10  /nm  thick  film  there  are 
significant  differences  between  the  two  corresponding  spectra  that  we  attribute  to  the  fact 
that  in  device  quality  a-Si:H  the  red  light  photons  have  an  absorption  depth  of  SOO  nm 
while  the  green  light  photons  have  an  absorption  depth  of  70  nm.  in  the  set  of  samples 
studied  here  the  band  gap  varies,  as  seen  in  Figure  1,  but  both  excitations  correspond  to 
above-optical  band  gap-absorption,  and  the  ratio  between  these  two  absorption  depths 
does  not  vary  significantly  within  this  set.  In  Figure  2  the  feature  around  ISO  cm ',  which 
is  known  to  be  associated  with  the  TA  bond  bending  mode,  splits,  in  contrast  to  many  and 
in  agreement  with  many  previous  reports.  I'he  most  interesting  observation  is  that  the 
relative  intensity  of  the  two  subpeaks  is  different  for  the  two  excitations.  Since  this  mode 
is  believed  to  involve  Si  triades*  and  Si-rings'^  this  finding,  which  has  not  been  reported 
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previously,  can  be  attributed  to 
different  relative  weights  of  these 
network  configurations'  in  the 
fllm  layer  adjacent  to  the  surface 
(the  green  excitation),  and  in  the 
bulk  (the  red  excitation).  We 
further  note  that  this  observation 
indicates  that  the  intensity  ratio 
of  the  peaks  ^  ^^11 

defined  in  order  to  be  used  as  a 
criterion  for  the  network  disorder. 

The  second  feature,  around  300 
cm'',  is  known  to  be  associated 
with  the  LA  mode  and  it  seems 
(by  the  above  argument)  to  be 
enhanced  in  the  surface  region. 

The  third  feature,  around  410  cm' 

',  which  is  known  to  be  associated 
with  the  LO  mode,  is  responsible 
for  the  asymmetric  broadening  of 
the  most  conspicuous  feature  in 
the  spectrum,  the  peak  of  the  TO  stretching  mode,  which  is  located  around  470  cm '.  The 
strength  of  the  LO  mode  in  the  film’s  layer  adjacent  to  the  surface  (a  larger  amplitude 
for  the  green  excitation)  may  be  interpreted  as  due  to  a  higher  concentration  of  Si  rings 
structures'*'^  on  the  surface  and  its  vicinity.  Such  higher  concentration  of  rings  can  be  the 
result  of  the  bond  bridging  that  takes  place  on  the  surface. 

Examining  the  TO  peak  position  in  both  spectra  shows  a  consistent  behavior  of 
both  (i>TO  and  The  decrease  of  Wfo  the  increase  of  A-ix)  with  respect  to  those  of 
c-Si  are  larger  for  the  green  excitation.  Considering  the  above-mentioned  criteria  for  the 
short  and  intermediate  range  order  of  the  Si  network^'*  we  may  conclude  that  the  top  70 
nm  layer  (the  depth  of  the  green  light  absorption)  is  less  ordered  than  the  bulk  material. 
Note  that  this  is  also  true  for  a  10  um  sample  for  which  a  film  thickness  dependence 
study  can  be  wrongly  interpreted  as  indicating  a  homogeneous  network  order  throughout 
the  film. 

The  above  conclusions  are  following  the  premise  that  the  different  behavior  of  the 
red  and  green  excitations  in  Figure  2  cannot  be  attributed  to  a  resonance  Raman 
scattering.  The  basic  argument  for  this  premise  is  that  the  two  excitations  used  here,  are 
well  removed  from  the  resonance  that  occurs\  in  a-Si:H,  for  photon  energies  larger  than 
3  eV.  If  there  is  any  contribution  of  the  "tail"  of  the  resonance  that  affects  our  results  it 
is  clearly  being  offset  by  the  penetration  effect  argued  above.  The  evidence  for  that  is  the 
fact  that  is  our  observation  that  Ujo  decreases,  the  Ato  increases,  and  the  Ita/I-io  increases 
with  increasing  laser  frequency,  while  in  the  pure  resonant  case  the  reverse  is  true. 
Furthermore,  while  the  switching  in  the  TA  subpeaks  amplitudes  can  be  easily  attributed 
to  the  difference  between  a  bulk  and  a  surface  response  (see  above)it  cannot  be 
attributed  to  a  resonant  effect.  As  we  show  below,  the  interpretation  of  the  above  data 
in  terms  of  the  film  inhomogeneity  is  also  consistent  with  the  following  findings  regarding 
the  thickness  dependence  of  the.  features,  mentioned  above. 

Let  us  turn  then  to  the  thickness  dependence  of  wjo.  As  it  is  clearly  shown  in 
Figure  3,  (i>to  increases  within  the  first  couple  of  microns  of  film  thickness,  and  then 
saturates  for  both  excitations.  This  increase  is  relatively  faster  for  the  green  excitation  in 
comparison  with  that  of  the  red  excitation.  We  interpret  these  results  as  indicating  that 
an  "equilibrium"  degree  of  network  disorder  for  the  corresponding  property  is  established, 
both  at  the  "free  surface"  vicinity  and  in  the  bulk,  only  for  depths  larger  than  a  couple  of 
microns.  If  we  assume  that  a  thin  sample  is  similar  in  its  properties  to  the  top  layer  of  a 
thick  sample  this  finding  is  consistent  with  the  findings  of  Figure  2.  The  range  of  variation 
of  the  network  order  is  also  consistent  with  the  homogeneity  scale  associated  with  other 
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Figure  2.  Typical  Raman  spectra  of  a-Si:H  films 
under  the  excitation  of  two  different  laser 
frequencies. 
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physical  properties  of  these  films, 
such  as  the  optical  band  gap 
(Figure  1)  and  the  phototransport 
properties'*.  For  the  derivation  of 
the  Ato  value,  so  that  the  result  is 
not  affected  by  the  contribution 
of  the  LO  mode,  we  have 
considered  in  this  study  the  value 
of  twice  the  right  hand  side  width 
of  the  TO  peak  as  done  in 
previous  studies*-'*’'*.  We  find 
that  Ato  is  quite  insensitive  to  the 
sample  thickness  but,  as  cleaYly 
seen  in  Figure  4,  it  is  larger  for 
the  green  excitation  than  for  the 
red  excitation.  We  note  in  passing 
that  the  independence  of  A^  on 
d  is  reminiscent  of  the 
independence  found'*  on  7$  or  on 
the  band  gap*,  in  a-Si:H  films 
made  by  glow  discharge 
decomposition  of  silane  when  the 
hydrogen  content,  C^,  is  less  than 
20  at  %  (see  below). 

The  main  conclusion  from 
the  results  shown  here  is  then, 
that  in  device  quality  materials, 
A6  is  a  constant.  Considering  the 
theoretical  predictions  and  the 
experimental  values'*'*  wc  cati 
conclude  that  the  hulk  A6  is  as 
narrow  as  it  can  be  (i.e.,  9"), 
while  at  the  vicinity  of  the  surface 
A6  has  not  relaxed  to  the  bulk 
optimized  value.  On  the  other 
hand  the  fact  that  the  surface  A6 
is  also  independent  of  the  history 
of  the  sample  growth,  indicates 
that  there  is  a  "surface  optimized" 
A  6  which  is  achieved  in  device 
quality  a-Si:H. 

Turning  to  the  1|w/Ito 
ratio,  which  we  denote  here  by 
Ata/Ato  to  indicate  that  the 
corresponding  data  is  obtained  by 
taking  the  ratio  of  the  areas 
under  the  corresponding  peaks, 
we  get  the  results  shown  in 
Figure  S.  Other  evaluations  of 
Ita/I-io  ratios  and  their  meaning 
will  be  discussed  elsewhere.  The 
overall  decrease  of  the  A^/A^a 
ratio  with  thickness  and  the  larger 
ratios  for  the  green  excitations 
are  consistent  with  previous 
results,  obtained  on  glass 
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Figure  3.  The  thickness  dependence  of  the  TO  peak 
position  for  the  two  applied  laser  excitations. 
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Figure  4.  The  thickness  dependence  of  the  TO  peak 
width  under  the  two  applied  laser  excitations. 
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Figure  5.  The  ratio  of  the  areas  under  the  TA  and 
the  TO  peaks,  under  the  two  applied  laser 
excitations. 
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subst^ates'^  and  with  the  general  decrease  of  network  disorder  from  the  surface  into  the 
bulk. 

DISCUSSION  AND  CONCLUSIONS 

The  data  presented  seems  to  indicate  quite  clearly  that  device  quality  a-Si:H  is  an 
inhomogeneous  system  in  the  order  of  the  Si  network  along  the  film’s  growth  axis.  The 
inhomogeneity  has  two  depth  scales,  one  that  is  less  than  70  nm  and  the  other  that  has 
a  tail  of  a  couple  of  microns  from  the  free  surface  into  the  bulk.  These  two  length'scales 
should  be  considered  when  analysis  of  a-Si:H  devices  is  carried  out.  The  most  important 
conclusion  of  the  present  work  is  that  the  Raman  spectra  reported  in  the  literature,  which 
are  normally  taken  using  514  nm  wavelength  excitation,  reflect  the  degree  of  network 
order  within  the  first  length  scale.  Deeper  in  the  bulk  the  amorphous  network  has  a 
considerably  larger  degree  of  order. 

Turning  to  the  various  features  of  disorder  and  their  meaning  let  us  start  with  the 
Ato  thickness  independence  shown  in  Figure  4.  The  fact  that  Ato  has  been  found  to  be 
independent  of  for  other  films  prepared  by  glow  discharge such  that  the  variation 
of  Eg  is  the  same*^  as  in  our  films  (see  Figure  1 ),  yields  the  conclusion  that  the  dihedral 
angle  distribution  disorder  is  optimized  under  the  preparation  conditions  of  device  quality 
a-Si:H.  The  fact  that  this  is  in  sharp  contrast  with  the  dependence  found  in  sputtered 
films'*'^'"’  and  the  fact  that  A  iq  Is  of  the  order  expected  for  the  minimum  piossible  A6  is 
indicating  that  no  more  improvement  can  be  made  in  narrowing  A6  by  variation  of 
deposition  conditions.  On  the  other  hand  the  constant  difference  between  A^  values  for 
the  two  excitation  frequencies  u.sed  can  be  interpreted  as  due  to  the  sharp  difference 
between  the  close  vicinity  of  the  free  surface  (70  nm)  and  the  bulk  of  a-Si:H.  While  not 
mentioned  previously  in  the  literature,  this  finding  is  not  to  surprising  since  the  bond 
reconstruction  near  the  surface,  due  to  the  termination  of  the  tetrahedral  network,  is 
quite  different  from  that  of  the  bulk.  The  present  results  show  then  that  the  A6  width  at 
the  surface  vicinity  is  also  optimized,  but,  as  to  be  expected  the  surface  optimized  value 
is  larger  than  the  bulk  optimized  value. 

The  results  for  the  thickness  dependence,  which  are  similar  to  those 

obtained  by  other  researchers  (for  glass  substrates^'*),  reveal  the  variation  of  this  ratio  on 
a  .scale  of  a  couple  of  microns,  ju.st  as  the  Ej,(d)  variation  shown  in  Figure  1.  We  find 
sigain  a  clesir  difference  between  the  close  vicinity  of  the  surface  iind  the  bulk.  Since  the 
TA  mode  is  associated  with  the  presence  of  Si  triad.s'  this  conclusion  is  consistent  with 
the  variation  of  the  photoelectronic  properties  which  depend  on  the  concentration  of 
dangling  bonds.  One  can  indeed  show  that  there  are  two  length  scales  involved  in  the 
latter  concentration.  The  first  is  associated  with  the  surface  recombination  velocity'*  and 
the  other  with  a  broader  decrease  of  this  concentration  toward  the  interior  of  the  bulk^. 
The  agreement  of  the  above  observations  with  those  found  in  the  study  of  the  Raman 
spectra  as  a  function  of  T^,  for  materials  of  the  same  optical  band  gap"' is  an  indication 
that  it  is  the  Eg  dependence  that  should  be  used  for  comparisons  of  features  in  the 
Raman  spectra.  In  particular,  since  the  dihedral  angle  order  is  shown  here  to  be 
optimized,  it  appears*  that  the  prime  reason  for  the  variation  of  Eg(d)  is  due  to  the 
viiriation  of  C„.  It  is  also  well  known  that  C„  decreases  from  the  surface  into  the  bulk''. 
Correspondingly,  there  is  the  decrease  of  the  A|a/A|,,  ratio  with  d,  and,  as  discus.sed 
above,  the  corresponding  decrease  of  the  dangling  bond  concentration.  This  conclusion 
is  also  consistent  with  the  dependence  on  d  as  given  in  Figure  3,  since  again,  the 
results  are  consistent  with  the  ii>to  dependence  on  T^.  The  explanation  for  the  latter 
results  was  given  by  the  model  of  the  quasi  interstitial  hydrogen.  This  model,  as  proposed 
by  Hishikawa  et  al’*,  suggests  that  the  main  role  of  the  hydrogen  atom  is  to  act  as  an 
interstitial  impurity  in  the  silicon  network,  yielding  repulsive  forces  between  itself  and  the 
neighboring  silicon  atoms.  Therefore  the  much  larger  hydrogen  concentration  at  the 
surface  and  its  gradual  decrease  into  the  bulk  account  well  for  the  observation.  A  model 
based  on  network  relaxation,  on  the  other  hand,  cannot  explain  simultaneously  the 
independence  of  A^  ond  the  dependence  of  on  d. 
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a  reprfait  froM  Ai^lkd  Opties 


T-matrix  approach  for  calculating  local  fields 
around  clusters  of  rotated  spheroids 


William  Vargas,  Luis  Cruz,  Luis  F.  Fonseca,  and  Manuel  G6mez 


A  T-matriz  fonnalism  ia  uaed  to  calculate  local  electric  fielda  around  cluatera  of  prolate  spheroids  in  the 
long-wavelength  regiine.  The  calculationa  are  performed  aa  a  function  of  interparticle  distance  as  well  as 
an^e  of  orientation.  The  ohaerved  red  shifts  in  the  reaonant  wavelengths  of  the  characteristic  peaks  are 
shown  to  obey  an  exponential  relationahip  as  a  function  of  interparticle  aeparation  and  a  sinusoidal 
relationahip  as  a  function  of  angle  of  rotation  of  the  spheroid.  The  behavior  of  the  cluster  is  rfiin— «< 
and  the  two  effecta  of  aeparation  and  rotation  are  compared. 


1.  Introduction 

Optical  properties  of  systems  of  particles  whose  size  is 
smaller  than  the  wavelength  of  the  incident  electro¬ 
magnetic  radiation  have  been  studied  for  a  long 
time,*-*  but  the  effect  of  the  clustering  of  such  parti¬ 
cles  is  not  yet  well  understood.  Recent  interest  has 
developed  in  this  subject  because  clustering  is  essen¬ 
tial  to  the  understanding  of  many  physical  processes 
of  interest.  Moreover,  the  current  calculations  of 
local  scattered  electric  fields  are  directly  applicable  to 
fields  such  as  surface-enhanced  Raznan  scatte^n'ig 
and  light  scattering  and  absorption  &om  metaiiic 
colloids  in  alkali  halide  crystals.*-* 

Clusters  of  spheres  have  been  analyzed,  primarily 
in  the  radiation  zone*-^*  or  with  the  use  of  an 
electrostatic  approach.**-**  The  purpose  of  this  study 
is  to  calculate  the  local  electro^namic  field  around 
clusters  of  two  prolate  spheroids  in  the  long-wave¬ 
length  regime  as  a  function  of  the  separation  ^tween 
them  and  of  their  relative  orientations  by  using  a 
T-matriz  approach,  which  describes  electromagnetic 
scattering  for  a  general  wave  that  is  incident  on 
objects  of  arbitrary  shape.  Althou^  the  formalism 
can  be  used  to  calculate  important  parameters  in  the 
radiative  regime,  such  as  scattering  cross  sections, 
the  local  fiel^  are  the  migor  concern  here. 

The  T-matriz  method  permits  an  electrod}mamic 
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approach  that  considers  vectorial  multipolar  fields,  in 
contrast  to  other  recent  electrostatic  approaches  that 
consider  scalar  multipolar  potentials.**  The  ezten- 
sion  to  two  scatterers  is  achieved  by  using  an  effective 
T  matrix  for  more  than  one  scatterer,  which  retains 
all  the  advantages  of  the  T  matrix  for  the  single 
scatterer  and  alM  permits  calculations  for  scatterers 
made  up  of  clusters  of  metallic  particles  of  arbitrary 
shape.  The  rotation  of  the  individual  T  matrices 
corresponding  to  any  given  particle  is  achieved  by 
using  a  transformation  of  the  spherical  harmonics 
under  finite  rotations.  This  stu^  extends  previous 
research  on  the  calculations  of  the  scattered  electric 
field  in  the  vicinity  of  clusters  of  metallic  scatterers.** 
In  this  vicinity  two-particle  clusters  of  various  config¬ 
urations  gave  enhancements  of  up  to  10*  and  demon¬ 
strated  the  existence  of  a  criti^  distance  greater 
than  the  touching  distance  between  the  scatterers  at 
which  maximum  enhancement  was  obtained. 

We  demonstrate  that  because  of  clustering  interac¬ 
tions  the  calculated  spectrum  shows  two  peaks  for  the 
case  of  clusters  of  two  small  prolate  spheroids  aligned 
in  the  direction  of  the  incident  electric  field,  in 
contrast  to  the  spectrum  of  the  isolated  small  prolate 
spheroid  that  has  only  one  peak.  The  prindj^  peak 
decreases  in  intensity  and  shifts  exponentially  as  a 
function  of  interparticie  distance  to  the  value  of  the 
resonance  for  the  isolated  spheroid,  whereas  the 
secondary  peak  has  a  vanishing  intensity  that  de¬ 
creases  exponentially  as  a  function  of  interparticie 
distance.  A  decoupling  interparticie  distance  d  of 
46,  where  6  is  the  seminugor  axis  of  the  prolate 
spheroids  that  make  up  the  cluster,  is  observed.  At 
that  distance  the  behavior  of  the  cluster  near  any  one 
of  the  two  constituent  particles  equaL  that  of  one 
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ilolated  spheroid.  It  is  shown  thst  there  is  also  a 
shift  in  wavelength  tp  the  red  as  a  result  of  the 
rotation  of  the  q;>heroids  and  that  this  shift  varies 
sinusoidaUy  as  a  function  of  orientation. 

For  all  the  calculations  heron,  silver  was  taken  as  a 
model  metal  by  using  the  complex  frequency-depen¬ 
dent  dielectric  function  report^  in  the  literature  for 
the  bulk  material, althou^  the  method  is  applica¬ 
ble  to  any  material  for  which  the  complex  dielectric 
constant  is  known.  For  the  range  of  edzes  used  in 
these  calculations,  an  electronic  mean  free  path  re¬ 
duced  by  the  boundary  of  the  scatterer  or  quantum 
size  effects  could  introduce  some  modlfleations  to  the 
data  for  silver  used  here. 

Section  2  presents  the  T-matrix  formalism  devel¬ 
oped  for  isolated  particles  and  its  modification  for 
metallic  clusters.  The  calculations  are  done  as  expan¬ 
sions  in  vector  spherical  harmonics  whose  num¬ 
ber  n  corresponds  to  terms  in  the  multipolar  expan¬ 
sion  of  the  electromagnetic  field.  Because  our  main 
objective  is  to  study  clusters,  spherical  harmonics  are 
us^  as  a  convenient  basis.  A  rotation  of  the  T 
matrix  corresponding  to  a  sin^e  scatter  is  performed, 
permitting  consideration  of  geometries  where  the 
spheroids  in  the  cluster  have  arbitrary  relative  orien¬ 
tations. 

Section  3  presents  the  calculation  of  the  local 
electric  fields  for  these  clusters.  The  clusters  are 
analyzed  as  a  function  of  separation  between  sphe¬ 
roids  and  as  a  function  of  ang^e  of  orientation  with 
respect  to  the  incident  wave.  The  rate  of  change  of 
the  wavelength  shift  as  a  function  of  interparticle 
distance  and  angle  of  orientation  is  then  discussed. 
In  Section  4  conclusions  are  presented. 


2.  Formalism 
A  Single  Scatterers 

The  T-matrix  formalism  developed  by  Waterman'^ 
takes  into  account  multipolar  contributions,  which 
are  essential  for  any  valid  calculation  of  local  fields  of 
sin^e  nonspherical  particles  as  well  as  for  all  clusters 
even  when  they  are  in  the  long-wavelength  regime. 
The  formalism  also  takes  into  account  phase-retarda¬ 
tion  effects  caused  by  the  size  of  the  scatterers,  which 
are  important  for  particles  and  clusters  whose  sizes 
are  comparable  with  the  wavelength  of  the  incident 
field. 

In  this  method  the  scattered,  internal,  and  incident 
fields  are  expanded  in  terms  of  the  corresponding 
elementary  fields  that  are  a  basis  set  of  solutions  for 
the  vector  Helmholtz  equation, 

V  X  V  X  0.  (1) 

The  incident  So  and  internal  e,  fields  are  expanded  in 
terms  of  the  basis  that  is  regular  at  the  origin  Re 
and  the  scattered  e,  field  in  terms  of  the  nonregular  ^ 


one. 

So  - 

fl 

*1  -  S  I>J.  Re |rl  <  r^, 

H 

e.  “  2  ir|  >  raa,  (2) 

II 

where  are  known  coefficients,  D„  and  F„  are 
unknown,  is  the  radius  of  the  maximum  sphere 
inscribed  in  the  scatterer,  and  Rmw  i*  the  minimum 
sphere  inscribing  the  scatterer.  . 

The  elementary  wave  functions  are  expressed  as 

x)’l*iy^(/‘)h„(Ar)],  (3) 

where  t  *  1,  2,  a  «  even(e)  or  odd(o),  n  =  1,  2, . . . , 
m  *  0, 1, . . . ,  n, 

(2n  +  l)(n  -  m)! 

4n(n  +  l)(n  +  m)!  ’ 

ytmn(^)  *  coa(m<|>)P,"(cos  6), 
yamnl^)  -  sin(m4>)P„"*(cos  e). 

The  index  t  «  1,  2  describes  the  type  of  excitation, 
magnetic  or  electric;  is  the  Neumann  symbol 
defined  as  co  =  1  and  2  otherwise;  n  is  the  order 
of  the  multipole;  and  a  gives  the  parity  of  the 
elementary  functions.  The  regular  form  of  the  basis 
functions  are  obtained  by  substituting  the  Hankel 
functions  with  the  Bessel  functions. 

The  surface  currents  on  the  scatterers  are  used  to 
express  the  expansion  coefficients  of  the  internal  field 
with  those  of  the  scattered  and  incident  fields,  respec¬ 
tively,  by  the  following  relationships: 

F=-»Re(g')G,  (4) 

A  =  iS'G,  (5) 

where  G  represents  the  vector  of  the  expcmsion 
coefficients  of  the  internal  field  and  Q'  represents  the 
transpose  of  the  Q  matrix.  For  a  particle  with  a 
complex  dielectric  Iimction,  Q  is  given  by 

J*{[V  ^  Re'FJ^r)!  X  qr„.{fcor) 

-hRe>r,(*r)x{Vxir,.{fcor)l),  (6) 

where  and  **  *  s  is  the 

surface  of  the  scatterer,  and  'l'n(ikr)  is  substituted  by 
Re  '^„{kr)  wherever  Re  Q  appears.  In  our  case  e«Kt  * 
1  (vacuum)  and  eiat  is  the  corresponding  value  for 
silver. 

Equation  (5)  permits  the  calculation  of  the  internal 
field,  at  least  when  |r|  <  r„i„.  Eliminating  G  from 
Ik)8.  (4)  and  (5),  we  obtain  a  relation  between  the 
coefficients  of  the  scattered  and  incident  fields: 

F  -  TA,  (7) 
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w^iere  T  is  the  T  matrix  of  the  sin^  acatterer  defined  as 
T--9-‘Re§.  (8) 

With  Eq.  (2)  the  local  fields  can  be  calculated  for 
|r|  <  rguBbyusix^theQmatrizandfor  |r|  > 
using  the  T  matrix,  but  ^e  region  rom  >  |r|  >  Taunis 
forbidden  becarise  the  convergence  of  the  expansions 
of  Eq.  (2)  is  not  ensured.  Some  research  has  been 
done  to  improve  on  this  shortcoming;  Biingi  and 
Seliga*^  proposed  a  mathematics  procedure  that 
used  the  T-matrix  framework  to  calculate  fields  in¬ 
side  this  forbidden  zone.  Barber  et  al.^^  used  a 
method  based  on  the  T-matrix  approach  to  evSuate 
local  fields  not  only  at  the  tips  but  at  any  place  on  the 
surface  of  the  spheroids. 

b.  Rotation  of  the  T  Matrix 

For  the  rotation  of  the  T  matrix,  two  coinciding 
coordinate  s3rstems  are  considered.  With  respect  to 
the  origin  of  both  systems,  the  Euler  angles  a,  and 
y  are  to  rotate  one  of  the  systems  relative  to  the 
other.  The  rotation  of  the  basis  functions  is  accom¬ 
plished  by  transforming  the  spherical  harmonics 
under  finite  rotations.  These  transformation  proper¬ 
ties  are  given  in  the  quantum  mechanics  literature 
(e.g..  Ref.  20). 

Explicitly,  if  and  are  the  vector  spheri¬ 
cal  functions  in  the  rotated  and  nonrotated  coordi¬ 
nate  ^sterns,  respectively,  then  one  can  write 

Pi  ‘Y)'l^'mm'«i 
m'v 

where  the  index  n  is  invariant  imder  rotations  and 
the  rotation  matrix  D  is  a  function  of  the  Euler  angles 
that  define  the  relate  orientation  of  the  two  coordi¬ 
nate  system.^ 

The  incident  and  scattered  electric  fields  in  Eq.  (2) 
are  then  expanded  with  respect  to  the  two  coordinate 
qrstems.  A  relationship  between  the  expansion  coef¬ 
ficients  for  Co  and  e,  in  the  two  systems  can  be 
established  by  using  ^e  orthogonality  of  the  vector 
spherical  functions.  Using  the  definition  of  the  T 
matrix,  we  finally  obtain 

T'-D-‘(a,P,7)TD(a,p,7).  W 

where  T'  is  the  T  matrix  for  the  rotated  coordinate 
system. 

C.  Ousters 

The  T-matrix  formalism  has  been  extended  to  sys¬ 
tems  with  more  than  one  scatterer  by  Peterson  and 
Str6m*  by  using  the  translation  theorems  for  the 
vector  spherical  functions.**  The  translation  proper¬ 
ties  oftb.  and  Re  ibn  are  summarized  by  Raf-  22: 

Re  ♦„(r  +  a)  -  2  B,,,,-,  (a)Re 

t'h* 

+  a)  -  2  <r«.»»(a)Re  'F,.,'(r),  |a|  >  lr|, 

T'n* 

+  a)  -  2  |a|  <  |r|. 

t'h' 


where  and  are  the  elements  of  the 

translation  matrices  as  defined  in  Ref.  9. 

Peterson  and  Strtm  obtained  a  T  matrix  for  the 
cluster  of  two  partides  in  terms  of  the  T  matrices  of 
each  sin^e  scatterer, 

T(l,  2)  =  R(a,)(T(l)[l  -  o(-a,  +  a*) 

X  T(2)a(-aj  +  ai)T(l)]-‘ 

X  [1  +  o(-ai  -e  a2)T(2)^ai  -  aj)]) 

X  R(-ai)  +  R(a2)(T(2)[l  -  of-aj  +  aj) 

X  T(l)^-a,  +  a2)T(2)]-‘[l  +  CT(-a2  +  aj) 

X  T(l)R(a2  -  ai)llR(-a2),  (10) 

where  ai  and  a2  are  the  distances  from  the  origin  to 
the  center  of  scatterers  1  and  2,  respectively. 

The  scattered  field  can  be  expressed  in  terms  of  the 
incident  field  by  using  the  T-matrix  with  the  follow¬ 
ing  relationship: 

B,  -  TEq. 

Most  of  the  results  presented  herein  will  be  ex¬ 
pressed  in  terms  of  the  total  electric  field  given  by 

E,  =  E.  +  Efl.  (11) 

It  is  important  to  remember  that  the  local  fields  can 
only  be  calculated  by  starting  from  a  minimum 
circumscribing  sphere  around  the  scatterers  to  en¬ 
sure  conver{^nce  of  the  spherical  wave  expansions. 

With  these  two  transformations,  rotation  and  trans¬ 
lation,  it  is  then  possible  to  consider  clusters  with 
variables  separation  and  arbitrary  orientation  of  the 
constituent  metallic  spheroids. 

3.  Discussion 

The  quantities  discussed  here  are  the  normalized 
local  field  intensities  |E(/eo|*  in  the  nei^borhood  of 
the  clusters  of  two  prolate  spheroids.  These  have 
been  calculated  by  inserting  the  individual  T  matrices 
of  each  spheroid,  Eq.  (8),  into  Eq.  (10),  thxis  yielding 
the  effective  T  matrix  for  the  cluster  as  a  function  of 
interparticle  distance.  In  the  c^  of  a  different 
orientation  between  the  spheroids,  their  rotation  has 
been  taken  into  account  by  applying  Eq.  (9)  and 
substituting  the  rotated  T  matrix  into  Eq.  (10). 
Finally,  the  total  electric  field  vector  has  been  c^cu- 
lated  firom  Eq.  (11),  from  which  the  enhancement  of 
the  intensity  of  the  field  was  obtained  by  taking  the 
square  of  the  magnitude  of  that  vector. 

To  illustrate  the  general  results  from  this  type  of 
calculation,  we  use  two  identical  prolate  spheroids 
with  an  aspect  ratio  a/b  >  0.9  and  o  *  5  nm,  where  o 
and  b  are  the  semiminor  and  seminugor  axes  of  the 
prolate  spheroid,  respectively.  The  formalism  per¬ 
mits  the  treatment  of  spheroids  of  any  eccentricity, 
but  because  of  computational  time  and  memoiy  limi¬ 
tations  only  the  described  size  is  considered. 
Althou^  the  clusters  of  spheroids  used  here  are  in 
the  long-wavelength  regime,  the  calculations  are  elec- 
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trodlynailuc  in  nature  because  the  formalism  is  devel¬ 
oped  from  the  vector  Helmholtz  Eq.  (1);  results  can 
and  have  been  obtained  for  other  sizes  and  aspect 
ratios.  “ 

The  calculations  were  done  1^  taking  the  direction 
of  propagation  of  the  incident  wave  front  to  be  along 
the  y  axis.  The  polarization  of  the  electric  field  is 
parallel  to  the  direction  of  the  interparticle  distance, 
which  is  the  x  axis.  By  testing  other  directions  of 
incidence  and  polarizations,  we  have  determined  that 
the  latter  orientation  gives  the  hipest  scattered- 
intensity  enhancement  (of  the  order  of  1(H)  for  the 
configuration  of  the  clusters  considered;  we  therefore 
choose  these  as  our  working  parameters.  The  calcu¬ 
lations  of  the  electric  field  are  performed  as  a  function 
of  interparticle  separation  and  relative  orientation  of 
the  spheroids. 

Two  resonant  peaks  are  exhibited  by  the  spectra  of 
all  the  geometrical  configurations  consider^.  It  is 
worth  noting  that  this  behavior  is  due  to  the  interpar¬ 
ticle  multipolar  interactions,  because  only  one  peak  is 
observed  when  the  spheroids  are  sufficiently  far  apart 
(when  d,  the  interparticle  distance,  is  of  the  order  of 
46). 

The  two  resonant  peaks  experience  a  characteristic 
red  shift  as  a  function  of  decreasing  distance  between 
spheroids  for  all  the  clusters  considered.  This  shift 
reflects  the  interaction  among  the  spheroids  in  the 
cluster,  which  becomes  stronger  as  the  distance  be¬ 
tween  the  centers  of  the  particles  decreases,  thereby 
resulting  in  maximum  shift  for  touching  spheroids. 
Figure  1  shows  the  calculated  spectrum  for  the 
cluster  of  two  spheroids  whose  semimqjor  axes  are 
aligned  with  the  incident  electric  field.  Because  this 
is  the  more  elongated  cluster,  it  exhibits  the  highest 
enhancement  shown  in  this  paper  (of  the  order  of 
1(H).  The  closer  the  spheroids  the  hi|d>er  the  interac- 


Fig.  1.  Intenaity  of  the  total  electric  field  at  the  point  of  obeerva- 
tion  P  aa  a  function  of  wavelength  for  a  cluater  of  aligned  apheroida 
with  a/b  B  0.9  and  o  •  5  nm  for  three  interparticle  diatancea 
(d  •  11.112  nm  »  26.  touching  apheroida.  11.612  nm,  and  12.112 
nm).  The  direction  of  the  incident  radiation  ia  along  they  axia  and 
the  polarization  ia  parallel  to  the  x  azia. 


tion  between  them  and  the  more  multipolar  terms  are 
required  in  the  calcula^on.  As  the  inteiparticle 
dis^ce  decreases,  the  low-energy  peak  shifts  more 
rapidly  toward  ^e  red  than  the  hi^-energy  peak. 
T^  tehavior  is'related  to  higher  multipolar  terms, 
other  than  the  dipolar,  which  become  more  important 
as  the  spheroidal  constituents  of  the  cluster  approach 
each  other.  As  the  interparticle  distance  riiminiBhpy, 
the  low-ener^  peak  decreases  rapidly  in  energy  while 
its  intensity  increases.  The  hi^-energy  peak  disap¬ 
pears  when  d  *  36,  the  decoupling  distance,  whereas 
the  low-energy  peak  shifts  to  the  resonant  wave¬ 
length  corresponding  to  the  irolated  particle. 

The  effect  of  rotation  has  been  anal}rzed  by  rotating 
only  one  spheroid  and  by  examining  the  local  field  at 
opposite  points  on  the  exclusion  sphere  on  the  ajris 
that  joins  the  centers  of  the  particles,  because  for  this 
case  the  cluster  is  asymmetric.  The  analysis  is  also 
performed  by  rotating  both  spheroids,  thus  maintain¬ 
ing  the  mirror  symmetry  of  the  cluster.  In  both 
cases  the  interparticle  ffistance  is  held  constant. 
Figures  2(a)  and  2(b)  show  the  behavior  of  the  local 
field  observed  at  the  two  previously  indicated  points 
on  the  X  axis  when  only  one  spheroid  is  rotated. 
Figure  2(a)  describes  the  field  near  the  nonrotated 
spheroid  and  Fig.  2(b)  describes  that  near  the  rotated 
spheroid.  The  difference  between  the  two  observa¬ 
tion  points  is  in  the  intensity  of  the  local  field, 
whereas  the  position  of  the  resonance  peaks  remains 
constant.  This  Constance  is  to  be  expected  because 
the  resonant  peaks  are  phenomena  associated  with 
surface  plasmons  of  the  cluster  that  depend  on  the 
configuration  of  the  whole  cluster  and  not  on  the 
point  of  observation.  Figure  3  shows  the  spectrum 
of  the  simultaneously  rotated  spheroids.  The  config¬ 
uration  6  =  90°  coincides  with  the  cluster  of  Fig.  1, 
which  exhibits  the  highest  enhancement  in  this  paper. 
The  rotation  of  the  spheroids  does  not  affert  the 
number  of  p>eaks  observed  in  the  spectrum,  although 
their  energies  are  shifted  toward  the  red  as  a  function 
of  increasing  angle  measured  with  respect  to  the  y 
axis.  For  this  case  the  red  shift  of  the  peaks  is  also 
attributed  to  the  increase  in  interaction  between  the 
spheroids,  because  for  an  isolated  spheroid  the  shift 
in  energy  as  a  function  of  angle  of  rotation  is  negligi¬ 
ble  for  slightly  eccentric  spheroids.  The  maximum 
enhancement  is  obtained  in  Figs.  2(a)  and  2(b)  when 
the  rotated  spheroid  is  at  an  angle  of  6  =  90°.  This 
occurrence  is  reasonable  because  the  semimqjor  axis 
of  the  rotated  spheroid  is  perpendicular  to  the  polar¬ 
ization  of  the  incident  electric  field  at  6  «  0°  and  is 
parallel  at  8  =  90°,  thereby  enhancing  the  field  by  its 
geometry  and  directing  the  energy  more  efficiently 
toward  the  other  radiating  spheroid.  As  we  ex¬ 
pected,  the  largest  red  shift  and  enhancement  occurs 
in  Fig.  3  where  both  spheroids  are  simultaneously 
rotated. 

The  resonant  wavelengths  for  the  cases  shown  in 
Figs.  1  and  3  have  been  studied  as  functions  of 
interparticle  distance  d  and  orientation  6,  respectively. 
Figures  4(a)  and  4(b)  display  the  behavior  of  the 
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Fig.2.  (a),(b)InUnaityoftbetaUielMtriclltUatth«poiiitPM* 
function  of  wavalength  for  aduster  of  i|4ieroida  witha/fr  >  0.9 and 
a  -  5  nm  for  diffemst  angleo  of  orionUtion  of  th«  •pharoids  with 
raapoet  to  the  vartieBi  ana.  Tha  diractiona  of  tha  inddant  radia¬ 
tion  and  the  polarixation  ara  aa  ahown. 


low-energy  peaks  as  a  function  of  orientation  ani^e 
and  of  interparticle  separation.  Figures  5(a)  and 
5(b)  provide  this  same  analysis  for  the  high-energy 
peak.  A  least-squares  fit  to  the  data  was  made  for  all 
four  cases  to  determine  the  behavior  of  the  red  shifts 
as  a  function  of  interparticle  separation  and  also  as  a 
function  of  orientation.  The  spread  of  calculated 
resonance  wavelengths  about  the  smooth  curves  is 
due  to  the  fact  that  all  resonant  wavelengths  analyzed 
in  Figs.  4  and  5  were  obtained  from  a  spline  inter^la- 
tion  of  the  calculated  duster  spec^  The  uncer¬ 
tainty  introduced  by  this  interpolation  in  determin¬ 
ing  the  resonant  wavelengths  is  of  the  order  of  ±2 
nm. 

In  the  case  of  both  the  low-  and  hi|d>*«n«rBy  P«aks 
for  the  rotated  spheroids,  the  shift  toward  the  red 
shows  a  sinusoidal  behavior,  as  is  illustrated  in  Figs. 
4(a)  and  5(a).  A  least-squares  fit  of  this  data  gives  a 
fiinctional  behavior  for  the  shift  that  is  governed  by 


Fig.  3.  Intensity  of  the  total  electric  field  at  the  point  P  as  a 
function  of  wavelength  for  a  cluster  of  simultaneously  rotated 
q>herotds  with  a/6  ■>  0.9,  a  »  5  nm  for  different  angles  of  rotation 
of  the  two  spheroids.  The  directions  of  the  incident  radiation  and 
the  polarisation  are  as  shown. 


the  simple  relationship  X.(0)  *  A  -f  B  sin^6):  For 
Fig.  4(a),  A  »  376.7  nm  and  B  *  47.8  nm;  for  Fig.  5(a), 
A  «  351.8  nm  and  B  «  7.2  nm. 

From  Figs.  4(b)  and  5(b),  it  appears  that  tbe 
position  of  ^e  peaks  have  an  eiqoonential  behavior  as 
a  function  of  increasing  interparticle  distance.  A 
curve  fitting  gives  the  functional  relationship  kid)  = 
k'  +  hk  ejqp{-(d  -  do)/‘n]-  For  the  case  of  the  low- 
energy  peak,  k'  was  t^en  to  be  the  resonant  wave¬ 
length  for  the  decoupled  cluster,  i.e.,  tbe  resonant 
frequency  of  an  isolated  spheroid.  Here  66  is  the 
difference  between  tbe  resonant  wavelength  of  the 
touching  spheroids  and  the  resonant  wavelength  with 
no  coupling  and  do  is  the  touching  distance  25.  The 


Fig.  4.  Shift  in  Uw  resonant  wavolongth  for  the  low-anargjr  paak 
as  a  fiinction  of  (a)  ths  roution  angis  6  of  tha  simultaiMoualy 
rotatad  spharoids,  and  (b)  intarpartiela  aaparaUon  d  for  tha  aligned 
spbaroidB.  Tha  solid  curve  raprasanU  a  laost-aquaras  fit  to  tha 
e^culatsd  points. 
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Fig.  5.  Shift  in  the  resonant  wavelength  for  the  high-energy  peak 
as  a  function  of  (a)  the  rotation  angle  6  of  the  aintultaneoualy 
rotated  spheroids,  and  (b)  interparticle  separation  d  for  the  aligned 
spheroids.  The  solid  curve  represents  a  least  squares  fit  to  the 
adeulated  points. 

only  fitting  parameter  used  is  tf,  which  corresponds  to 
a  decay  factor.  For  Fig.  4(h),  k  *  361.5  nm,  8X.  =  94 
nm,  and  the  fitted  paramet^  is  -q  «  0.57  nm.  For 
the  case  of  the  high-energy  peak  it  is  also  necessary  to 
fit  X';  the  resulting  fitted  curve  for  Fig.  5(b)  3rields 
X’  =  350  nm,  8X  =  13.6  nm,  and  t)  *  0.32  nm.  From 
the  analysis  it  is  evident  that  the  low-energy  peak 
suffers  a  considerably  larger  shift  in  wavelength  than 
the  hi|di*energy  peak  as  a  function  of  both  angie  and 
interparticle  distance. 

The  decoupling  distance  for  the  cases  discussed 
here  occur  at  of  »  36.  For  all  other  cases  studied  by 
us  the  decoupling  distance  has  been  found  to  be  d  ^ 
46,  which  is  twice  the  distance  of  touching  spheroids. 

4.  Conclusions 

A  T-matrix  formalism  has  been  used  to  calculate  the 
local  electric  field  near  clusters  of  metallic  spheroids, 
and  the  influence  of  clustering  on  their  spectrum  has 
been  investigated  as  a  function  of  interparticle  dis¬ 
tance  d  and  relative  orientation  6  in  the  long- 
wavelength  regime.  The  calculations  show  two  reso¬ 
nant  pedes  in  the  spectrum  caused  by  interparticle 
interactions,  which  is  in  contrast  with  one  peak  in  the 
spectrum  of  the  single  spheroid.  The  high-energy 
p^  disappears  and  the  low-energy  peak  becomes  the 
peak  of  the  isolated  spheroid  for  interparticle  dis¬ 
tances  of  46.  This  fact  should  prove  useful  in  the 
analysis  of  clusters  of  more  than  two  ellipsoidal 
particles. 

As  a  result  of  the  interaction  among  the  spheroids 
that  form  the  cluster,  enhancements  of  up  to  1(F  are 
observed.  Hi^er  enhancements  should  be  observed 
for  more  eccentric  ellipsoidal  constituents  of  the 
cluster.  The  enhancement  is  shown  to  be  sensitive 
to  the  relative  orientation  of  the  spheroids  and  the 
interparticle  distance.  Both  resonant  peaks  are  red 
shifldl  as  a  function  of  reduced  interpa^cle  distance 


and  of  increasing  an^  of  rotation  with  respect  to  the 
yaxis.  ^ 

The  position  of  the  hi^ier-energy  peak  is  less 
sensitive  to  interparticle  distance  than  the  lower- 
energy  one,  and  rotations  are  less  effective  in  shifting 
the  peaks  than  the  separation  of  the  spheroids. 
Through  a  least-squares  fit,  we  demonstrated  that 
the  shifting  of  the  low-energy  peak  obeys  an  exponen¬ 
tial  relationship  as  a  function  of  the  varying  interpar¬ 
ticle  distance  that  deca}rs  to  the  value  of  an  isolated 
spheroid  for  values  of  d  *  36,  whereas  the  red  shift 
ob^s  a  sinusoidal  relationship  as  a  function  of  the 
an^e  of  rotation  of  the  spheroids. 

For  experimental  systems  composed  of  a  collection 
of  scatterers  with  a  random  distribution  of  interparti¬ 
cle  distances  and  orientations,  these  calculations  will 
be  important  in  predicting  the  estimated  effective 
width  of  the  observed  resonances.  The  predicted 
width  is  expected  to  be  greater  than  the  ones  shown 
here  because  such  a  system  will  have  contributions 
from  clusters  with  a  wide  variety  of  interparticle 
distances  and  relative  orientations.  A  rou^  esti¬ 
mate  suggests  that  the  resonant  peak  of  that  s}rstem 
should  have  a  width  of  the  same  order  as  the  differ¬ 
ence  between  the  resonant  energy  of  the  cluster  of 
touching  spheroids  and  the  resonant  energy  of  a 
decoupled  duster. 

Recently  some  efforts  have  appeared  in  the  litera¬ 
ture  that  treat  clusters  of  spheroids  by  using  ellipsoi¬ 
dal  harmonics,*®  but  to  our  knowledge  only  calcula¬ 
tions  restricted  to  the  far  zone  have  ^n  performed. 
We  do  not  expect  that  the  use  of  ellipsoidal  harmonics 
will  simplify  in  any  significant  manner  the  calculation 
of  the  field  in  the  near  zone,  because  a  cluster  does  not 
have  either  spherical  or  ellipsoidal  symmetry;  there¬ 
fore,  we  do  not  anticipate  any  advantage  of  one 
expansion  over  the  other. 
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In  this  presentation  we  repon  results  obtained  by  the  first  combined  application  of  the  Metal  Oxide 
Semiconductor  and  the  Photocarrier  Grating  configurations.  This  combination  enabled  the  first 
simultaneous  study  of  the  two  carriers  mobfihy-iifetime  products  and  their  light  intensity  exponents  as 
a  funaion  of  the  position  of  the  Fermi  level,  in  undoned  a-Si:H.  We  found  that  anticorrelations  and 
correlations  prevail  between  these  two  sets  of  quantities.  The  conclusion  we  derive  from  these 
behaviois  is  Aat  the  "defea  pool"  model  accounts  for  the  phototranspon  data  much  bener  than  any 
other  model. 


1 .  INTRODUCTION 

One  of  the  methods  [1,2]  used  to  determine 
the  energetic  location  of  the  defea  states  in 
hydrogenated  amorphous  a-Si:H,  is 

finding  the  dependence  of  measureable 
phototranspon  properties  on  the  energy 
separation  between  the  conduction  band  edge 
in  the  bulk.  £^,  and  the  equilibrium  (or  "dark") 
Fermi  level,  Ep.  Predictions  for  the 
dependencies  of  the  two  carriers  mobility- 
lifenme,  p.T,  products  on  A£  *  •  Ep  have 

been  given  recently  by  vaiioSs  researchers  for 
both,  the  "standard"  [3,4]  and  the  "defea  pool" 
[4],  models  of  a-CI.II  Vruile  the  models 
proposed  were  for  undoped  a-Si:H.  the 
expenmental  results  were  limited  thus  fer  to 
data  obtained  on  doped  materials  [2.4,5]  in 
which  a  different  distribution  of  recombination 
centers  may  exist  [6,7].  In  the  present  work  we 
repon  the  fiist  determination  of  the  above 
dependencies,  in  undoped  a-Si:H,  applying  a 
combination  of  Metal-Oxide-Setniconduaor 
(MOS)  and  PhotoCarrier  Grating  (PCG) 
configurations.  Previously  either  the  MOS 
configuration  [8]  or  the  PCG  configuration  [9] 
have  been  used  for  the  study  of  phototranspon 
parameters.  The  novel  application  of  the 
combined  configuration  enables  a  simultaneous 
study  of  both  the  majority  and  minority  carrier 


phototranspon  properties  as  a  funaion  of 
AEj  =  Eg(s)-Ep,  where  Ej.(s)  is  the  cnerg>’  of 
the  conduction  band  edge  at  the 
semiconduaor-oxidc  interface.  In  panicular  we 
are  able  to  present  the  first  simultaneous  results 
of  the  dependencies  of  the  light  intensity 
exponents  on  AE^  and  to  show  that  these 
dependencies  on  yield  a  more  convincing 
proof  for  the  "defea  pool"  recombination- 
centers  model  than  "quantitative"  fits  of  the 
theoretical  predictions  to  the  measured  p.T 
products  [4].  Since  the  majority  carriers  in 
undoped  a-Si:H  are  the  electrons  we  define 
their  light  intensity  exponent  by  the  relation 
(u.T)g  K  where  is  the  electrons  p,T 
produa  and  G  is  the  earners  generation  rate. 
Correspondingly  for  the  holes,  we  define  the 
exponent  S  by  the  relation  (p,T)jj'-«  G'^.  The 
main  conclusion  of  the  comparison  of  our 
experimental  data  with  the  presently  available 
theoretical  results  is  that  the  "defea  pool"  is  the 
best  description  given  thus  far  for  the  deep 
recombination-level  distribution  in  a-Si:H. 
This  conclusion  applies  for  both  undoped  and 
slightly  doped  (a  fraction  of  a  ppm)  a-Si:H. 

2.  THEORETICAL  BACKGROUND 

Following  the  many  data  which  indicate  the 
presence  of  three  charged  states  of  a  dangling 
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bond  [10],  the  presently  accepted  models  of  the 
deep  states  (Le.  the  recombination  centers)  in 
a-Si:H  do  consider  these  states  [2*4]. 
Correspondingly,  the  mlKt  common  model 
used  for  interpreting  the  various  data  is  the 
simplest  possible  model  which  is  based  on  the 
presence  of  these  states.  In  this  model,  known 
as  the  “standard"  model  [2-4]  of  a-Si:H,  it  is 
assumed  that  there  are  D***  and  D'*  dangling 
bond  centers,  the  common  energy  level  of 
which  lies  below  Ep,  and  the  D'  centers,  the 
energy  level  of  which  lies  above  Ep.  The  most 
significant  feature  of  this  model  is  the 
conelated  occupation  statistics  [11]  of  these 
states,  i.e.,  that  the  occupations  of  the  D*^  and 
D°  centers  also  determine  the  occupation  of  the 
D‘  centers.  Very  recently  a  numerical  study  of 
the  "standard"  model  [3]  has  shown  that  in 
undoped  a-Si:H,  (M.T)g  will  inaease  with 
decreasing  AE,  while  (p.t)l  will  decrease  with 
AE.  It  was  filler  found  [3]  that  y  decreases 
from  1  to  1/2  with  decreasing  AE,  while  S  is 
independent  of  AE  having  the  value  S  *  0. 
Qualitatively  these  results  are  the  same  as  those 
expected  for  a  single  recombination  level  in 
general  [12.13]  and  for  parameters  which  are 
appropriate  for  a-Si;H  in  particular  [2].  It  is 
important  to  note  that  the  S=0  prediction  is 
independent  of  the  particular  parameters 
chosen  for  the  above  two  models  [13]. 

The  more  recently  suggested  model  for  a-Si:H 
is  the  so-called  "defea  poof  model  in  which 
one  assumes  three  independent  "sets"  of  levels, 
so  that  each  set  consists  of  a  D^/D*^  level  and  a 
D‘  level  [6].  Thus  far  prediaions  based  on  this 
model  were  made  [4]  for  the  AE  dependencies 
of  (|XT)g  and  CP'T)jj.  but  not  for  y  and  S.  We 
argue  however  here  that  while  y  will  have 
qualitatively  the  same  AE  dependence  as  in  the 
"standard"  model,  it  is  veiy  likely  that  in  the 
"defect  pool"  model  S  has  a  non-zero  value 
which  decreases  with  AE.  Basically,  the 
argument  is  based  on  the  very  simple  picmre 
[12]  of  the  single  recombination-center  level. 
In  this  picture  the  variation  of  y  with  increasing 
AE  is  controlled  by  the  charge  neutrality 
condition  which  is  fulfilled  by  the  faa  that  the 
electron  concentration  is  equal  to  the 
concentration  of  the  positively  charged  deep 
recombination  centers  [12]  centers  in  our 


case).  We  argue  [13],  by  symmetry,  that  the 
same  condition  holds  for  the  holes  and  the 
negatively  charged  recombination  centers  (D- 
centeis  in  our  case)  provided  (as  assumed  in  the 
"defea  pool"  model)  that  no  correlatinn  exists 
between  the  occupation  of  the  two  npes  of 
centers  [4].  Hence,  S  will  vary  in  a  way  similar 
to  y,  Le.  it  will  decrease  from  S  =  1/4  to  S  =  0 
with  decreasing  AE.  We  find  then  that  the 
S  ^  0  observatioEt  is  a  signiature  of  a  system 
which  has  at  least  two  occupation-independent 
levels.  A  more  detailed  discussion  and  the 
justification  for  considering  recombination  m 
the  deep  states  rather  than  in  the  band  tail 
states  will  be  given  elsewhere  [13]. 

3.  EXPERIMENTAL 

The  MOS  structure  used  in  the  present  stud\' 
was  almost  the  same  as  the  one  used  previously 
[8,14]  for  the  study  of  the  photoconductivity 
and  its  dependence  on  AEj  in  a-Si:H.  Since 
this  structure  has  been  described  and  discussed 
before,  its  description  will  not  be  repeated  here. 
Let  us  mention  only  that  our  l|xm-thick  fiJm  of 
device-quality  undoped  a-Si:H  has  been 
deposit^  [15]  using  rf  glow  discharge 
decomposition  of  silane.  The  new  expcnmental 
configuration  which  was  used  here  was  the 
application  of  the  PCG  technique  to  the  MOS 
structure.  We  should  point  out  that  the 
common  PCG  measurement  [9]  aaually  yields 
the  ambipolar  diffusion  length  L,  from  which 
one  derives  using  the  well  known  relation 
[15]  between  them  ((p.T)^  *  L^).  We  can 
determine  then  both  and  (m-t)jj  as  a 

funaion  of  AEj. 

4.  EXPERIMENTAL  RESULTS  AND 
DISCUSSION 

Since  the  relation  between  the  applied  gate 
voltage  ,Vq,  and  AE^  has  been  discussed  in  the 
past  [8]  we  can  scale  Vq,  which  is  in  the  range 
62Vq^-6V,  with  the  band  bending  scale  at  the 
surface,  A^,  which  is  in  the  range 
O.^AEj^.5eV.  We  note  that  it  is  essentially 
the  forward  bias  (Vq>0)  range  which  yields 
this  variatioii  [8,13].  An  independent 
experimental  proof  for  the  above  Vq-AEj 


scaling  is  derived  from  the  faa  that  the 
presently  found  variation  of  7  between  1  and 
OJ  (see  below)  with  the  above  variation  of  Vq 
is  in  excellent  agreement  widi  die  predicdom 
made  for  this  parameter  in  all  the  available 
theories  [1'3]  and  all  the  available  experimental 
results  which  either  used  doped  materials  [1^] 
(i.e.  the  7  dependence  on  A£)  or  the  MOS 
configuration  [8,14]  (Le.  the  7  dependence  on 
AEr).  Now  tlut  we  have  set  the  A£  -  AF 
scale  we  can  compare  our  results  with  me 
predictions  of  foe  "standard"  [3,4]  and  foe 
"defea  pool”  [4]  models. 


Sc-Sf  (eV) 

Figure  1.  The  measured  two  carriers  p,T 
products  as  a  function  of  AE^.  The  results 
are  compared  with  foe  predicdons  of  Rel  4  for 
the  "standard"  model  (dashed  curves)  and  foe 
"defea  pool"  model  (full  curve^. 


Our  experimental  results  for  foe  dependencies 
of  (p.T)  and  (p.T)jj  00  A£.  are  shown  in 
Figure.  L  For  comparison  with  the  theoretical 
prediaions  we  also  show  foe  calculated  results 
for  the  above  two  models  [4],  for  a  defea 
concentration  of  10^^  cm'^  It  is  seen  that  foe 
measured  dependence  of  (p.T).  on  is  in  a 
mudi  better  agreement  wim  me  "defea  pool" 
model  than  wim  the  "standard"  model 
prediaions.  By  examining  ofoer  "standard" 
model  predictions  [3]  it  appean  that  this 
conclusion  goes  be^nd  me  choice  of  me 
parameters  used  in  the  calculations,  since  me 
(p.T)  dependence  on  AE.  is  much  sharper  for 
the  "Kfea  pool"  model  On  the  omer  hand  we 
see  that  the  dependence  of  (p>t)|j  is  in  better 


qualitative  agreement  wim  the  "standard" 
model  These  two  behaviors  indicate  clearly 
that  me  true  situation  is  more  complicated 
me  "standard"  model  Le.  that  mere  is  at  least  a 
partial  occupation>independence  of  me 
recombination  levels  involved  (see  below).  We 
must  note  however  that  the  agreement  of  me 
(p.T)p/(|AT)jj^  value  wim  me  "defea  pool"  model 
premction  is  simply  a  consequence  of  me 
parameters  chosen,  as  can  be  seen  from  the 
predictions  for  me  "standard"  model  when 
different  parameters  have  been  used  [3]. 
Hence,  while  me  agreement  wim  me  "defea 
pool"  predictions  is  better,  in  view  of  me  lack  of 
parameter-effea-analysis  it  is  not  too 
convincing.  An  important  observation  is  me 
faa  that  the  present  results  are  almost  identical 
wim  me  results  obtained  by  variation  of  AE  by 
doping  [44,16].  This  has  mree  consequences. 
First,  it  lends  fimher  support  to  our  AF^»AF. 
scaling.  Second,  me  slight  doping  necessary  to 
move  Ep  from  E--0.8  to  E.-04  eV  docs  not 
yield  a  consicwrable  change  in  the 
recombination  level  struaure  in  comparison 
wim  me  undoped  material  and  third  (which  is 
very  important  for  me  present  experimental 
approach)  me  "vicinity  to  me  surfece"  nature  of 
the  combined  PCG'MOS  configuration  does 
not  effea  me  results. 
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Figure  Z  The  measured  light  intensity 
exponent  as  a  function  of  AE^. 

In  Figure.  2  we  show  me  AE^  dependencies  of 
me  light  intensity  exponents.  The  results  for  7 
in  Figtire.  2  are  in  excellent  agreement  wim 
mose  of  previous  experimental  [8,14]  and 


theoretical  [2,13]  works  on  undoped  a-Si:H, 
T^ere  7  was  found  to  decrease  quite  sharply 
from  7  *8  1  to  7  0.4  in  the  range  of  OJ  ^ 

AES  0.8  eV. 

So  for  aU  the  data  discussed  above  is 
consistent  qualitatively  with  both  the  one  level 
model  [2,12]  and  the  "standard"  model  [23]  of 
a-Si:H.  Let  us  turn  then  to  the  results  of  the 
S(AE  )  dependence  (reported  here  for  the  first 
time  tor  undoped  a-Si:H).  we  see  in 
Figure.  2,  S  decreases '  appreciably  with 
decreasing  AE^.  This  result  is  in  agreement 
with  the  experimental  finding  of  the  behavior  of 
S  when  boron-  doping  was  used  [16].  While 
these  and  the  present  experimental 
observations  of  the  behavior  of  S  are  consistent, 
the  variation  of  S,  and  its  high  non-zero  value, 
are  inconsistent  qualitatively  with  the 
expectations  implied  [13]  by  the  one  level  [12] 
or  the  "standard'‘[3]  model  (see  above).  Hence 
the  one-level  model  and  the  "standard"  model 
cannot  account  for  the  S  data,  and  consequently 
we  should  consider  the  next  model  in 
complexity  for  a-Si:H,  i.e.  the  "defea  pool" 
model.  This  model,  as  discussed  above  (Sec.  2), 
is  most  likely  to  yield  a  non-zero  value  for  S. 
Hence,  the  behavior  of  S  is  a  much  more 
convincing  tool  than  the  (ixt)^  and  (|i.T)jj 
results  [4]  for  the  determination  of  the 
recombination  level  structure  in  a-Si:H. 

5.  CONCLUSIONS 

We  have  followed  the  variation  of  the  four 
phototranspon  properties  (p.T)g,  (p-t)^^,  7  and 
S  as  a  funaion  of  the  energy  separation.  AE^  in 
undoped  a-Si:H.  We  found  that  the  qualitative 
SCAB  )  dependence  excludes  the  one  level 
model  or  the  "standard”  model  as  representative 
models  for  the  recombination  level  structure  in 
a-Si:H  while  it  does  not  exclude  the  validity  of 
the  "defea  pool”  model  This  conclusion  is 
further  supported  by  the  better  agreement  of 
the  AEj  dependencies  of  (p.T)g  and  (P'T)}^  with 
the  calculated  predictions  of  the  "defea  pool" 
model. 
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Study  of  the  density  of  states  in  a-Si:H  using  the  Si/electrolyte  ^stem 
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®Racah  Institute  of  Physics,  The  Hebrew  Univeisity,  Jerusalem  91904,  Israel. 
^Department  of  Physics,  University  of  Puerto  Rico,  Rio  Piedras,  PR  00931. 


Localized  states  in  a-Si:H  are  studied  by  pulsed  measurements  on  the  a-Si:H/electrolyte  (S.E) 
system.  The  S/E  inter&ce  is  essentially  blocking  to  current  flow  and.  as  a  result,  surface  space-charge 
layers,  ranging  from  large  depletion  to  very  strong  accumulation  conditions,  can  be  induced  and 
studied.  Measurements  in  the  depletion  range  under  illumination  yield  directly  the  total  densin  of 
occupied  states  in  the  entire  energy  gap.  This  is  useful  in  obtaining  a  quick  and  reliable  assessment  of 
the  quality  of  the  amorphous  films.  In  high-grade  films  we  find  that  the  total  densir>'  of  occupied  states 
is  around  10^®  cm'^  The  data  in  the^ccumulation  range,  on  the  other  hand,  provide  useful 
information  on  unoccupied  states  near  the  conduction  band  edge.  The  S/E  system  is  utilized  also  to 
apply  a  sweep-out  technique  for  an  accurate  determination  of  p,-,  the  product  of  the  electron  mobihrv 
and  lifetime,  even  when  this  value  is  very  low. 


1.  INTRODUenON 

Consideraole  effort  has  been  devoted  to 
derive  the  density  of  states  spectra  in  a-Si:H 
films  [1].  In  this  paper  use  is  made  for  this 
purpose  of  pulse  measurements  on  the 
semiconduaor/electrolyte  (S/E)  system.  Such 
measurements,  which  proved  to  be  very  useful 
in  the  study  of  crystalline  semiconduaors  [2], 
have  been  found  to  be  equally  effective  when 
applied  to  a-Si:H  films.  The  blocking  nature  of 
the  S/E  interface  allows  one  to  induce  by  an 
applied  bias  space-charge  layers  at  the  a-Si;H 
surface,  ranging  from  large  depletion  to  vciy 
strong  accumulation  conditions.  In  this 
manner,  the  entire  energy  gap  in  the  space 
charge  region,  together  with  its  localized  states, 
can  be  swung  below  and  above  the  Fermi  level. 
At  the  same  time,  one  can  measure  the  surface 
space-charge  density  as  a  function  of  the 
barrier  height  V^.  In  a-Si:H.  resides 
predominantly  in  the  localized  states  (except  in 
strong  accumulation  conditions),  so  that  the 
measurements  yield,  at  least  in  principle,  the 
energy  distribution  of  the  density  of  states. 

High-grade,  device-quality  films  are 
usually  close  to  intrinsic  (resistivity  —10'°  ohm- 
cm).  Such  high-resistivity  films  cannot  be 


handled  by  our  measurement  technique  because 
the  surface  space-charge  capacitance  cannot  be 
charged  within  the  short-duration  applied  pulse 
biases.  Accordingly,  most  of  the  results  to  be 
reported  here  have  been  obtained  under 
illumination,  for  which  the  photo  resistivm  is 
typically  10^  ohm-cm  or  less.  Our  results  in  the 
depletion  range  are  in  remarkably  good 
agreement  with  the  theoretical  curve  that  takes 
into  account  the  presence  of  localized  states. 
The  data  thus  yield  directly  the  total  density  of 
occupied  states  in  the  entire  energy  gap.  This  is 
very  useful  in  obtaining  a  quick  and  reliable 
assessment  of  the  quality  of  the  amorphous 
films.  In  high-grade  materials  we  find  that  the 
total  density  of  occupied  states  is  around  lO'** 
cm'^.  The  results  in  the  accumulation  range,  on 
the  other  hand,  provide  useful  information  on 
the  unoccupied  states  near  the  conduction-band 
edge.  Analysis  of  the  data  in  this  case,  however, 
is  more  difficult  because  surface  states  are 
apparently  also  involved. 

The  blocking  nature  of  the  a- 
Si:H/elecTrolyte  interface  is  utilized  also  to  apply 
a  sweep-out  technique  for  an  accurate 
determination  of  p,T,  the  product  of  the  electron 
mobility  and  lifetime.  The  technique  is  effective 
over  a  very  wide  range  of  p,T.  In  a  rather  poor- 
quality  film,  for  example,  we  have  measured  a 


V 


fiT  value  as  low  as  SxlO"®  cm^A^. 

2.  EXPERIMENTAL 

Device  quality  a-Si:H  Sims  were  prepared 
by  rf  glow-discharge  decomposition  of  silane. 
First,  a  thin  n+  layer  was  deposited  on  a 
conducting  glass  substrate,  follow^  by  a  1  pjn 
thick  intrinsic  film.  The  n+  layer  provides  an 
ohmic  contact  between  the  conduaing  glass 
and  the  intrinsic  film.  The  conduaing  glass  was 
cut  into  squares  of  about  0.5  cm^  in  area  and  a 
contact  attached  to  the  conducting  glass.  The 
unre  lead,  contaa  area  and  the  entire  sample 
were  masked  by  epoxy  cement,  except  for  a 
small  area  (—  2  mm^)  of  the  film’s  surface  to  be 
exposed  to  the  elearolyte.  The  sample  and  a 
platinum  elearode  were  immersed  in  an 
indifferent  elearolyte  such  as  Ca(NOj)2  or 
(NH4),S04.  The  sample  was  illuminated 
through  the  (transparent)  electrolyte  by  a  —2 
mW  He-Ne  laser.  The  measurement  technique 
has  been  described  elsewhere  [2],  and  will  be 
reviewed  only  briefly  here.  A  shon  (0.1-40 
pLsec)  voltage  pulse  applied  between  the  Pt  and 
the  sample  is  used  to  charge  up  the 
semiconduaor  space-charge  capacitance.  The 
voltage  drop  aaoss  these  elearodes,  measured 
just  after  the  termination  of  the  pulse, 
represents  to  a  very  good  approximation  the 
change  5Vj  in  barrier  height  aaoss  the  film’s 
space-charge  layer  induced  by  the  applied  pulse. 
TTie  change  SQ.^.  in  space-charge  density  is 
obtained  from  the  voltage  developed  aaoss  a 
series  capacitor,  again  at  the' termination  of  the 
pulse.  Pulses  of  varying  amplitude  are  applied 
singly,  one  per  data  point  taken.  In  this  manner 
damage  to  the  amorphous  film  is  minimized. 
We  found  that  applying  a  large  number  of 
pulses  degrades  the  material,  drastically 
reducing  its  resistivity  and  introducing  large 
trap  densities. 

The  quiescent  barrier  height  is 

determined  quite  accurately  (to  within  20  meV) 
from  measurements  in  the  depletion  range. 
The  entire  Q^j.  vs.  curve  can  then  be 
constructed,  on  the  basis  of  the  vs.  BVj 
data  [2].  In  what  follows  the  surface  electron 
density  Nj  »  where  q  is  the  electronic 


charge,  rather  than  Qj-  is  plotted  against 

The  dark  and  photo  resistance  of  each 
sample  (between  the  a-Si:H/elearolyte  interface 
and  the  n+  contaa)  was  derived  from  current- 
voltage  charaaeristics  measured  at  the  onset  of 
the  applied  pulse,  before  the  surface  space- 
charge  capacitance  can  be  charged  through  the 
sample’s  resistance  (see  below).  In  this  manner, 
the  blocking  interface  is  effectively  shorted. 


3.  RESULTS  AND  DISCUSSION 


Figure  1.  Surface  electron  density  Nj  vs.  barrier 
height  Vj  in  an  illuminated  sample  of  a-Si:H. 


Typical  results  of  vs.  V^.  obtained  for 
an  illuminated  a-Si^H  sample  from 
measurements  with  a  pulse  of  3-M.sec  (squares) 
and  30-p.sec  (circles)  duration,  are  displayed  in 
the  semilog  plot  of  Fig.  1.  In  the  depletion 
range  (Vj<0),  the  N.  values  are  negative  ana, 
because  of  the  logarithmic  scale  used,  the  plots 
are  those  of  their  absolute  magnitudes.  The  free 
electron  concentration  under  illumination  n|j| 
in  the  amorphous  film  has  been  derived  from 
the  measured  photoconductivity  on  the 
assumption  that  the  electron  mobility  p.  is  10 
cm^/Vsec  (3.4).  The  corresponding  theoretical 
dependence  of  Nj  on  V^,  as  obtained  from  a 
solution  of  Poisson’s  equation  in  the  absence  of 
localized  states,  is  shown  by  the  dashed  curve. 
The  experimental  points  are  seen  to  lie  well 
above  this  curve  indicating,  as  expeaed  for 
amorphous  films,  that  the  space-charge  layer  in 
both  the  depletion  and  accumulation  ranges  is 
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dominated  by  localized  states.  Note  that  is, 
as  expected,  nearly  zero  since  illumination  tends 
to  flatten  the  ban^  at  the  surface. 

When  a  depletion  layer  is  formed  (by  an 
applied  negative  pulse),  flee  and  trapped 
elearons  are  expelled  from  the  surface  region, 
leaving  behind  the  positively-cbarged  localized 
states  that  make  up  the  immobile  space-charge 
in  the  depletion  layer.  In  an  unilluminated 
sample  only  shallow  states,  down  to  0.3-0.4  eV 
below  the  conduction  band  edge,  are  able  to 
thermally  emit  their  trapped  electrons  into  the 
conduction  band  within  the  measurement  time 
(the  pulse  duration)  [2].  Under  illumination, 
on  the  other  hand,  the  situation  is  quite 
different.  Hole-electron  pairs  are  continuously 
generated  by  the  light  The  electrons—are 
expelled  from  the  surface  region,  while  the 
holes  are  attraaed  to  the  surface  where  thq^  can 
recombine  with  the  trapped  electrons.  In  this 
manner  elearons  in  occupied  states  throughout 
the  energy  gap  can  be  discharged  and  expelled 
from  the  space-charge  layer,  irrespective  of  the 
depth  of  the  states.  Referring  to  Fig.  1,  we  see 
that  for  a  measurement  time  of  3  psec  (squares) 
only  a  fraction  of  the  trapped  electrons  are 
expelled,  but  when  the  pulse  duration  is 
extended  to  30  psec  (circles),  praaically  all 
occupied  states  in  the  depletion  layer  are 
discharged  and  expelled  (see  below).  The  solid 
curves,  passing  in  each  case  thiough  the  points, 
were  calculated  from  Poisson’s  equation  for  the 
case  in  which  localized  states  of  densities  (as 
marked)  are  present.  The  agreement  between 
theory  and  experiment  is  seen  to  be  remarkably 
good,  adding  considerable  confidence  to  our 
analysis.  The  total  density  of  occupied  states  in 
the  entire  energy  gap,  as  derived  from  the  data 
for  the  30  psec  measurement  time  is  about 
10'®cm‘^,  which  is  the  expected  value  for  the 
films  studied. 

When  accumulation  layers  are  formed, 
unoccupied  states  are  filled  up.  Here  again,  the 
experimental  points  are  seen  to  be  initially  well 
above  the  dashed  curve  calculated  for  the  case 
of  no  localized  states  present,  indicating  that 
localized  states  dominate  the  space  charge  layer 
in  this  range  as  well.  Surfrce  states  are  probably 
also  involved.  At  stronger  accumulation. 


however,  the  calculated  curve  aosses  the 
experimeatal  points  and  climbs  above  them. 
This  is  not  understood  at  present,  and  more 
work  is  in  progress  to  account  for  such  a 
behavior.  Further  studies  arc  needed  also  to 
distinguish  between  surface  and  bulk  states. 


Figure  2.  Expelled  trap  densii\'  vs.  pulse 
duration  (measurement  time)  for  three  a-Si;H 
samples  under  illumination. 
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Figure  3.  Dark  and  photo  conductivitv'  of 
a  degraded  a-Si:H  sample  vs.  applied  voltage 
pulse.  Solid  curve  in  the  negative  pulse  range 
represents  eq.  (3)  for  pT  =  5  x  10’*  cm“/V. 


The  expelled  trap  density  under  depletion 
conditions,  as  derived  frxim  measurements  such 
as  those  shown  in  Fig.  1,  is  plotted  in  Fig.  2 
against  pulse  du^ non.  The  lower  two  curves 
were  obtained  fo:  ?h-grade  a-Si;H  films,  while 
the  upper  curve  -  .  obtained  after  a  film  has 
been  degraded  by  the  application  of  many 


voltage  pulses.  In  all  cases,  the  expelled  trap 
density  increases  with  pulse  duration,  but  tends 
to  saturate  at  a  pulse  duration  of  30-40  psec. 
We  interpret  the  saturation  level  as 
representing  the  total  density  of  occupied  states 
in  the  energy  gap.  We  are  in  the  process  of 
developing  a  model  involving  the  hole  trapping 
kinetics  in  order  to  account  for  the  shape  of  the 
experimental  curve. 

The  dark  and  photoconductivity  of  a 
degraded  a-Si:H  sample  vs.  applied  pulse 
voltage  is  displayed  in  Fig.  3  on  a  semilog  plot. 
The  measurements  were  taken  at  the  onset  of 
the  pulse,  before  the  space-charge 

capacitance  at  the  a-Si:H  interbce  can  charge 
up.  For  positive  pulses  (elearolyte  positive  with 
respect  to  the  n-^-  contact),  the  conductivi^  in 
both  the  dark  and  under  illumination  is  seen  to 
increase  with  applied  voltage.  This  is  due  to 
electron  injection  from  the  n+  layer.  In  the 
negative  polarity,  the  dark  conductivity  is 
independent  of  voltage,  as  it  should.  The 
photoconductivity,  on  the  other  hand,  is  seen  to 
decrease  with  pulse  amplitude,  approaching  the 
dark  conductivity  at  large  negative  voltages.  We 
attribute  this  behavior  to  electron  sweep-out  by 
the  applied  fields,  as  shown  by  the  following 
considerations.  The  continuity  equation  for  the 
case  of  a  blocking  contaa  can  be  written  as; 

dN/dt  =  L  -  N/7  -  p,EN/d ,  (1) 

where  N  is  the  total  density  of  photo  electrons 
in  the  sample  per  cm^  L  is  the  number  of  hole- 
electron  pairs  generated  by  the  light  per  cm^sec, 
E  is  the  applied  held  and  d  is  the  film’s 
thickness.  Strictly  speaking,  the  last  term  in  eq. 
(1)  should  be  p.En^](d),  where  n^](d)  is  the 
electron  concentration  at  the  n-t-  contact.  For 
simplicity,  we  have  replaced  n[ji(d)  by  the 
average  electron  concentration  N/d,  which  is  a 
fairly  good  approximation.  Space-charge 
effects  are  negligible  in  this  polarity,  even  if  all 
photo  electrons  are  swept  away  from  the 
sample  (leaving  behind  the  photogenerated 
holes).  Hence  the  field  is  very  nearly  uniform  at 
the  value  V/d.  Solution  of  eq.  (1),  with  the 
boundary  condition  that  N  assumes  the  zero- 
field  density  Nq  =  Lt  at  the  onset  of  the  pulse 
(t  *  0).  yields 


N=No/(l+ftTE/d)+ 

No«Pl-(l/T-i-p,E/d)t]/(l-i-d/M.TE).  (2) 

Now  d/p,E  is  the  electron  transit  time  through 
the  film  and  is  typically  lO"*  sec  for  an  applied 
voltage  of  1  V.  The  measurement  is  taken 
about  1  fjisec  following  the  pulse  onset,  so  that 
the  second  term  in  eq.  (2)  can  be  neglected  to  a 
very  good  approximation.  Re-writing  eq.  (2)  in 
terms  of  the  photoconduaivity  cr  =  qixN/d  we 
then  have 

cr  =  iTq/(1  -t-  p-xV/d^).  (3) 

The  solid  curve  fitted  to  the  photoconductivm. 
data  at  negative  voltages  in  Fig.  3  is  a  plot  of  eq. 
(3)  for  p.T  =  SxlO"*  cm^A'.  The  fit  with  the 
experimental  points  is  seen  to  be  ver>  good. 
The  rather  low  value  of  p-x  is  reasonable  for  the 
degraded  sample  used  in  the  measurement.  The 
method  can  be  applied  effectively  even  if  p.x  is 
an  order  of  magnitude  lower.  It  should  be 
noted  that  p.  in  this  context  is  the  trap- 
controlled  mobility  (4). 
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ABSTRACT 

A  comparative  study  of  the  deposition  temperature  (T,)  dependence  of  the  mobility-lifetime 
(/re)  products  of  the  charge  carriers  in  glow-discharge  and  rf  sputter-deposited  a-Si:H  is 
described  and  discussed.  The  T,-dependeooe  of  die  /ns  and  the  mqority  carrier  light-intensity 
exponents  of  the  two  types  of  films  are  strikingly  similar.  These  observations  then  lead  to  the 
amdusion  that  the  structure  of  the  recombination  levels  and  the  recombinatioD  processes  are  in 
accord  with  the  ‘‘defect  pool”  model,  in  contrast  to  previous  suggestions.  The  differences 
between  the  two  types  of  films  thus  ^pear  to  be  limited  to  the  concentrations  of  dangling 
bonds. 


INTRODUCTION 

Various  techniques  for  depositing  a-Si:H  have  been  studied  over  the  years,  some  of 
which  were  claimed  to  yield  hi^er  quality  films  than  others  [1].  The  glow  discharge  (GD) 
decomposition  of  silane  has  been  exhaustively  investigated,  and  most  of  the  applications 
employ  this  deposition  method  [2].  The  study  of  films  deposited  by  other  techniques  was 
usually  limited  to  the  description  of  their  measurable  /noperties  and  the  comparision  of  these 
properties  with  those  of  undoped  “device-quality”  GD  films  [3].  All  of  the  theoretical 
progress,  devoted  to  the  understanding  of  the  electronic  structure  of  a-Si:H,  was  achieved  by 
comparing  the  theoretical  predictions  with  the  experimental  results  on  GD  films  [4-9].  The 
question  of  whether  the  (^r^transpoit  data  on  material  produced  differently  indicate  different 
electronic  structures,  transport  mechanisms,  ce  recomUnahon  processes  has  been  addressed 
only  rarely.  One  suc^  stutty  of  the  phototransport  data  associate  with  the  majority  carriers  in 
GD  and  rf  sputtered  (RFS)  was  reported  by  Jousse  et  al.  [10],  who  monitor^  the  effects  of 
boron  doping  on  the  photoutmsport  parameters.  They  concluded  that  in  RFS  materials  the 
recombination  process  is  control!^  only  by  the  valence  bandtail,  but  in  GD  films  the 
recombinatioo  involves  a  defect  center  with  either  a  single  or  two  positively  correlated  energy 
states.  This  paper  reexamines  these  condusions  in  view  of  new  phototransport  data  and  tte 
recent  theoretit^  develc^nnents  [6,9]  in  understanding  the  behavior  of  the  minority  carrier 
phototransport  properties  in  OD  device-quality  material  [6-9].  RFS  films  were  course 
extensively  shidied  [11],  and  in  contrast  to  GD  deposition,  enable  a  T, -independent  control  of 
hydrogen  incofpmation. 

The  results  on  two  sets  of  sam/des,  a  OD  and  an  RFS  deposited  set  are  described  in  this 
pepCT.  The  variable  deposition  parameter  used  to  obtain  sequential  films  with  continuously 
varjnng  properties  in  bofli  sets  was  Ts.  Indeed,  this  parameter  has  been  employed  extensively 
for  numerous  purposes,  e.g.,  in  finding  correlations  between  various  structural  [13]  and/or 
electronic  pre^wrties  P].  This  parameter  was  conspicuous  in  estaUishing  the  relation  between 
the  deep  stales  density  and  the  width  of  die  valence  bandtail  in  OD  a-Si:H  [1^.  In  addition,  it 
is  probably  the  most  convenient  variable  depoution  parameter,  and  its  effect  on  the  films  is 
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coDoeiMiially  very  dear  and  uDCODtested:  A  higher  T,  enhances  the ‘‘order**  and  density  of  the 
distorted  tetruhedial  network  on  the  one  hand  [13]  and  the  effraon  oTlhe  hydrogen  on  the 
other  hand  [14],  Since  the  fonner  process  improves  the  transport  and  the  pesavadoo  of  the 
dangling  bends  wlule  die  latter  increases  their  density,  it  is  not  surprising  diat  the  charge  carrier 
lifetimes  peak  around  T, »  2S0*C  However,  the  de&uls  of  the  dynamical  processes  and  dieir 
dependence  on  the  substrate  temperature  are  obviously  com^ex.  For  example,  the  relative 
roles  of  T,  in  hydrogen  diffusion  [17]  and  incorporatian  [IQ  ate  still  unsettled.  Similarly,  its 
direct  effect  on  the  network  rdaxadon  and  indirect  effect  on  the  relaxation  via  hydrogen 
effusion  [19]  are  also  cot  well  undmatood.  This  paper,  wfaidi  is  devoted  .to  phototransport 
data,  mal^  no  attempt  to  correlate  the  data  with  the  stmttuial  properties  of  the  materials  since 
the  transport  properties  of  semiconductors  ate  very  sensitive  to  sm^  variations  in  the  structural 
and  compositioiial  details  of  dre  material.  Hence,  a  very  comprehensive  study  will  be  needed 
in  order  to  discuss  that  conelatian.  However,  the  basic  rratuie  of  the  recombiri^an  mechanism 
does  not  have  to  depend  on  those  details.  This  basic  mechanism  is  the  subject  of  this  paper. 


EXPERIMENTAL  PROCEDURE 

The  GD  films  that  were  dqxjsited  at  Solarex  [21]  at  250*C  and  at  a  typical  dqxxitiQn  rate 
of  4  A/s  onto  Coming  7QS9  glass  substrates  were  <^5000  A  thick.  The  RI^  nims  prepared  at 
the  Ames  Laboratory  [22]  were  deposited  by  sputtering  a  6”  diameter  polyciystalline  Si  target 
located -1.5**  above  tltt  heated  substrate  at  13  j6  MHz  and  an  rf  power  of  400  W.  Thepartial 
pressures  of  the  Ar  and  H2  gases  were  10  and  1  mt.  resp.  The  thickness  and  deposition  rates 
of  the  films  ranged  from  0.8  to  1.5  }ax\.  and  1.5  to  2.9  A/s,  resp.  The  contacts  were  coplanar 
KiCr  or  Au  contacts  separated  by  0.4  mm.  While  die  optoelectronic  properties  of  GD  films 
barely  change  upon  armealing  at  2S0*C  [16],  those  of  the  RFS  samples,  which  were  annealed 
at  that  temperamre  in  air  for  2  h,  indeed  improved. 

The  phototransport  properties  were  studied  by  the  photocarrier  grating  (PCG)  technique 
[23-25]  using  a  H>Ne  laser,  i.e.,  with  iiluminadoD  that  can  be  considered  to  be  Sinifoimly" 
absorbed  in  the  films.  The  results  were  obtained  under  ambipolar  conditions  [24]  so  the 

minority  carrier  (hole)  mobility-lifetime  product  could  be  derived  directly  from  the 
measured  ambipolar  diffusiao  length  L  via  ^  reladcm 

=  (q/2kDL2  (1) 

The  majority  carrier  (electron)  (/rc)e  product  was  derived  dirccUy  from  the  measured  photocm- 
ductivity  Opj,  and  the  rdadon 

(Kc)e  =  o^(qG)  (2) 

where  O,  the  fdiotogeneradai  rate  of  electron-hde  pairs,  was  determined  by  ajMuming  unit 
quantum  efficiency  ^  generation  by  the  absorbed  photems.  As  mentioned  above,  however,  the 
most  significant  parameters  for  interpretation  of  the  data  are  the  so-called  light  intensity 
oqponents  defined  by  die  relations 


(>«),  «  0^^ 

(3) 

(Ki)h  "  0-“- 

(4) 

The  f4x  resi’llB  shown  below  were  obtained  at  10  mW/cm^  illuminarion,  and  the  exponents 
were  derived  from  the  behavior  in  the  0.8  -  10  tnW/cm^  range.  Theanalysisof  these  data  was 
described  in  previous  rqiorts  [203]- 
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conceptually  very  clear  and  uncontested:  A  higher  T,  enhances  the  “oixier*’  and  density  of  the 
distorted  tetrahedral  network  on  the  one  hand  [13]  and  the  effusion  of  the  hydrogen  on  the 
other  hand  [14].  Since  the  fonner  process  imixdves  the  transport  and  the  passivation  of  the 
dan^ng  bonds  while  die  latter  increases  their  density,  it  is  not  surprising  that  the  charge  carrier 
lifetimes  peak  around  T,  %  250*C  However,  the  details  of  the  dynamical  processes  and  their 
dependence  on  the  substrate  temperature  are  obviously  complex.  For  example,  the  relative 
roles  ofT,  in  hydrogen  diffusion  [17]  and  incorporation  [18]  are  still  unsettled.  Similarly,  its 
direct  effect  on  the  network  relaxation  and  indirect  effect  on  the  relaxation  via  hydrogen 
effusion  [19]  are  also  not  well  understood.  This  paper,  which  is  devoted  to  photoiransport 
data,  mates  no  attempt  to  correlate  the  data  with  the  stractural  properties  of  the  materials  since 
the  uansport  properties  of  semiconductors  are  very  sensitive  to  sm^  variations  in  the  structural 
and  compositioiial  details  of  the  material.  Hence,  a  very  comprehensive  study  will  be  needed 
in  order  to  discuss  that  correlation.  However,  the  basic  nature  of  the  recombination  mechanism 
does  not  have  to  depend  on  Chose  details.  This  basic  medianism  is  the  subject  of  this  paper. 


EXPERIMENTAL  PROCEDURE 

The  GD  films  that  were  deposited  at  Solarex  [21]  at  250*C  and  at  a  typcal  deposition  rate 
of  4  A/s  onto  Coming  7059  glass  substrates  were  ~5000  A  thick.  The  RFS  films  prepared  at 
the  Ames  Laboratory  [22]  were  deposited  by  sputtering  a  6"  diameter  polycrystalline  Si  target 
located  -1.5”  above  the  heated  substtate  at  13  MHz  and  an  if  power  of  400  W.  The  partial 
pressures  of  the  Ar  and  H2  gases  were  10  and  1  mL  rest).  The  thickness  and  dep>osition  rates 
of  the  films  ranged  from  0.8  to  1.5  ;«n  and  13  to  2.9  A/s,  rcsp.  The  contacts  were  coplanar 
NiCr  or  Au  contacts  separated  by  0.4  mm.  While  the  optoelectronic  properties  of  GD  films 
barely  change  upon  aiuiealing  at  250*C  [16],  those  of  the  RFS  samples,  which  were  annealed 
at  that  temperature  in  air  for  2  h,  indeed  improved. 

The  phototransport  properties  were  studied  by  the  photocarrier  grating  (PCG)  technique 
[23-25]  using  a  Hc-Ne  laser,  i.e.,  with  illumin^on  that  can  be  considered  to  be  ‘^formly” 
absorbed  in  the  films.  The  results  were  obtained  under  ambipolar  conditions  [24]  so  the 

minority  carrier  (hole)  mobility-lifetime  product  {}n\  could  be  derived  directly  from  the 
measured  ambipolar  diffusion  length  L  via  the  relation 

(Kt)h  =  (q/2kr)L2  (1) 

The  majority  carrier  (electron)  (/n:)^  product  was  derived  directly  from  the  measured  photocon¬ 
ductivity  Op),  and  the  relation 

=  V^qG)  (2) 

where  G,  the  photogeneration  rate  of  electron-hole  pairs,  was  determined  by  assuming  unit 
quantum  efficiency  of  generation  by  the  absorbed  photons.  As  mentioned  above,  however,  the 
most  significant  parameters  for  interpretation  of  the  data  are  the  so-called  light  intensity 
exponents  [26],  defined  by  the  relatians 

On),  «  QT-i  (3) 

(;«5h  “  (4) 

The  }4X  results  shown  below  were  obtained  at  10  mW/cm^  illumination,  and  the  exponents 

were  derived  from  the  behavior  in  the  0.8  - 10  mW/cin^  range.  The  analysis  of  these  was 
described  in  previous  reports  [20,23]. 
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RESULTS 

Figure  1(a)  displays  the  well-known  behavior  of  (/rc)^  in  GD  a-Si:H  [10,14].  It 
increases  significantly  with  increasing  T,  up  to  «2(X}*C,  and  decreases  considerably  at  Tg  2 

300*C.  The  behavior  of  (jn\  is  similar  (Hgure  1(b)),  but  it  decreases  at  T,  ^  350*C.  The 

behavior  of  y  (Figure  1(c))  appears  to  be  correlated  with  (fn)^  in  some  ranges  of  T,  but  anti¬ 
correlated  to  it  in  other  ranges  of  the  deposition  temperature.  On  the  other  hand,  the  behavior 

of  S  (Figure  1(d))  aiipears  to  be  anticoirelated  with  over  the  entire  range  of 

Due  to  their  poor  photoconductivity.  oi  the  as-deposited  RFS  films  could  not  be 

measured  by  the  PCG  technique.  Although  Opi,  improved  after  annealing,  it  was  still 
insufficient  to  enaUe  the  deteimination  S  (£q.  (4)),  which  required  attenuated  illumination. 

Figure  2(a)  displays  the  dependence  of  (jn)^  aa  T,;  it  is  strikingly  similar  to  the  dependence  in 
GD  films  (Fig.  1(a)).  On  the  other  hand,  the  peak  of  (}n\  vs  T,  is  much  narrower  peak  in  the 

RFS  samples  than  in  the  GD  films.  The  behavior  of  y  (Figure  2(c))  is  interesting.  A  priori  the 
behavior  of  the  annealed  films  appears  to  be  anticorrdaled  with  that  of  the  unannealed. 

However,  comparison  with  -YCr,)  of  the  GD  films  (Rg.  1(c))  indicates  that  all  of  the  T,- 
dependences  can  be  matched  by  s^piopriate  displacements  along  the  T,  axis.  This  approach  is 
adopted  in  analysing  and  interpreting  the  results. 

DISCUSSION  AND  CONCLUSIONS 

As  mentioned  above,  the  general  qualitative  dependence  of  on  T,  is  well  understood: 
At  low  Tg,  two  factors  are  responsible  for  the  high  dbfect  concentration  and  consequent  short 
lifetime  of  both  types  of  charge  carriers:  (i)  H  dUTusion  during  deposition  is  too  slow  to 
piassivate  many  of  the  dangling  bonds  and  (ii)  the  bandtails  are  broad,  due  to  the  unrelaxed 
network,  which  yields  shallower  recombination  centers  [14].  At  high  T,,  hdyrogen  effusion 
increases  the  density  of  dangling  bonds  [14,15]  and  may  also  lower  the  network  order  [19] 
which  is  improved  1^  its  inoorporatiai.  These  effects  are  confirmed  by  various  measurements 
[14],  notably  the  constant  photocurrent  mediod  (CFM)  [3].  Hence,  the  results  on  the  GD  films 
shown  in  Figs.  1(a)  and  1(b)  can  be  understc^  simply  within  this  picture.  The  detailed 

behavior  of  (/rr),  at  18O  s  Tg  ^  300'C  can  be  tftributed  to  the  competition  between  the 

competing  processes.  This  may  also  explain  the  variations  in  the  (jn)^  behavior  in  films 
deposited  under  different  conditions  [3]. 

While  die  jn  behavior  cannot  disclose  more  specific  infoimation  on  the  basic  recombina¬ 
tion  medianism,  the  light  intensity  exponents  y  and  S  (as  shown  for  photooonductcxs  in  gener¬ 
al  [26]  and  for  a-Si:H  in  particular  [9,20])  can.  Turning  to  the  former,  we  see  that  y  i.  0.9. 

‘The  Rose  model”  o(  bandtail  recombinatioo  [26]  does  not  allow  y  >  1.  and  suggests  that  the 
diaracteiistic  width  of  the  eqxxieDtia]  conduction  bandtail  is  given  by 

Eep-  [y/(1-Y)3W.  (5) 

If  recombination  occurs  only  in  the  bandtail  states  1  ^  g  e  0.9  implies  that  ^  2  0.2  eV.  This 
width  is  far  larger  than  any  previous,  more  direct  estimate  of  (»  25  meV)-  Indeed, 
apfdication  of  the  ”Ro8e  model”  [26]  for  a-Si:H  led  to  the  conclusion  that  the  recombination 
occurs  not  in  a  bandtail  but  in  a  de^  localized  state  [27].  Another  iniexpretahon  of  the  the  Rose 
model  was  suggested  for  the  reoombinaiioD  domiDating  the  minority  earners  (holes)  [2^.  In 
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Figure  1.  The  d^xtsUion  tempemure  T,  dependence  of  the  photatransport  parameters  of  the 
GDftbns:  The  mt^lttty^lifetime  pxproducs  of  the  electrons  (a)  and  holes  (b),  and  the  corres¬ 
ponding  dependences  of  the  Ught  intensity  exponents  y(c)  and  S(d)  (see  eqs:  (3)  and  (4)). 


Iliis  caoe  the  model  yidds 

[(l-2SV2S]kT.  (d) 

Unlike  the  cue  of  Ebo.  th*  vdueo  ibown  in  Rg.  1(d)  ere  not  in  iharp  oontmt  with  our  know¬ 
ledge  of  E^o  P]-  Fbr  example,  for  S  *  0.15  we  get  E,o  =  r7/3)kT  at  room  tempemture. 
However,  thia  modd  is  still  inappropriate  since  S  decreases  (i.e.,  increases)  as  T 
mcieases  from  200  to  300*C,  while  CPM  data  clearly  jwdiMty  a  nanowing  vdenoe 

TheincreasiagEyo  also  sharply  contrestswlA  the  above  picture  of  enhanced  network 
Older  at  higher  T,. 
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The  foregoing  discussion  leads  to  an  interpretation  in  teims  of  deep  recombination  centers 
which,  in  a-SirH,  must  consist  of  various  dangling  bonds.  The  simple  single  recombination 

level  model  is  easily  ruled  out  in  the  range  y  >  1.  since  these  values  require  at  least  two  levds, 

the  occupation  of  which  is  a  primi  uncorrelated  [26}.  For ’f<  I  the  results  shown  in  Fig.  1(c) 
can  inde^  be  interpreted  as  due  to  a  single  level.  However.  Fig.  1(d)  rules  this  out  since  a 
single  level  system  was  shown  to  yield  S  =  0  [9].  For  a*Si:H  it  is  reasonable  to  consider  the 
correlated  occupancy  or  “staDdard”  model  of  dangling  bonds  [5.7,10]  which  assumes  a  singly 
(D^  or  unoccupied  (D^)  level  lying  below  the  Fenni  level  Ep  and  a  doubly  occupied  (D*)  level 
above  Ep.  The  correlated  occupancy  of  these  levels  was  recently  argued  to  yield  the  same 
behavior  as  a  single  recombinatioD  level  [9].  Indeed,  the  ‘'standard”  model  [^  yields  S  =  0  for 
this  picture  as  well.  In  summary,  while  Figs.  1(a)  >  1(c)  are  consistent  with  the  “standard” 
model,  Fig.  1(d)  is  not.  Thus,  another  independenlly  occupied  level  must  be  added  to  account 
for  all  of  the  results.  One  possibility  is  to  add  such  a  level  to  the  “standard”  model.  Calcula¬ 
tions  which  essentially  consider  the  valence  b^dtail  states  as  another  recombination  center 
indeed  yielded  posiuve  values  of  S  in  accordant  with  Fig.  1(d)  [7,8].  Yet  these  calculations 

yield  S  >  0  only  when  y  >  ^  [8].  in  contrast  to  the  data  shown  in  Fig.  1(c)  in  the  most 
interesting  range  of  T^  namely  200  £  T,  ^  300*C 

In  contrast  to  the  above  models,  the  recent  “defect  pool  ”  model  of  essentially  uncorrelated 
dangling  bond  bands,  where  a  D*  band  lies  below  a  band  which  lies  below  a  band  [5],  is 
consistent  with  all  of  the  data  shown  in  Figs.  1(a)  -  1(d).  The  uncorrelated  nature  of  the  bands 

enable  S*0  values,  while  their  number  and  energetic  order  relieve  the  correlation  between  y  > 

1  and  S  >  0.  Indeed,  for  the  light  intensities  used  in  this  work  the  values  of  y  =  0.95  and  S  = 

0. 18  arc  to  be  expected  from  this  model  [6].  Therefore,  the  S  >  0  values  exclude  the  earlier 
model  [10]  for  GD  material  and  show  that  the  the  “defect  pool”  model  is  to  date  the  simplest 
which  can  account  for  all  of  the  phototransport  data. 

The  annealed  RFS  films  (Fig.  2)  exhibit  close  similaribes  with  the  GD  samples  (Fig.  1). 
In  particular,  the  two-peak  behavior  of  Figs.  1(a)  and  2(a)  indicates  a  similar  electronic 
structure  and  recombination  mechanism  for  electrons  in  the  two  types  of  a-Si:H.  It  is  futher 

interesting  that  the  as-deposited  values  of  the  low-T,  films  are  much  lower  than  the  values 

of  the  aimealed  or  GD  samples.  On  the  other  hand,  the  Qex)^  values  of  all  the  higb-T,  samjdes 
are  close.  This  implies  that  for  sufficiently  high  T,,  but  lower  than  that  yielding  poor  photo¬ 
transport  properties,  the  H  content  and  network  disorder  are  independent  of  the  initial  deposi¬ 
tion  process. 

The  behavior  of  (;n:)),  in  the  RFS  films  is  similar  to  that  in  the  GD  films  around  the  p^ 
at  T,  =  250*C.  but  that  pe^  is  much  narrower  in  the  RFS  samples  (Figs.  1(b)  and  2(b)).  If 

(jrc)^  is  controlled  mainly  by  the  D*  centers,  then  the  results  indicate  that  their  content  is  more 
dependent  on  the  deposition  method  than  the  density  of  the  electron  traps,  which  are 

presumably  the  D9  and  D^  centers.  The  similar  behavimr  of  [}a)^  suggests  diat  the  y  values  of 
the  annealed  RFS  films  (Fig.  2(c))  are  essentially  the  same  as  those  of  the  GD  films  (Fig.  1(c)) 

if  the  curves  are  shifted  to  match  the  y  >  1  regions  in  the  two  types  of  samples.  They  >  1  value 
of  the  RFS  film  at  T,  s  150*C  again  rules  out  the  baodtail-only  [10],  the  single  recombination 

level,  or  the  “standard”  model  without  bandtails.  For  y  <  1  the  data  are  inconsistent  with  the 
*7>andtail  only”  model  (see  above)  but  consistent  with  all  of  the  other  models  in  which  the 

approach  of  Ep  to  the  dominant  recombination  level  (say.  the  D°)  will  decrease  (ja)^  and 
increase  (/it)i,  and  y.  However,  the  striking  similarity  of  the  T, -dependence  of  and  y  in 
the  annealed  RFS  and  GD  films  strongly  suggests  that  the  “defect  pool" model  is  ai^icable  to 

the  aimealed  RFS  films  as  well.  In  particular,  the  decreasing  y(T,)  from  y  >  1  and  its  detailed 
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behavior  provide  strong  evidence  for  the  similaiity  of  the  reocnni^nation  processes  in  the  two 
Qrpes  of  naterials. 

In  conclusion,  the  T, •dependence  of  die  pbototransport  parameters  of  glow-discharge 
(GD)  a-Si;H  indicated  the  presence  of  at  least  two  significant  occupation-independent 
recombination  levels.  The  results  are  omsistent  with  the  “defect  pool"  model.  The  closely 
similar  behaviw  of  the  majority  carrier  phototxanspoit  parameters  of  GD  and  annealed  rf 
sputter-deposited  (RFS)  lllms  suggest  similar  recombinadon  processes.  The  differences  are 
apparently  due  to  high^  dangling  bond  contents  as  well  as  different  relative  densities  of  the 
various  types  of  dangling  bonds  in  the  RFS  films. 
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EXCITING  FREQUENCY-DEPENDENT  RAMAN  SCATTERING 
IN  a-Si,.,C^H  ALLOYS*,  Gerardo  Morell.  R.S.  Katiyar,  S.Z.Weisz, 
Deparment  of  Physics  ,  University  of  Puerto  Rico, San  Juan,  PR, 
00931,  USA;  and  Isaac  Balberg,  The  Racah  Institute  of  Physics, 
The  Hebrew  University,  Jerusalem  91904,  Israel 


We  have  measured  the  Raman  spectra  of  a-Si,.,C,:H  alloys  using 
three  different  frequencies  of  energies;  2.71, 2.41,  and  2.14  eV.  These 
alloys  are  characterized  by  their  Tauc  optical  gap,  that  falls  between 
1.8  dnd  2.4  eV.  As  more  carbonjs  incorporated  into  the  a-Si:H  matrix 
the  width  of  the  Si-Si  TO  band  as  well  as  the  TA/TO  band  intensities 
ratio  decrease  indicating  an  increased  order  in  the  network.  However, 
the  degree  of  ordering  in  all  samples  is  inhomogeneous.  The  surface 
layer  is  more  disordered  than  the  bulk,  as  can  be  concluded  from  the 
changes  in  the  spectra  of  any  given  sample  as  a  function  of  exciting 
energy.  Spectra  taken  with  the  ll4  eV  line  shows  a  systematically 
higher  TO  mode  frequencj-  of  about  5  cm  '  and  a  smaller  width  of 
about  10  cm  '  than  the  ones  measured  with  the  2.41  or  2.71  eV  lines. 
Hydrogen  dilution  during  the  growth  process  is  seen  to  have  a  limited 
effect  on  improving  the  degree  of  ordering  in  the  alloys.  Also,  an 
intensity  enhancement  is  observed  in  some  samples  that  can  be 
explained  as  a  resonance  effect. 


•Work  supported  in  pan  by  EPSCoR-NSF  gram  EHR-9108775  and  by 
U.S.  Army  Research  Office  gram  No.  DAALO3-8y-G-01 14. 
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Infrared  and  Raman  Studies  on  a*Ge,.,Sn,:H  Thin  Films 

E.  Ching-Prado,  RS.  Katiyar,  W.  Muiloz,  0.  Resto,  and  S.Z.Weisz 
Department  of  Physics,  University  of  Puerto  Rico 
Rio  Piedras,  P.R  00931  USA 


Abstract 

Infrared  and  Raman  measurements  were  carried  out  on  a-Ge,.xSn^;H  thin  films  grown  on 
Si  substrates.  The  hydrogen  concentration  is  found  to  decrease  with  increasing  tin  content. 
Preferential  formation  of  monohyjryde  groups  (GeH)  is  observed  with  increasing  Sn 
concentration.  Tin  dependence  of  the  GeH  and  GeHj  stretching  modes  are  explained  due  to  the 
lower  electronegativity  of  Sn  in  comparison  to  that  of  Ge.  The  study  reveals  that  the  stability- 
ratio  electronegativity  of  germanium  must  be  revised  or  its  value  should  be  smaller  than  3.59  for 
the  a-Ge:H  system.  Also,  the  GeH  wagging  and  GeH-  roll  modes  show  that  there  is  a  large  bond 
angle  variation  with  increasing  Sn  concentration.  Additionally,  the  changes  in  the  frequency, 
width  and  intensity  of  the  Ge-Ge(TO)  like  phonon,  clearly  indicate  that  the  structural  disorder 
increases  with  increasing  tin  content.  Auger  electron  microprobe  and  Raman  study  with  different 
excitation  laser  lines  show  differences  near  the  surface  in  comparison  to  the  bulk,  which  is 
produced  by  the  partial  segregation  of  tin  atoms,  formation  of  Sn-clusters  in  the  surface,  and 
reduction  of  the  germanium  concentration  near  the  surface,  as  well  as  preferential  anachment  of 
hydrogen  to  Ge  than  to  Sn. 


